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ABSTRACT 


The  problem  of  electromagnetic  wave  propagation  along  a  dielectric 
rod  of  elliptical  cross  section  Is  considered.  The  field  components  and 
the  dispersion  relations  of  the  principal  modes  are  obtained.  The  prin¬ 
cipal  modes  degenerate  to  the  well  known  HE  modes  of  the  clrculeu: 

mn 

dielectric  rod  as  the  eccentricity  of  the  elliptical  rod  approaches 
zero.  It  Is  foiind  that  there  are  two  non- degenerate  principal  modes 
which  possess  no  cut  off  frequencies.  They  are  called  the  dominant 
principal  modes. 

In  contrast  to  the  case  of  a  circular  dielectric  rod,  the  boundary 
conditions  for  the  elliptical  rod  cannot  be  satisfied  by  using  a  single 
product  term  consisting  of  a  radial  and  a  periodic  Mathieu  function  of 
a  specific  order  to  describe  the  field  conponents  in  both  regions  (the 
region  inside  the  rod  auid  the  region  outside  the  rod) .  It  is  generally 
believed  that  sui  infinite  series  of  such  product  terms  must  be  used  to 
describe  the  field  components  in  both  regions.  In  the  present  investi¬ 
gation,  it  is  shown  that  the  boundary  conditions  may  be  ftilfilled  if  the 
field  components  in  one  of  the  two  regions  are  represented  by  a  single 
product  term  consisting  of  a  radial  and  a  periodic  Mathieu  function  of 
a  specific  order.  The  field  components  in  the  other  region  are  then 
represented  by  an  infinite  series  of  such  product  terms.  The  problem 
is  therefore  sufficiently  sinplifled  to  permit  analysis. 

The  propagation  characteristics  (the  propagation  constant,  the  field 
distribution  and  the  attenuation  constant)  of  the  dominant  principal 
modes  are  given  theoretically  and  experimentally.  It  is  found  that  the 
analytic  and  e3q>erlmental  results  are  in  very  good  agreement.  The  ft's 
of  a  dielectric  rod  cavity  resonator  supporting  the  dominemt  principal 
modes  are  also  given. 
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CHAPTER  I 


INTRODUCTION 


The  concept  of  guiding  electromagnetic  waves  either  along  a 
single  conducting  wire  with  finite  surface  impedance  or  along  a  dielec¬ 
tric  rod  is  not  new.  As  early  as  lQ99}  Sommerfeld  (l)  conceived  the 
idea  of  guiding  a  clrciilarly  symmetric  TO4  wave  along  a  conducting  wire 
with  small  surface  resistivity.  In  1910>  Hondros  and  Debye  (2)  demon¬ 
strated  theoretically  that  it  is  possible  to  propagate  a  TM  wave  along 
a  lossless  dielectric  cylinder.  However,  due  to  the  large  field  extent 
outside  the  wire  eind  the  relatively  high  attenuation  of  this  surface 
wave,  the  "open-wire"  line  remained  a  novelty  for  almost  half  a  century. 
Recent  developments  in  the  generation  and  application  of  millimeter  and 
sub-millimeter  electromagnetic  waves,  the  availability  of  very  low  loss 
dielectrics,  and  the  development  of  fiber  optics,  have  renewed  interest 
in  the  surface  waveguides.  There  have  appeared  numerous  papers  and 
reports  concerning  various  forms  of  surface  waveguides  and  the  feasibi¬ 
lity  of  these  guides  as  practical  transmission  lines. 

Before  discussing  the  purpose  and  the  scope  of  the  present  inves¬ 
tigations,  a  survey  of  previous  work  on  surface  waveguides  is  in  order. 

1.1  Survey  of  the  Llteratvire 

The  surface  wave  guiding  structures  are  capable  of  supporting  waves 

intimately  bounded  to  the  surface  of  the  structxire.  These  waves  have 

exponential  decay  characteristics  in  regions  away  from  the  surface  and 

i&z 

are  governed  by  the  usual  propagation  function  e  along  the  axis  of 

the  structvire,  where  z  is  the  axial  coordinate  and  p  is  the  propaga¬ 
tion  constant.  For  real  values  of  p  such  waves  persist  at  arbitrarily 
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large  distances  from  the  source.  The-  steady  state  Solutions  with 

■  -iujt 

harmonic  time  dependence  e  §re  the  only  ones  considered  here. 

Of  primary  Interest,  .are  the  val^s  of  &  as  a' function  of  the  fre- 
qxtency  and  of  the  properties  of  the  guiding  system. 

The  surface  wave  guiding  system  can  take  many  forms.  The  one 
intensively  studied  in  the  past  was  a  surface  wave  structure  of 
infinite  extent  imbedded  in  an  Infltiite  uniform,  medium.  The  problem 
then  consisted  of  finding  the  solution  that  satisfied . the  homogeneous 
field  equations  and  the  boundary  condi^;lons  with  thp  source  at  infinity. 

Typical  sxirface  wave  structures  may  be  classified  into  three 
types .  The  first  type  Is  the  dielectric  coated  conductor,  such  as 
dielectric  coated  conducting  plane  emd  wire.  The  second  type  Is  the 
Interface  of  two  dielectric  media,  such  aS  dielectric  rods,  dielectric 

tubes,  or  dielectric  strips.  The  third  type  consists  of  various  open 

♦ 

periodic  struct\ares,  such  as  unbounded  helix,  corrugated  plane  or  cylin¬ 
der.  Sketches  of  these  three  types  of  surface  wave  structures  are  shown 
in  Figure  I-l.  • 

Among  the  various’  structures  mentioned  above’,  only  those  intimately 
related  to  the  propagation  of  surface  waves  along  an  elliptical  dielec¬ 
tric  cylinder  will  be  discussed  fiirther,  namely,  the  Sommerfeld-Goubau 
wire,  the  circular  dielectric  rod,  the;  dielectric  tube,  and  the  ellip- 
tical  dielectric  rod.  Related  topics  such  as  the  interaction  of  two 
surface  waveguides  emd  the  problem  of. excltatldn  of  surface  waves  will 

also  be  mentioned  briefly.  .  .  ’ 

•  *•  *  . 

•  '  *  *  . 

(a)  The  Sommerfeld-Goubau  Wire, 

The  possibility  of  propagating  a  surface  electromagnetic  wave  along 

a  circular  conducting  wire  wae. first  dempastrated  .theoretically  by 

•  .  •  •  • 


(a)  dielectric  coated 
conducting  plane 


(b)  dielectric  coated 
conducting  wire 


(c)  dielectric  slabs  (d)  dielectric  rod 


( e )  corrugated 
plane 


(f)  unbounded  helix 


Fig.  I-l.  Typical  ourface  wave  gvildlng  structures. 

Tj^je  1,  the  dielectric  coated  conductor,  (a),(b); 

Type  2,  the  Interface  of  two  dielectric  media,  (c),(d)j 
Type  3f  the  open  periodic  structure,  (e),(f). 
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Sommerfeld  (l)  in  1899*  The  wave  was  a  circularly  symmetric  TM  mode 
with  components  E  ,  E  and  was  loosely  bovuid  to  the  surface  of 

the  wire.  In  a  numerical  example  he  showed  that  the  damping  at  high 
frequency  for  this  type  of  wave  was  too  pronounced  to  use  as  a  com¬ 
munication  wave.  Consequently  the  practical  uses  of  this  type  of 
transmission  line  were  very  limited.  In  1909  Hondros  (3)^  a  student 
of  Sommerfeld,  showed  that  an  asymmetric  field  distribution  was  also  * 
possible.  But  the  wave  was  so  strongly  attenuated  that  it  could  not 
be  observed  experimentally. 

Recently  in  1950  Goubau  (^4.)  reinvestigated  the  properties  of 
the  Sommerfeld  line  and  studied  its  suitability  as  a  practical  com¬ 
munication  line.  Kis  investigation  showed  that  a  circularly  symmetric 
surface  wave  might  be  guided  by  a  conducting  wire  of  small  diameter 
with  the  same  low  attenuation  as  that  of  the  conventional  coaxial  con¬ 
ductor  guide.  However,  the  field  extended  radially  to  a  considerable 
distance  outside  the  wire  before  its  strength  decayed  to  a  negligible 
value;  so  that  any  small  imperfection  of  the  surface  or  any  small  cur¬ 
vature  along  the  wire  would  cause  radiation  loss.  The  practical  use  of 
this  siirface  waveguide  was  therefore  limited.  In  an  effort  to  reduce 
the  radial  extension  of  the  field,  Goubau  (5)  proposed  the  coating  of 
the  conductor  with  a  thin  sheath  of  dielectric, or  corrugating  the  wire. 
This  reduction  of  radial  field  extent  was  achieved  with  the  penalty  of 
higher  atteniiation  due  to  dielectric  loss  or  corrugation.  This 
Increased  attenviation  due  to  dielectric  loss  together  with  the  original 
attenuation  due  to  loss  in  the  conductor  has  been  calculated  by 
Goubau  (5).  It  should  be  noted,  however,  that  the  first  theoretical 


analysis  of  electromagnetic  wave  propagation  along  a  conducting  wire 
with  a  cylindrical  insulating  sheath  was  given  by  Harms  (5)  in  1907. 

Since  Goubau's  report,  numerous  papers  concerning  this  type  of 
single  wire  line  have  been  published.  Among  the^fe  are  the  papers  by 
Barlow  and  Karbowlak  (6)  in  1953  on  the  measurement  of  radial  field 
distribution;  Shelbe,  King  and  Van  Zieland  (7)  in  195^  on  the  measured 
losses  of  the  "Goubau  Line";  and  Roberts  (8)  on  the  excitation  of  the  . 
single  wire  line.  Kiely  (9)  also  reported  on  the  effect  of  fog  and 
rain  drops  on  the  attenuation  characteristics  of  the  wave  propagating 
along  a  long  single  wire  line. 

(b)  Circxilar  Dielectric  Rod 

Hondros  and  Debye  (2)  in  1910  analyzed  theoretically  the 'guiding 
of  a  circularly  symmetric  IM  wave  along  a  solid  lossless  dielectric 
cylinder  and  thereby  removed  the  cause  of  the  strong  attenuation  d\ae 
to  the  conductor  (l).  In  1915  Zahn  (lO)  and  his  two  students,  R'iiter 
and  Schriever  (11,12),  confirmed  the  existence  of  such  a  IM  wave  experi¬ 
mentally.  Carson,  Mead  and  Schelkunoff  (13)  noted  in  their  paper  that 
Southworth  in  1920  also  accidentally  observed  such  a  wave  in  a  troxjgh 
of  water.  When  the  generation  of  high  frequency  electromagnetic  waves 
(about  10  cm)  became  feasible,  Southvrorth  (l4)  described  some  experi¬ 
mental  work  dealing  with  phase  velocity  and  attenuation  of  the  circu¬ 
larly  symmetric  TM  wave  on  the  circular  dielectric  gxilde. 

Not  until  1936  were  the  propagation  properties  of  asymmetric 
waves  on  a  round  dielectric  rod  considered.  A  rather  complete  mathe¬ 
matical  analysis  of  this  problem  was  given  by  Carson,  Mead  and  Schel¬ 
kunoff  (13) •  It  was  noted  in  their  paper  that  in  order  to  satisfy  the 
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boundary  conditional  a  hybrid  wave  (l.e.|  the ' coexistence  of  longltu- 
.  dljftal  eJbectric  and  magnetic  fields)  must  be  assumed.  In  other  wordS| 

:  ^ynhtttrlc  ^  and  TM!  moden'  were  Inextricably  coupled  to  each  other 
: along  a  clrculsu;  dielactric  rod*  They  also  shioved  that|  l)  pure  TE 
and  TM  waves  couU  only  exist  in  the  circularly  symmetric  case,  and 
;  2) -  there  exlated  one  and  only :one  mode,  namely  the  lowest  order  hybrid 
wave  called  the  HEq  i  modci . which  possessed  no  cutoff  frequency*  and 
could  propagate  at  all  frequencies.  All  other  circularly  symmetric  or 
non^ symmetric  modes  had  cutoff  frequencies.  The  dispersion  relations 
Of  these  nodes  were  also  obtained  in  their  paperi  but  no  numerical 
results  were  given.  • 

Since  then  the  development  of  metal  tube  waveguides  as  transmis¬ 
sion  systems  conqpletely  over-shadowed  the  development  of  dielectric 
waveguides.  This  is  largely  due  to  the  fact  that  the  field  is  con¬ 
tained  entirely  within  the  metal  tube  guide.  For  the  dielectric  guldei 
hdwfevcri  the  field  is  not  entirely  contained  which  leads  to  greater 
transmission  loss  due  to  radiation  when  bends  and  discontinuities  are 
present.  A  large  number  of  papers  have  been  published  on  the  subject 
of  propagation  of  electromagnetic  waves  in  a  hollow  metal  tube.  Borgnis 
and  Papas  (15)  gave  a  very  comprehensive  treatment  on  this  subject. 

In  19^5  Mallach  (16)  published  his  res\ilts  on  the  use  of  the 
dielectric  rod  as  a  directive  radiator.  He  showed  e:qperimentally  that 
the  radiation  pattern  obtained. by  the  use  of  the  asymmetric  HE^  mode 
produced  only  one  lobe  in  the  prifaclpal  direction  of  radiation. 

*Th±8  cutoff  frequency  does  not  have  the  conventional  definition  as 
that  for  the  metal  waveguide  modes  (see  p.295  of  reference  15).  It  is 
here  defined  that  the  cutoff  freqxiency  for  the  s\irface  waveguide  mode 
is  the  frequency  below  which  the  dielectric  rod  ceases  to  act  as  a 
binding  medi\im  and  the  wave  is  no  longer  guided  by  this  svurface  wave 
structure . 


Immediately  after  Mallach's  paper,  Wegener  (17)  presented  a  dissertation 
In  which  the  asymmetric  HE.,^  mode,  together  with  the  lowest  order  cir¬ 
cularly  symmetric  TE  and  TM  modes  were  analyzed  in  detail.  Not  only 
were  the  numerical  results  of  the  propagation  constants  for  these 
waves  obtained,  but  also  their  attenuation  characteristics.  T^peirently 
he  was  not  aware  of  Carson,  Mead  and  Schelkunoff ' s  work.  A  few  experi¬ 
mental  points  were  also  included  in  his  work  to  substantiate  his 
theoretical  results.  Elsasser  (l8)  in  19^9,  independent  of  Wegener's 
work,  published  his  computation  on  the  attenuation  properties  of  these 
three  lowest  order  waves  by  the  pertvirbatlon  method  (15)*  In  a  com¬ 
panion  paper,  Chandler  (19)  verified  experimentally  Elsasser 's  results 
considering  the  dominant  HE^  mode.  He  found  that  the  guiding  effect 
was  retained  even  when  the  rod  was  only  a  fraction  of  a  wavelength  in 
diameter.  Since  the  greater  part  of  the  guided  energy  was  outside  the 
dielectric,  very  little  loss  was  observed.  For  the  first  time  the 
cavity  resonator  technique  was  introduced  to  measvire  the  attenxiation 
constant  of  the  HE^j^  mode.  The  resonator  technique  was  very  s\il table 
for  investigating  very  low  loss  uniform  waveguides.  It  should  be  noted, 
however,  that  the  formula  relating  the  Q  of  the  resonator  and  the 
attenuation  constant  a  in  Chandler's  paper  is  only  applicable  for 
very  small  2a/ ,  where  a  is  the  radius  of  the  rod  and  X^  is  the 
free  space  wavelength  (see  Chapter  V) . 

King  (20)  in  1952  proposed  the  so-called  "dielectric  image  line" 
as  a  practical  surface  wave  guiding  device.  The  "dielectric  image  line" 
was  made  of  a  semicircular  dielectric  rod  mounted  on  a  conducting 
sheet  and  was  designed  specifically  for  the  domlneuit  HE^  mode.  He 
indicated  that  the  conducting  sheet  not  only  could  act  as  a  supporting 
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devlce,  but  also  as  a  polarization  anchor  for  this  dominant  mode.  A 
detEdled  study  on  the  attenviation  and  the  reullal  field  decay  cheirac- 
terlstlcs  of  the  mode  guided  by  this  Image  line  vas  repoxted  by 
Schleslnger  and  King  (21)  In  19^3*  Again  the  cavity  resonator  method, 
used  by  Chandler,  vas  used  for  the  attenuation  constant  measurement. 

As  of  now  the  "dielectric  image  line"  Is  still  the  best  and  the  most 
practical  device  for  supporting  the  dominant  dielectric  mode. 

(c)  Circular  Dielectric  Tiobe 

A  natural  generalization  of  the  analysis  of  the  propagation  of 
electromagnetic  waves  on  a  dielectric  rod  would  be  that  for  the  circu¬ 
lar  dielectric  tube.  The  earliest  theoretical  emalysis  was  carried  out 
by  Zachoval  (22)  in  1932.  He  considered  the  propagation  of  a  circulsurly 
symmetric  TM  wave  along  a  lossless  circular  dielectric  tube.  Two  yeeirs 
later  Llska  (23)  verified  2iachoval's  work  e:q)erimentally .  A  more  com¬ 
plete  treatment  on  the  theory  of  dielectric  tube  waveguides  was  given 
by  Astraham  (2^)  in  19^9,  in  which  both  symmetric  and  asymmetric  propa¬ 
gating  waves  were  considered.  He  also  substantiated  his  theoretical 
results  by  experimental  data.  Independently,  Unger  (25)  in  1959 
reported  his  investigation  on  the  same  subject  and  shoved  that  a  dielec¬ 
tric  tube  with  a  thin  wall  coxild  sx^pport  the  dominant  mode  with  very 
little  loss.  But  the  radial  field  extent  was  rather  Isurge.  One  of  the 
most  promising  applications  of  dielectric  tube  waveguides  may  be  fomd 
in  the  field  of  millimeter  wave  cavity  resonator  and  beam  coxqpllng  struc- 
t\ire  (26) . 
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(d)  Elliptical  Dielectric  Cylinder  *.*‘  ' 

The  first  attenpt  to  find  the  dispersion  relation- of  aa  eledtro-  , 

magnetic  -wave  guided  by  an  elliptical  cylinder  struetxire /was  made  by 

Karbowiak  (27)  in  195^*  He  considered  the  elliptical'  cross-Section 

Sommerfeld  line  and  the  elliptical  cross-section  Goubau  line.  The  wave 

equation  was  formulated  in  elliptical  coordinates  and  solutions  were 

obtained.  However,  he  matched  the  boundary  conditions  only  at  one 

point  on  the  boundary  surface;  therefore  his  results  can,  at  best,  be  . . 

considered  an  approximation  for  very  small  eccentricity.  Another 

atten5)t  to  solve  the  problem  of  surface  wave  propagation  along,  an  ellip- . 

tical  dielectric  rod  was  made  by  King  and  Wiltse  (28) .  Again  they  for-  . 

mulated  the  problem  in  elliptical  coordinates  and  obtained  solutions  of 

Maxwell's  equations  in  this  coordinate  system.  But  in  matching  the 

fields  on  the  boundary,  similar  over-simplifi'cations' of  the  boiindary  . 

•  *  •  •  •.*  . 

conditions  were  made.  The  "approximation",  of  these  two '.'approaches  can 
be  best  Illustrated  by  the  following  example.'.  For  the 'sake  of .  clarity, 
only  the  matching  of  the  axial  electric  field  on  the  boundary  -will  be., 
considered.  . 

Karbowiak' s  method.  The  e;q)ression' for  the  axial  electric  .field 
in  region  1  which  is  the  region  inside  the  dielectric  rod  is.- 

•  •  • 

where  is  an  arbitrary  constant.  The  egression -for  the  social  elec¬ 

tric  field  in  the  surrounding  medium  is 

Ezo  =  cej(il,r2)  (2) 

where  is  an  axbitrsury  constant.  It  shoiHd  be  noted  that  all  these 

'•‘'The  no-tatlons  of  the  Mathieu  or  the  modified  Mathieau  functions  stre 
defined  in  Chapter  II . 
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Mathieu  functions  and  modified  Mathieu  functions  are  functions  of  the 

characteristics  of  the  medium.  The  boundary  condition  dictates  the 

continuity  of  the  axial  electric  field,  i.e.,  at  5  =  t  ,  E  =  E- 

°  ^1  o 

we  have 


(3) 


2  2 

It  should  be  noted  that  ce^(Tj,y^)  euad  functions  of 

Tj  .  The  only  way  that  equation  3  can  be  satisfied  is  by  assuming 
2  2 

®®j^(n>Ti)  =  ce*(Ti,r2)  >  which  is  not  true  except  when  the  eccentricity 
is  zero.  This  was  the  assun^jtion  made  by  Karbowiak. 


King  and  Wlltse's  method.  King  and  Wiltse  realized  the  invali¬ 
dity  of  Karbowiak' s  assun^jtion  and  proposed  to  attack  the  problem  in  a 
slightly  different  way.  Tliey  assumed  that  the  expression  for  the  axial 
electric  field  in  the  dielectric  rod  is 


where  the  are  the  arbitrary  constants;  and  the  expression  for  the 

axial  electric  field  in  the  surrounding  medium  is 


E,  -  \  EeE  (5,r|)  ce*(„T|)  (5) 

O 

Where  L  is  an  arbitrary  constant.  Satisfying  the  boundary  condition 
n 

at  I  =  f  »  we  have 
^  ^0 


p 

They  then  multiply  both  sides  of  the  equation  by  ce^(T),T'^)  and 


integrate  with  respect  to  from  0  to  2k  ,  obtaining 


(7) 


2n  2n 

/2  2  12  2 

ce  (il,r,  )dTj  and  M  =  ce  (ii,r, )  ce*(Ti,r5)dT|  . 
ni  njnj-nc 


This  was  how  they  eliminated  the  summation  sign.  It  can.  be  seen  that 
an  identical  result,  i.e.,  equation  7,  can  be  obtained  by  the  use  of 

2 

equations  1  and  2  .  Mxiltiplying  both  qides  of  equation  3  by  ce^(n,y'^) 
and  integrating  from  0  to  2k  ,  one  obtains  equation  7. 

Therefore  the  validity  of  King  and  Wiltse's  solution  is  also  ques 
tionable . 


(e)  Related  Topics 

Unlike  the  waves  in  the  metal  tvibe  waveguides,  there  are  no 
evanescent  modes  on  an  open  surface  waveguide.  It  is  not  possible  to 
express  any  arbitrary  field  distribution  in  terms  of  the  propagating 
modes  alone.  Hence,  there  must  exist  a  different  type  of  wave,  namely 
the  radiated  wave  (29)  if  any  source  is  present  in  a  finite  region. 

As  a  matter  of  fact,  it  should  be  noted  here  that  the  presence 
of  the  surface  wave  was  actually  first  postulated  by  Sommerfeld  (30) 
in  1909  when  he  was  considering  the  now  classical  problem*  which  bears 
his  name.  He  foxind  theoretically  that  there  existed  not  only  a  radiated 
wave  due  to  the  oscillating  dipole,  but  also  a  surface  wave  which 
traveled  along  the  interface  of  the  two  dielectrics.  Since  then,  a 
great  number  of  papers  and  reports  have  been  published  concerning  varia¬ 
tions  of  this  problem.  The  most  recent  investigations  have  been  reported 

*The  Sommerfeld  problem  is  discussed  very  clearly  and  thoroughly  in 
Stratton  ( 3I) • 
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by  Roe  (29),  Whitmer  (32),  Tai  (33),  Brick  (3*+),  Walt  (35),  Cxillen  (36), 
and  Brown  (37),  to  mention  only  a  few. 

The  problem  of  interaction  between  two  parallel  uniform  surface  • 

waveguides  is  also  an  interesting  one.  Since  the  wave  equation  is  not 
separable  in  the  bipolar  coordinates,  approximate  methods  must  be 
employed.  Quite  a  few  authors  used  the  electrostatic  approximation*  in 
the  earlier  years.  However  this  approximation  was  not  satisfactory  at 
very  high  frequencies.  Most  recently  Armand  (38)  and  Marcuse  (39) 
treated  the  problem  of  interaction  between  two  parallel  Goubau  wires 
without  resorting  to  the  electrostatic  approximation.  They  formulated 
the  problem  by  assuming  the  interaction  of  only  one  single  mode  on  each 
wire,  namely  the  circularly  symmetric  IM  mode.  They  indicated  the 
presence  of  space  beats  and  the  energy  exchange  phenomenon.  Numerical 
exan^iles  were  also  given. 

1.2  Purpose  and  Scope  of  the  Present  Investigation 

In  order  that  the  dielectric  rod  may  be  a  low  loss  svurface  wave 
device,  one  m\ist  choose  a  small  value  of  ka  where  k  is  the  free 
space  wave  number  and  a  is  the  re^ius  of  the  dielectric  cylinder.  In 
the  millimeter  wavelength  range,  the  radius  of  the  dielectric  cylinder 
becomes  inconveniently  small.  Fortunately  it  has  been  found  experimen¬ 
tally  (4l)  that  if  the  circular  rod  is  flattened,  (i.e.,  if  the  circular 
rod  is  rendered  to  an  elliptical  rod  of  the  same  area) ,  the  attenuation 
of  the  dominant  mode  may  be  reduced  considerably,  provided  that  the 

■*The  electrostatic  approximation  is  as  follows:  In  calculating  the 
structure  of  the  field  one  would  neglect  quantities  of  the  order  of 

0(d  \/k^)  ,  where  d  is  the  distance  between  the  wires,  kQ  is  the 
wave  nmber  for  free  space,  eind  k  is  the  propagation  constant  of 
the  wave. 


electric  field  of  the  dominant  node  Is  parallel  at  the  center  of  the 
rod  to  the  minor  axis  of  the  elliptical  rod.  The  use  of  very  thin 
fibers  of  various  cross-section  as  optical  wavegtiides  or  as  mode  selec¬ 
tors  in  optical  masers  has  also  received  considerable  attention.  [For 
example >  see  reference  (40)].  Ftirthermore,  it  is  noted  that  so  far 
there  exists  no  satisfactory  way  of  analyzing  the  problem  of  surface 
wave  propagation  along  a  dielectric, rod  of  elliptical  cross-section. 

It  is  therefore  the  purpose,  of  the.  present  investigation  to  develop  a 
method  to  analyze. .this  problem  theoretically,  to  examine  in  particular 
tte  propagation  characteristics  of  the  dominant  nodes,  and  to  perform 
experiments  to  .verify  the  analytic  results. 

•  The  investigation  is  divided  into  six  parts,  and  the  results  are 
correspondingly  presented  in  Chapters  II,  III,  IV,  V,  VI  and  VII.  In 
C,hapter.II  the  fundament.ai  theory  of  wave  propagation  along  an  ellipti- 
cal  dielectric  rod  is  given.  A  method  is  developed  to  assure  that  the 
solutions  of  the  wave  equation  satisfy  all  the  boundary  conditions  on 
ttie  s\u:face  of  the  dielectric  rod.  The  characteristic  equations  for 
the  principal  modes  are  given  sO  that  the  variation  of  guide  wavelength 
with  frequency,  the  dielectric  constant,  and  the  physical  dimensions  of 
the  guide  may  be  obtained.  It  is  shown  analytically  that  there  exist 
two  non-degenerate  modes  which  possess  no  cutoff  frequency.  They  are 
called  the  dominant  .mbdes,  and  it-  is  the  propagation  characteristics  of 
these  that  will  be  considered  in  detail  in  the  subsequent  chapters.  It 
is  also  shown  that  all  the  principal  modes  on  an  elliptical  dielectric 
rod  degenerate  smoothly  to  the  well  known  modes  on  the  circular  dielec¬ 
tric  guide,  as  the  eccentricity  of  the  elliptical  rod  approaches 


zero. 


Numerical  results  of  the  characteristic  equations  for  the 
dominant  modes  are  obtained  and  discussed  in  Chapter  III.  Sketches 
of  the  field  configurations  are  also  given.  The  decay  characteristic 
of  the  axisil  electric  field  is  computed. 

In  Chapter  IV  the  attenuation  properties  and  the  power  distribu¬ 
tion  characteristics  of  the  dominant  modes  are  analyzed  theoretically 
with  the  assuji^Jtion  that  the  dielectric  loss  is  small.  Numerical 
results  are  computed.  It  is  foiond  that  the  attenuation  constant  of 
the  dominant  mode*  propagating  along  an  elliptical  dielectric 

rod  is  much  less  than  that  of  the  dominant  mode  along  a  circular 

dielectric  rod  having  the  same  cross-sectional  area.  Physical  inter¬ 
pretation  of  these  results  is  also  presented. 

The  Q's  of  an  elllptlceil  dielectric  rod  cavity  supporting  the 
dominant  modes  are  given  in  Chapter  V.  It  is  shown  that  very  high  Q 
cavity  may  be  constructed  using  thin  elliptical  dielectric  rod.  Also 
derived  is  a  formula  relating  the  Q  of  a  cavity  and  the  attenuation 
constant  of  a  transmission  line  supporting  the  same  mode.  This  for¬ 
mula  is  more  general  than  the  one  given  by  Davidson  and  Slmmonds  (4l) 
in  that  it  is  also  valid  for  the  hybrid  modes.  This  relation  is  very 
important  whenever  the  cavity  resonator  method  (19)  is  used  to  measure 
the  attenmtlon  constant. 

To  verify  the  theoretical  resvilts  a  systematic  experimental 
investigation  on  the  propagation  characteristics  of  the  two  dominant 
modes  was  performed.  A  detailed  description  of  the  measuring  apparatus 
and  technique  is  presented  in  Chapter  VI.  Experimental  data  are  then 


*The  meaning  of  this  symbol  is  given  in  Chapter  II. 
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compared  with  theoretical  results,  and  they  are  in  very  good  agreement. 

Summary  and  conclusions  are  given  in  Chapter  VII.  The  advanteiges 
of  using  a  flat  elliptical  dielectric  rod  instead  of  a  circ\ilar  dielec¬ 
tric  rod  ua  a  microwave  guide  are  pointed  out.  It  is  also  indicated 
that  the  analytic  method  used  here  may  he  applied  to  other  similar  prob¬ 
lems.  The  problem  with  source  present  is  also  discussed  briefly. 
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CHAPTER  II  -  THEORY  OF  ELLIPTICAL  DIELECTRIC  WAVEGUIDES 

The  problem  Is  formulated  In  terms  of  the  elliptical  cylinder 
coordinates;  the  appropriate  solutions  of  the  wave  equation  in  this 
coordinate  system  are  then  obtained.  The  difficulties  of  satisfying 
the  boundary  conditions  on  the  elliptical  surface  are  pointed  out.  A 
method  to  overcome  such  difficulties  is  introduced.  Various  notations 
and  classifications  of  the  principal  propagating  modes  are  defined. 

Upon  matching  the  boundary  conditions  by  the  indicated  method,  a  set  of 
characteristic  equations  and  explicit  forms  for  all  field  components 
corresponding  to  various  modes  €ire  obtsdned.  The  existence  of  the 
dominant  modes  having  no  cutoff  frequency  is  demonstrated.  Finally,  it 
will  be  shown  that  as  the  eccentricity  approaches  zero,  all  principal 
propagating  modes  degenerate  to  the  well  known  circular  modes. 

2.1  Formulation  of  the  Problem 

The  surface  wave  propagation  along  an  infinitely  long,  straight, 
isotropic,  and  homogeneous  dielectric  cylinder  of  elliptical  cross  sec¬ 
tion  Imbedded  in  sua  infinite  dielectric  medium  of  dielectric  constant 

€  and  msignetlc  permeability  n  ,  is  considered.  The  dielectric 
o  0 

cylinder  has  a  dielectric  constant  emd  a  magnetic  permeability 

.  We  assume  that  ,  the  free  space  magnetic  permeability;  ■ 

,  suid  that  the  conductivity  in  both  media  is  zero.  We  fxirther.  . 
assume  that  the  exciting  soxurce  is  so  far  away  that,  in  the  region  of  ■' 
interest,  the  surface  waves  dominate  the  rewiiated  waves  from  the  source’. 

To  analyze  the  source-free  dielectric  surface  waveguide  of  elllp-’  .• 
tlcal  cross  section,  the  elliptical  cylinder  coordinates  (|,tj,8),  as 
shown  in  Figure  II-l,  are  introduced.  The  elliptical  cylinder  coordinates 
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Flg.II-l  (a)  Cross-section  of  elliptical  surface  wavegvilde. 

and  sire  the  foci  of  the  ellipse.  The  dlstsince 
betveen  foci  is  the  focal  distance,  2q  . 

(h)  Degenerate  form  of  ellipse  when  e  =  1.  As  e  -►  1 
semi-minor  axis  Q  ,  and  semi -major  suds  -»  q  . 

(c)  Degenerate  form  of  ellipse  when  e  =  0.  As  e  -►  0 

q  -*  0,  -*  oo  ,  semi -major  suds  -►  semi-minor  axis  -►  r  . 
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are  related  to  the  rectangxilar  coordinates  (x',y',z')  through  the 
following, 

x'  =  q  cosh  I  COB  tj 
y’  =  q  sinh  |  sin  t] 
z  •  =  z 

(0^|<CD,  0  ^  ^  2j() 

where  q  Is  the  semlfocal  length  of  the  ellipse.  The  contour  surfaces 
of  constant  |  are  confocal  elliptic  cylinders,  and  those  of  constant 
q  are  confoceLL  hyperbolic  cylinders.  The  elliptic  cylinders  and  hyper¬ 
bolic  cylinders  have  foci  at  x'=  q,  y'=  0  and  x'=  -q  ,  y'*  0  .  The 
semi-minor  axis  is  equal  to  q  sinh  .  The  eccentricity  e,  defined 

as  the  ratio  of  the  semlfocal  distance  to  the  semi-major  axis,  is  given 
by  e  =  l/cosh  . 

One  of  the  confocal  elliptic  cylinders  with  I  =  is  assumed  to 
coincide  with  the  boundary  of  the  solid  dielectric  cylinder,  and  the 
z-axls  coincides  with  its  longitudinal  axis. 

2.2  Maxwell's  Equations  and  Their  Solutions  in  Elliptical  Cylindrical 
Coordinates  ' 

It  is  well  known  that  the  harmonic  form  of  Maxwell's  equations  in 
a  source-free  medliim  characterized  by  e  and  n  are  given  by 


V  X  E  =  i  H  ( la) 

V  X  H  *  =  -  ia)€  E  (lb) 

V  •  H  =  0  (Ic) 

V  .  E  =  0  (Id) 


where  E  and  H  sure  the  electric  field  vector  and  the  metgnetic  field 
vector  respectively.  The  harmonic  time  dependence  of  e"^^ 


for  all 
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fleld  quantities  Is  assiomed.  The  rationalized  MKS  system  Is  used 

throughout  this  work.  We  shall  now  confine  our  treatment  to  waves 

propagating  ed.ong  the  positive  z-axLs.  In  complex  representation 

these  assun^jtlons  restilt  In  a  multiplication  of  all  wave  functions  by 
-ioit  ipz 

6  0  }  li*0*^ 


E(|,Ti,z,t)  =je^E^(|,Ti)  +e^E^(|,Ti)  e  e^^*  (2) 

=je^H^(5,7j)  +e^H^(6,T))  +e^H^(|,Tl))  e  (3) 


where  e  ,  e  and  e  are  tmlt  vectors  In  the  z  directions  res- 

-6'  -Ti  -z 

pectlvely,  emd  P  ,  the  propagation  constant  of  the  wave  In  the  z  direc¬ 
tion,  Is  to  be  determined  from  the  boundary  conditions. 

In  elliptical  cylinder  coordinates,  equations  la  and  lb  become 


-  la)€  p  E 

z 

'  ^  V  ■  “i- 

(4) 

-  la)€  p  Eg 

=  li  <».)  -  ipp 

(5) 

-  lai€  p  E^ 

(6) 

1 0)  H 

°  It  V  ■  ^ 

(7) 

1(1)^  pHg 

=  |j  (E^)  -  IE  p 

(8) 

1  oi  |i  p 

=  '  lepEj  (E^) 

(9) 

,  2  ,  2  .1/2  2  2  /2jt.2  .  .  4 

where  p  =  q(slnh  i  +  sin  tj)  ,  k  =  co  pic  =  (— )  and  X  Is  the 

wavelength  of  a  uniform  plane  wave  In  the  medium.  The  above  eqtiatlons 

4  through  9  can  be  combined  to  give  the  field  components  E^,  E^,  H^, 

In  terms  of  E  and  H  only;  we  have, 
z  z 


E. 


■h-r  W 

(k  -  e^p  I. 


dE 

2 

“5? 

dE 


z  .  •  z 


1  C  •  '^“z 

=  — g - ^  "jip  -  ia>^  V 

/k  -  6- )p  L  J 


H 


i 


H. 


<k  -  r)p 

1  ■  ' 

(k^-  p^)p 

(k^-  e^)p 


dE^ 

dE. 


-  i  uj€  —5^  + 


Sh 

! 

"5; 

■11 

^H, 

"§i 


.ia,€-^,-  ip-^ 


(10) 

(U) 

(12) 

(13) 


Taking  the  derivative  of  equation  12  with  respect  to  t)  and  the  derl- 

•  • 

vatlve  of  equation  13  with  respect  to  |  and  substituting  these 
expressions  into  equation  U^one  obtains  the  equation 


■  2  2 

^  E  .  S  E  2  ■  2  •  2  2  2 

- o  +  — o"  ■*■  [‘I  (k  -  0.  )(sinh  I  +  sin  t))  ]  E  ■  =■  0  . 

Si,®  .  .  k  . . . 


(Ik) 


Similarly,  teiklrig  the  derivative  of  equation  10  with  respect  to  t)  and 
the  derivative  of  equation'll  with, respect  to  |  anti  substituting 
these  expressions  into  equation -7,  one  gets 


Equations  l4  and  15  are  the  wave  equations.  It  should  be  noted  that 
these ‘two  wave  equations  are  of  the  same  form,  therefore-  it  is  only 
necessary,  to  solve  one  of  them.  If  =  0  a  TM  wave  results;  if 
E  =  0  a  TE"  wave  results.  The  ^st  general  e}q)ressi6ns  for 'the  elec-- 

Z  * 

troniEignetlc  fields  consist  of  a  linear  combination- of  the-  solutions  of 
TE  and  TM  waves* 

.  Consider-  the  ..following  partial  differential  equation  ' 
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P^)(slnh^|  +  sin^Tj)]A  =  0  (16) 

drj  *- 

in  TiAilch  A  may  be  H  or  E  .  In  order  to  obtain  the  solutions  of 

z  z 

equation  16  one  sets 

A(|,tj)  =  R(|)  0(ti)  (17) 

and  substitutes  equation  17  into  equation  16.  Applying  the  usual 
separation  of  variables  procedure,  one  may  separate  equation  16  into 
the  following  two  ordinary  differential  equations 

-  -  +  (c  -  2r^cos  27j)@(h)  =  0  (l8) 

atj 

and 

-  (c  -  2t^cosK25)  R(|)  =  0  (19) 

dr 

2  2  2  2 

where  c  is  the  separation  constant  and  y  =  (k  -  3  )q  /4  .  Equation 
l8  is  the  Mathieu  differential  equation;  equation  19,  which  follows 
from  l8  by  the  transformation  tj  =  ±  i|  ,  is  the  modified  Mathieu 
differential  equation  (^2). 

For  physically  admissible  single-valued  electromagnetic  fields, 

A(ril)  must  be  a  periodic  function  of  tj  ,  of  period  jt  or  2jt  , 

2 

and  the  separation  constant  c  ,  in  this  case  a  function  of  y  ,  takes 

2  2 

on  an  infinite  set  ‘of  characteristic  values  for  every  y  .  When  y 

is  real  the  characteristic  values  are  all  real;  since  we  are  considering 

2 

solutions  in  a  lossless  medium,  only  real  values  of  c  and  y  are  of 

2 

interest.  Corresponding  to  y  =  0  there  are  two  independent  periodic 

solutions,  namely  sin  nij  and  cos  nij  with  the  separation  constant 

2  2 
c  =  n  where  n  is  an  integer.  It  can  be  shown  (43)  that  when  y 
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differs  from  zero,  a  characteristic  value  c  determines  one  and  only 
one  periodic  solution  which  is  either  even  or  odd  in  tj  .  The  charac¬ 
teristic  val\aes  c  ,  giving  rise  to  even  and  odd  solutions  eure  denoted 

•  2  2 

here  by  )  and  )  respectively.  The  subscript  n  identl- 

2  2 

fies  those  sets  of  characteristic  values  which  approach  n  as  y 

approaches  zero.  It  is  known  from  the  Sturmian  theory  of  second  order 

linear  differential  equations  that  solutions  associated  with  ) 

2 

and  b  (t  )  have  n  zeros  in  the  interval  0  ^  tj  ^  jf  (44)  . 

2 

For  arbitrary  positive  real  values  of  y  ,  the  periodic  solutions 
of  Matheiu's  equation  l8  are* (42, 45) 


(9(n)  = 


sen(Tl>r^) 


(even)  a^(r  ) 
(odd) 


(20) 


and  the  corresponding  sol’.tions  for  the  modified  Mathieu's  equation  19 
are* 


R(l) 


®1  ®2^ey^(|>T^)' 

^1 


(even) 

(odd) 


a  (/) 
n' 


(21) 


2 

For  arbitrary  negative  real  values  of  y  the  periodic  solutions  of 
Mathieu's  eqxiation  l8  are* 


d(n)  = 


(even) 

(odd) 


2 

(Snd  y  I)  when  n  even) 
(bn(ly^l)  when  n  odd) 

(bn(ly^l)  when  n  even) 

2 

(a  (ly  I)  when  n  odd) 
n 


and  the  corresponding  solutions  for  the  modified  Mathieu's  eqxiation  19 


*See  Appendix  A  for  the  definitions  and  series  expansions  of  these 
Mathieu  and  modified  Mathieu  functions. 
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are* 

R(l) 


rc^Ce*(|,lr^|)  +  C2Pek^(S,|r^|) 
[d^Se*( 1, 1 r^l )  +  dgGek^C 6, |  I ) 


(even) 


(odd) 


(“ndr^D 

(\(lr^l) 

(a^dr^l) 


When  n  even) 
when  n  odd) 
when  n  even) 
when  n  odd) 
(23) 


2  2 

)  and  ^jj(t  )  are  the  characteristic;  values  and  n  is  the  order 
of  the  fianction.  a^,  c^,  Cg,  and  dg  are  the  arbitrary- 

constants. 

The  proper  choice  of  the  above  solutions  to  represertt  the  electro¬ 
magnetic  field  of  an  elliptical  dielectric  cylinder  depends  upon  the. 
boimdary  conditions.  For  region  1,  which  is  the  space  inside  the  dielec¬ 
tric  rod,  all  field  components  must  be  finite.  For  region  0,  which  is  • 
the  space  outside  the  dielectric  cylinder,  in  order  that  energy  flow  only 
along  the  axis  of  the  cylinder,  all  field  coii?)onents,  must  approach  zero 
as  the  radial  argument  approaches  infinity.  Consequently  we  must  discard 
the  functions  Fey^(|,T  )  and  Gey^(6,T  )  ,  since  they  are  Infinite  at  • 
the  origin,  l.e.,  at  |  =  0  .  The  functions  Ce*(|,ly  j)  and  Se*('|,(.T  |) 
must  eilso  be  discarded  since  they  become  infinite  at  infinity.  Therefore 
the  solutions  of  the  wave  equations l4  and  15  are  as  follows; 


H  (|,TJ,Z,t) 
^1 


A  Ce  (|,rf) 
n  n  '  1 


ce 


-ioit 


iPlz 
e  • 


00 


A^Se^(|,r5  se^(Ti,ri)e“^-^.  (I  V|i0)  (24) 

n=l  .  “ 

(M,z,t)  =  ^  L/ek^(S,lrol)ce*(,,,lrol) 


00 


.  Y.  (coMM„)  (25) 

_ n=l _ : _ _ 1_ 

-•^We  follow  the  notation  adopted  by  McLachlan  (47)  and  Meixner  (44). 
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z.  t-j-nn  1  n  i 

1  n  *  0 

E  (I  i  I  ^  0)  (26) 

n  =  1 

E^  (l,H,^,t)  =  ^  P'Fek  (I,  ITqI)  ce*(Ti,|rQl)  e 

o  .  n=  0 

■*■  E  se*(n,  IT^I  )e  e  °  ,  (c»^5^|^). 

n  =  1 

(27) 

A,A',BiB',L,L',P  and  P'  are  coefficients  which  are  related 
n  n  n  n  n  n'  n  n 

2  1  21 

by  the  boxmdary  conditions  £uad  are  f motions  of  n,  cs,  ITj^l  >  and 

the  nature  of  the  exciting  sources,  but  independent  of  the  coordinates. 

and  are  respectively  (k^-  P^)q^/4  and  l(kQ-  /h\ 

2  2  2  2 

with  k,  =  (u  li€,  and  k  =  o)  u€  .  e.  is  the  dielectric  constant  of 
i  1  0  o  1 

the  cylinder  and  is  the  dielectric  constsmt  of  the  surrounding 
medium.  5  -  surface  of  the  dielectric  cylinder.  All  trans¬ 

verse  field  con^jonents  for  both  regions  can  be  derived  from  equations 
10  through  13,  using  eqmtions  2h  through  27.  Incidentally,  the  Hertz 

vectors  jt'  and  «"  (15)  rather  than  E  and  H  may  be  used  as  the 
z  z  z  z 

scalar  quantities  from  which  the  other  field  coii?>onents  may  be  derived. 

2.3  The  Bomdary  Conditions 

The  task  of  solving  an  electromagnetic  wave  boundary  value  prob¬ 
lem  is  to  find  finite  €Uid  single-valued  solutions  which  satisfy  the 
source-free  Meocwell's  eqmtions  and  the  bomdary  conditions.  The  bom- 
dary  conditions  are,  that  the  tangential  con^onents  of  the  electric  and 
magnetic  fields  must,  in  general,  be  continuous  through  any  surface.  If 
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the  region  of  interest  Is  infinite,  then  the  radiation  condition  (46) 
must  also  be  satisfied.  The  above  conditions  are  necessary  emd  suffi¬ 
cient.  In  the  present  problem,  the  continuity  conditions  in  the 
elliptical  cylindrical  coordinates  are 


and 


E  =  E 
z.  z 
1  0 


H  =  H 
z.  z 
1  0 


E  =  E 

H  =  H 

^1  To 


(1) 

(2) 

(3) 

(4) 


5  ”^o  ’  TJ-O  and  +00  >  z  >  -00  . 

In  order  to  illustrate  the  difficulties  encountered  in  satisfying 
the  above  boundary  conditions  for  the  elliptical  dielectric  cylinder, 
we  shall  first  consider  the  case  of  the  surface  wave  propagation  along 
a  circular  dielectric  cylinder.  The  required  axial  electromagnetic 
fields  both  Inside  and  outside  tlie  circular  dielectric  cylinder  are  (I3) 


00 


Ihz  -iu)t 


E  =  )  A  J  (C  r)  cos  n©  e  e 
z  ^  n  n  ’ 
n  =  0 


00 


ikz  -ia>t 


_o  .0  ,,  ,,o  V  „  ixz  - 

L  ®  ® 

n  =  0 

^  ..  /.i  \  .  „  i4z  -icut 

ir  =  )  B  J(^r)  sin  n©  e  e 
z  n  n ’  ' 

n  =  1 


00 


ikz  -i(ut 


H°  =  y  B°  K  (^°r)  sin  n©  e  e 
z  /Li,  n  n  ^ 
n  =  1 


where 


=  /k^-  k^  and  5° 


=  l/k^-  k^ 
0 


Ml 

U 

Ml 

0 

a) 

(5) 

(a  ^  r  < 

00) 

(6) 

(0  ^  r  ^ 

a) 

(7) 

(a  ^  r  < 

00) 

(8) 

2 

and 

,  2  2 

= 

ko  =  u)  , 

is  the  dielectric  constant  of  the  cylinder,  is  the  dielectric 
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constant  of  the  surrounding  medium,  and  £.  >  e  .  A^,  A°,  and 

^  *  4  o  n'  n'  n^ 

B°  eure  the  arbitrary  constants  and  a  is  the  radius  of  the  cylinder. 
The  boundary  conditions  are 

-  E°  (9) 

hJ  -  (10) 

=  E“  (U) 

(12) 

at  r  =  a,0^0^2«  and  -oo  ^  z  ^  oo  .  Substituting  equations  5  and 
6  into  equation  9,  one  obtains 

^  A^  J^(5^a)cos  nO  e^^^  ®  I]  ^  K^(5*^a)coB  n©  e^^^.  (I3) 

n  =  0  n  *  0 

Kxiltiplying  both  sides  of  equation  13  by  cos  a6  and  integrating  with 
respect  to  0  from  0  to  2it  we  have,  due  to  the  orthogonality  of  the 
trigonometric  functions, 

.A^Jn(cM  =  A°KjC°a)  .  (l4) 

It  should  be  noted  that  for  each  mode  (in  this  case  for  each  n  )  there 
shoifLd  be  only  one  propagation  constant.  Equation  l4  shows  that  the 
boundary  conditions  may  be  satisfied  for  each  n  separately,  due  to  the 
orthogonality  in  0  of  the  fundamental  solutions  and  the  fact  that  the 
angulsu:  function  (cos  nO  or  sin  n©)  is  independent  of  the  character¬ 
istics  of  the  medl\jm.  Similar  procedures  and  conclusions  can  be  applied 
to  the  bovuadary  conditions,  equations  10,  11  and  12. 

Consider  the  boundary  condition,  equation  1,  for  the  elliptical 
dielectric  cylinder.  Substituting  equations  2.2-26  and  2.2-27  into  1 
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'one  gets 


00 

E 

n  =  0 

00 

E 

r  =  0 


n  =  1 


Ipiz 


00 


o 

00 

r  =  1 


IPoZ 


(15) 


Equation  15  may  be  vritten  as  two  separate  equations,  one  corresponding 
to  the  even  type  modes.,  the  other  to  the  odd  type  modes.  These  equations 
are 

"  -  .2, 


n=  0 


E  P;Fek^(lo^l  ^o' ^  ce*(Ti,lr^l.)  e 

r  =  0 


(16) 


and 


00  •,  ■  p-  „  ip  z 

E  se  (TJ,rr)  6 

nnol  n  X 

n  =  1 


?  P  Gek  (s  lr^l)8ej(i,l,^|)  e 
r  =  1. 


(17) 


Suppose  one  multiplies  both  sides  of  17  by  se  (ti,t,  )  and  integrates 

ml 

with  respect  to  .  Tj  from  0  to  2n  .  Due  to  the  orthogonality  of  the 
Mathieu  functions  (see  Appendix  A),  equation  17  becomes 


2,  2 


5 

.2 1 N  _  o 


N  B  Se  U  .>rf)e  =  T  P  Gek  (|  ,  r  )  e 
mm  .m  r  r"’o.  I'o'^ 

r  =  1  , 


r 


3®*  (n^lr^l)  sejq,rj)dq 
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where  N  Is  the  normalization  constant, 
m 

B,  =  3  one  gets 
*^1  0 

CD 

N  B  Se  (|  ,Tf)  =  r  P  Oek  (»  ,ir^l) 
mm  m'^o^'l^  r  r'^o^l  o'^ 

r  *  X 

Equation  l8  Involves  the  arhitrary  constant  B  (m  =  1,  or  2,  or  3>  or 

m 

-t  j  • )  and  an  infinite  nvmiber  of  arbltrsonr  constants  P, ,  P.,  P,  •  •  •  P  . 

1'  2'  3  GO 

Similar  procedures  may  he  applied  to  the  remaining  boundary  conditions, 
equations  2,  3  and  4,  and  each  of  them  contributes  an  arbitrary  constant 
on  the  left  hand  side  of  the  eqviation  and  eui  Infinite  number  of  arbi¬ 
trary  constants  on  the  right  hand  side  of  the  equation.  For  exan^>le, 

using  2,  an  algebraic  equation  involving  A  (m  =  1,  or  2,  or  3>  or  *  •  •) 

m 

and  L-,  L.,  •  •  •  L  results;  using  3,  an  algebraic  equation  involv- 

123  00 

ing  A  and  B  (m  *  1,  or  2,  or  3,  or  •  •  •  )/  and  P, ,  P-,  P.,  •  •  ♦  p 
mm  1  2  3  00 

smd  Lg,  results;  using  k,  another  algebraic  eq\xation 

involving  A  and  B  (m  =  1,  or  2,  or  3,  or  •  •  •  )^  and  P  ,  P  ,  P  . 

mm  X  2  3 

•  •  •  P  and  L, ,  L.  •  •  •  L  results ,  Since  these  equations 
00  1  2'  3  00  ^ 

involve  eui  infinite  number  of  arbitrary  constants,  an  infinite  set  of 
linear  algebraic  equations  is  required.  This  means  m  must  be  equal 
to  0,  then  1,  then  2,  •  •  then  oo  .  It  can  therefore  be  seen  by  the 
method  outlined  above  that  in  matching  the  boundary  conditions,  an  in¬ 
finite  order  of  Mathleu  functions  must  be  used  to  describe  the  fields 
in  both  madia,  i.e.,  both  inside  and  outside  the  elliptical  dielectric 
rod. 

To  overcome  the  above  difficulties  one  assumes  that  only  one  term 
of  the  Mathleu  function  could  represent  the  field  configurations  in  one 
medium,  while  an  infinite  series  of  Mathleu  functions  in  the  other. 
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Equatloti  17  can  be  written  as 


IP-Z  00 


'■*  r=  1 


2 

^l)e  ° 


y  B  Se  (|  ,Tr?)se  (q,r?)e  ^  =  P  Gek  (5  Jr^  )se*(ii,  )e  °  . 
n  n  o  1  n  '  1  r  r  o  '  o'  r  o  ' 

n  =  1 


In  equation  19  one  term  of  the  Mathieu  function  has  been  vised  to  repre¬ 
sent  the  field  configuration  inside  the  dielectric  rod,  and  in  equation 
20  one  term  of  the  Mathieu  function  has  been  used  to  represent  the  field 
configviration  outside  the  dielectric  rod.  Consider  equation  I9.  Setting 

p 

and  multiplying  both  sides  of  I9  by  se^(Ti,Y‘p  and  integrating 
with  respect  to  t)  from  0  to  2jt  ,  one  obtains 

N  B  Se  (1  ,r?)  =  y  P  Gek  (1  ,1  r^l)  (m  =  0,1,2,  ••  •)  (21) 

nra  n  n  o  1  r  r  0  '  o'  rm  '  >  f 

r  =  1 

2k 

r  2  2 

where  N  =  |  se  (|  ^T,  )dii  ,  when  n  =  m 
nm  J  n'^o  1  '  ' 


when  n  / 


dn 

nn  I  r  o  o  mol 


For  each  value  of  n  there  exists  an  infinite  set  of  linear  algebraic 
equations  \rtiich  may  be  combined  and  simplified  to  give  an  equation  in¬ 
volving  only  two  eurbitrary  constants,  B^  and  ,  i.e., 

where  F„(p _ )  is  a  function  of  B 

n  rm  '^rm 


As  an  exan^le,  suppose  n  =  equation  21  becomes 


,  (m.0,1,,2,-0  (23) 


r  =  1 


2« 


/2  2 

f  m  =  1  , 
m  ^  1  • 


0 
=  0 


The  infinite  set  of  linear  algebraic  equations  from  equation  23  are 


t  P50ek^(5„,|rJ)^5^  . 


0  =  Pj^oek^Cgr^Ds^j-*  P3Gek^(j^,|r^|)Pjj  ♦ 


0  =  Pl,aeg5„,lr;l)^l5+  P3(iek3(S^,|r^  1)^35  + 

*  P30ek3(|^,|r^|)P33  + 


(pita) 


(Pko) 


Combining  the  above  set  of  equations  one  obtains 


(25) 


where 


= 


^33 

.  ^53 

•  •  • 

^35 

•  •  • 

»55 

•  •  • 

•  •  • 

•  •  • 

.  Pn 

^31 

^51  • 

^13 

^33 

^3 

'  ®15 

•  •  • 

^35 

•  •  • 

®55  •. 

•  •  •  • 

(26) 
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P^,  P^;  •  •  •  P^  can  all  be  expressed  in  terms  of  P^^  ;  for  exanple. 


and 


GeVSo'I’-ol) 
- rrr  ' 


Oeki(|,|r„l) 


(-Pi) 


(27) 


(28) 


The  above  infinite  determinants  may  be  solved  by  the  method  of  suces- 
sive  approximations  (47).  It  maybe  shovn  that  P^  >  P^  >  >  *  *  *  ^  ^oo 

in  this  exanple. 

The  method  described  above  will  be  \ised  in  a  later  section  to 
satisfy  the  boundary  conditions  and  to  obtain  the  characteristic  equa¬ 
tions  of  the  principal  propagating  modes.  As  one  may  anticipate,  the 
mode  classifications  are  much  more  Involved  in  the  case  of  the  ellip¬ 
tical  dielectric  cylinder  than  the  circuleu:  dielectric  cylinder.  This 
will  be  discussed  In  the  next  section. 


2.4  The  Notations  and  Classifications  of  the  Propagating  Modes 

For  a  circiilar  dielectric  waveguide  it  is  well  known  that  the 
pure  TE  and  TM  waves  can  exist  only  if  the  fields  are  independent  of  the 
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angiilar  coordinates.  These  circularly  symmetric  waves  are  designated 

hy  H  for  the  pure  TE  waves  emd  E  for  the  pure  TM  waves.  The 
'on  on 

subscript  o  signifies  the  ang\ilar  variations  and  n  signifies  the 
nth  root  of  the  chcuracterlstlc  equation.  The  coexistence  of  E  and 
H  waves  Is  required  to  satisfy  the  boundary  conditions  If  the  field 
Is  a  function  of  the  angular  coordinate  .*  These  asymmetric  waves  are 
then  designated  by  HE  if  the  cross-sectional  field  pattern  resem- 

Tnn 

bles  that  of  an  H  wave  and  by  EH  If  the  cross-sectional  field 

mn 

pattern  resembles  that  of  an  E  wave.  The  subscripts  m  and  n  denote 
respectively  the  number  of  cyclic  variations  with  d  and  the  nth  root 
of  the  characteristic  equation.  These  hybrid  asyimoetrlc  modes  discus¬ 
sed  above  are  doubly  degenerate  since  an  equally  valid  solution 
results  If  sin  mO  Is  replaced  by  cos  mO  ,  and  cos  niQ  by  -sin  ni0  . 

As  pointed  out  in  the  preceding  section,  no  pure  TE  or  TM  waves 
can  exist  on  an  ellipticeQ.  dielectric  rod.  All  inodes  must  be  hybrid. 

In  contrast  to  the  case  of  a  circular  rod,  the  fields  in  one  of  the 
regions  of  an  elliptical  rod  must  be  represented  by  a  set  of  infinite 
series  of  Mathleu  and  modified  Mathleu  functions;  while  In  the  other 
region  by  a  single  product  term  of  Mathleu  and  modified  Mathleu  func¬ 
tions.  We  shall  be  concerned  only  with  the  modes  hereafter  denoted  as 
the  principal  modes,  which  will  degenerate  to  the  well-known  hybrid 
HE  modes  when  the  eccentricity  of  the  ellipse  Is  zero.  The  principal 

THTl 

♦Physically  speaking^  the  presence  of  E^  in  a  predominantly  H  wave 
(i.e.,  the  HE  wave)  or  vice  versa  (i.e.,  the  EH  wave)  assures  the 
return  path  for  the  electric  or  magnetic  lines  of  force;  in  other 
words,  the  electric  and  magnetic  field  lines  miist  form  closed  loops 
In  the  case  of  the  surface  wave  propagation  along  a  dielectric  rod. 

The  existence  of  a  circularly  symmetric  pure  E  or  H  wave  along  the 
dielectric  rod  Is  a  special  case;  since  the  electric  and  magnetic 
lines  of  force  of  the  E  or  H  wave  have  already  formed  closed  loops. 
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modes  vlll  be  denoted  by  or  depending  upon 

whether  the  modes  are  even  or  odd.  The  axial  magnetic  and  electric 
fields  of  an  even  mode  are  represented  by  eua  even  and  odd  Mathieu 
function  respectively,  and  those  of  an  odd  mode  by  an  odd  and  even 
Mathieu  function  respectively.  The  subscript  m  is  the  order  of 
the  Mathieu  fmction  used  for  the  single  product  term,  and  n  is  the 
nth  root  of  the  characteristic  equation.  The  superscript  1  or  0 
indicates  the  region  Inside  or  outside  in  which  a  single  product  term 
was  used  to  represent  the  field  configuration.  The  symbol  HE  means 
that  the  cross-sectional  field  pattern  of  this  dielectric  rod  mode 
resembles  that  of  eui  H  wave  in  the  metal  guide.  The  symbol  EH  is 
used  if  the  cross-sectional  field  pattern  of  the  dielectric  rod  mode 
is  similar  to  that  of  an  E  wave  in  the  metal  guide. 


2.3  The  Field  Conponents  and  the  Detenninantal  Equations  of  the 

Principal  Modes 

Having  properly  classified  the  modes  we  are  now  in  a  position  to 
describe  the  field  components  of  the  principal  modes  and  to  apply  the 
boundary  conditions  in  order  to  obtain  the  characteristic  equations 
from  which  the  propagation  constants  may  be  determined. 

In  order  to  sin^jlify  the  notations  for  the  Mathieu  and  modified 
Mathieu  functions  without  any  ambiguities,  the  following  abbreviations 
are  used: 


Oeji)  . 

Seje)  = 

= 

Pek^d.lr^l) 

ce*(Ti)  =  ce*(,j,|rQ|) 

Oek^(l)  . 

se*(Ti)  =  se*(Tj,|T^^l)  • 
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The  principal  inodes  asa  separated  Into  foiir  major  types,  the 

mode,  the  mode,  the  mode,  and  the  HE^^^  mode, 

emn  'emn  ’  omn  ’  omn 


(a)  The  HE^^^  Mode 
e  mn 

Appropriate  solutions  of  Maxwell's  equations  for  both  regions  of 

the  elliptical  dielectric  rod  have  been  given  In  section  2.2.  According 

to  the  definition  of  the  HE^^^  mode,  specified  In  section  2.4,  the 

e  mn  ^ 

axial  coiq)onents  of  the  ma^etlc  and  electric  fields  sure  for  region  1 


(0  ^  M  gg) 


iP,z 

H  =  A  Ce  (I)  qe  (ti)  e 
in  u  SI 

lp,z 


(1) 

(2) 


and  for  region  0  (  g  ^  |  <  oo  ) . 


00 


H 


ip  z 

5]  L  Fek  (I)  ce;(n)  e  ° 


z  r  r 

o  r  =  0 


00  ip  z 

"  I!  se*(Tj)  e  ° 

o  r  =  1 


(3) . 

(4) 


A  ,  B  ,  L  ,  and  P  are  the  arbitrary  constants.  The  harmonic  time 
m'  m'  r  r 

dependence  e  has  been  Inylled  In  these  expressions  as  well  as  In 

the  subsequent  expressions  for  theJfleld  Intensities.  Substituting  1 
and  2  Into  2.2-10  through  2.2-13,  and  carrying  out  the  differentiation, 
one  obtains  the  transverse  field  conponents  for  region  1(0^  g  ^g  )# 


H  - 

■  ,.  2  „2 


H  (k^-  P  )p  pmm  m) 

K  =  — ce'(q)  +-giB  Se'(g)  se  (Tj)f  e  (6) 

•  T)j|^  (x  -  p^)p  ®  “  mm=  m'^ 
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E 


13 


{• 


^  A^Ce;(|)ce^(Tj)  +  B^Sej6)8e-(Ti)(  e  (8) 


2  2  2  2  1/2 

where  k^  =  o)  ^l€^  and  p  =  q  (sin  h  5  +  sin  q)  '  .  The  abbreviations 

se^(Tj)  =  ^  have  been  used.  Substituting  equations  3  and  ^ 

into  2.2-10  through  2.2-13,  and  carrying  out  the  differentiation,  one 
obtains  the  transverse  field  conqjonents  for  region  0  (|  ^|<oo) 


H 


13 


00  r  aie  ^  ^ 

Tj  L  V®^r^  PpOek^C  i)  se*{  tj)  J  e 


IS  z 


^o  (kf-S^)pr  =  0 

- -  E  [  V®^r^  1)  Pr°®^rU)  sej(  n)]  e 


H  = 


\  (k^-S^)pr‘^l 


^0  (k?-S^)p 


2.  "  E  [  ?  L^Fek^(  i)  ce*(  tj)  +  P^Gek^C  |)  se*(  tj)  ] 
-6  Jp  r  =  l 


(9) 

IPqZ 

(10) 

iPoZ 

(11) 


E 


13 


(k^-3^)p  r=Q 


E  [  -  ^  1/ek^  ( I )  ce*  ( q)  +  P j.Gek^  ( I )  se*  ( q)l  e 

I"  —  n  ^  ^  J 


iS  z 


(12) 


2  2 

where  k^  =  u)  .  The  prime  denotes  differentiation  with  respect  to 
I  or  q  as  the  case  may  be. 

If  m  is  odd,  r  must  also  be  odd  due  to  the  orthogonality  of  the 
Mathleu  functions.  In  other  words,  idien  m  is  odd  the  series  are 
summed  over  all  odd  values  of  r  ;  and  iidien  m  is  even  the  series  are 
summed  over  edl  even  values  of  r  .  The  problem  with  m  odd  will  be 
ansQyzed  in  detail. 
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Equatlng  the  tangential  electric  and  magnetic  fields  H 

Z  Z  T| 

and  E  at  the  houndanry  surface  i  =  I  *  we  arrive  at  the  following 
tj  *  »o 

equations 

IP  z  oo  16  z 

m  m  0  m  t~> .  r  r  o  r 


r=  1 
odd 

16  z  00  16  z 

mmom  ^.rror 

r  =  1 

odd 


(13) 


(U) 


6,  r  1 

■2^  ^  V®;^5o)®®„(n)]  ® 

c,  -  )  -1- 


(i^-e£) 


p  00 

o 


(k^-  6^)  r 
'  o  o' 

odd 


00  r  / 

,E  ^  L  v®^r^  ^o^  Pr°®^r^  ^o^ ® 


16  z 
^0 


(15) 


1  -  6-.  J  -L 


(1^  -  Pi) 


6  oo 


E  [-I'rF 

.  -  6  )  r=  1  ^o  ' 


16  z 
0 


-O 

O  0 


Odd 


(16) 


Setting  6  =  Pq  =  ?!  eliminating  ce*(n)  and  8e*(n)  by  equa¬ 

tions  13  and  l4,  we  obtain 


GO 


V'.<  «o>  '%<  l)  -  E  ,  «o)  <=°J ( l) 

r  =  1 


odd 

00 


V®m(5o)8e^(T,)  =  E,  ^ 

r  =  1 
odd 


(17) 


(18) 
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k^-  1  cue 

L  P  j 

k?-  P^  OD  0)6 


r-3>  ^i-r  w(‘„)ae.(,) 


k-.  p' 


2  .2 


L  k  “  D  -* 

O  ^ 


2  2 

kf  -  p^  00 


(~ - 5)  y  LFek'(|  )ce*(T|) 

\2  ^2'  i->  ^  P  r  r"*o'  r' 


k;-P  r=l 
odd 


The  method  discussed  in  section  2.3  will  now  he  used  to  eliminate  the 

1)  dependence  in  the  above  equations.  Multiplying  equations  17  and  20 

by  ce  (n)  and  equations  18  and  19  by  se  (tj)  >  and  integrating  with 
s  s 

respect  to  tj  from  0  to  2%  results  in 


ACe(4)€C=  Y  L  Fek  ( |  )  a 
m  m  ^o  ms  m  r  r  ^o  r 

r  =  1 

odd 

00 

B™Se„(5„)  «  S  =  r  P  Gek^(L)6 
n  lu  o  ns  n  ^ «  r  . .  r  o  r 
r  =  X 

odd 


k2  p2 


*-rW‘So)  - 


kf-  p2  oo€,  00 


(4—5)  -4  E  Pr“«‘‘;(5o)^ 

v-2  _  P  r  =  1  r  u  I 


2  o2 


V  ^  e^C^.  B Jl  -  (^)l  se^(l^) 

L  K  -  p  J 


Jj2  .  o2  00 

(-^ - )  Z  L  Fek' (6  )a 

\2  p2  p  r  =  l  V  r''^o^“r8 

*^o  ■  P  odd 


(24) 
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and  8  .  1,3,5,7,9,---  •  C^,  6^.  and  are 

defined  as 


€ 

ms 


C 


m 


S 


m 


rs 


rs 


6 

ms 


ms 


1 

when  m  =  s 

0 

when  m  /  s  ; 

f 

ce^(Ti)dB  j 

ffi 

0 

2« 

/ 

8e^(ti)^B  > 

(25) 

"0 

r 

ce*(i))  ceg(ii)dq 

(26) 

i 

! 

8e*(i))  3eg(Ti)dq 

(27) 

J 

0 

2x 

1 

8e^(n)  ceg(Tj)dTi 

(28) 

r 

ce^(T|)  8eg(Ti)dTj  . 

(29) 

0 


These  Integrals  have  been  evEduated  In  Appendix  A. 4.  For  each  value 
of  m  there  exists  four  infinite  sets  of  linear  algebraic  equations 
which  may  be  combined  to  give  the  following: 


=  L  PekJl  )M  (a  ) 
mmo  m  momrs 


B  Se  ( I  )  =  P  Gek  ( S  )N  O  ) 

m  m'^o  m  m'^o  m'  rs 


(30) 

(31) 


k2-  p2  1 

1  -  (-|— 2) 


coe 

C«m(5o)  Sn^^s^^ms)  ^ -T  = 


k^-  UJ5 


%  "P 


(32) 
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Aj|k^Ce'(i  )-B 
tn  p  ss  o  in 


2  2 


SeJSo^ 

mo  m  rs  ms 


(33) 


M  (a  ),  N  (3  ),  ft  (3  ),  and  R  (a  ,6  )  eure  obtained  according 

m'  rs"  m'  rs"  ^'^rs^  ms"  m'  rs'  ms'  ^ 

to  the  technique  developed  In  section  2.3*  As  Ein  example^  we  choose  the 
m  =  1  mode.  These  constants  sure  then  given  by 


“33  “53 

“35  “55 


“n 

“31 

“51 

“13 

“33 

“53 

“15 

“35 

“55 

(3»^) 


(35) 
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f 


In  -  X  - 

'  11  ^13 


^31  ^51 
^35  ^55 


^33  ^3 
^35  ^5. 


L15 


^31  ^51 
^33  ^53 
^37  ^57 


^33  ^53 
^35  ^55 


... 


s. 


(36) 


Ri(a  ,8,  )  = 

1'  rs'  Is' 


)  «  •  •  «  • 


®ll"  ®13 


«  • 

“31 

“51  **' 

0£^ •  •  • 

\ 

•  •  • 

33 

53 

“37 

ot  •  •  • 

57 

■  +  B.c  - 

•  •  « 

•  •  •  •  •  s 

—  +  .  •  • 

•  4- 

“33  ®53 

-15 

“33  “53 

“35  “55  - 

•  ••  ••• 

“35  “55 

•  •  •  •  •  • 

2x 

2it 

(37) 


where  =  J"  ce^(Ti)dT)  and  ~  J'  ®®i(t)'^T  •  These  infinite  deter- 

mants  may  be^solved  by  the  method  of  successive  approximations  (47). 

Equations  30  through  33  ®re  a  set  of  four  homogeneous  linear 

algebraic  equations  In  coefficients  A  .  B  ,  L  and  P  from. which 

m'  m'  m  m 

only  the  ratios  of  these  coefficients  can  be  determined.  These  ratios 
provide  the  covqpling  factors  between  the  different  coefficients.  For  a 
nontrivial  solution  the  determinant  of  this  set  of  equations  must  vanish. 
Therefore  we  have 
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Putting 

C08h^|^(k^-  (39) 

=  -q^co8h^5^(k^-  p2)  (40) 

after  some  extensive  algebraic  manipulation  equation  38  flnsCLly  yields 

the  transcendental  characteristic  equation  for  the  mode 

^  e  mn 


T  Ce*(|  )  -  Fek'(|  )  1 

1  ^0'  1  m  ^o^  ' 

Ce  ( 5  )  ^  Fek  (5  ) 

*  mo  y  mo 


fl  se;(g  1  %Oek;(gi 


,2  2w  2  ®o  .  2v 

(x  +  y  )(x  7-  +  y  ^ 

+  - r-r - i -  R  (a  ,b  )Q  (P  ,%  )  =  0  .  (41) 

4  4  m'  rs'  ms  ^  rs'  ' 

X  y 

The  propagation  constant  P  of  the  wave  can  be  calculated  from 

e  mn 

the  above  characteristic  equation,  together  with  the  fact  that  P^=  P^ 
or 

+  y^  =  k^q^coBh^|^(^  -  1)  ,  (42) 

o 

in  terms  of  the  frequency,  the  siee  of  the  g\ilde  ,  the  eccentricity  of 
the  guide  and  the  electromagnetic  constsuits  of  both  regions  and,  of 
coxirse,  the  order  of  the  mode. 

The  ratios  between  the  different  arbitrary  constants  are  iii5)ortant 
because  they  indicate  the  co\q)llng  between  the  ajiipllt\ide  and  phase  of 
the  different  field  components.  From  equations  30  through  33^  expres¬ 
sing  all  arbitrary  constants  int^terms  of  ,  we  obtain 


Ce  (I  ) 
m'^o' 


(>•3) 


\|  S  p\/(x^+  y^)(x2+y2  flj  ^ 


m  o 


-  R  (a  ,1 

m'  rs' 


P  Se  (5  )  B 

m _ m  ^  _m 

\  Gek  (5  )«  (P  ) 

The  other  three  typesodf  principeJ.  modes  ,  and 

f  may  he  analyzed  in  a  similar  manner  as  outlined  above;  only 

o  mn 

their  principal  results  will  he  given  here. 


(h)  The  Mode 

e  mn 

The  axial  field  components  for  the  mode  are,  for  region 

e  mn 

1  {0  ^  i^)  , 


OU  xp  X. 

1  r  =  1 


UU  xp  z 

=  S  B?  Se  (I)  se  (tj)  e 
^1  r  =  l  ^  ^ 

odd 


and  for  region  0  ( 5^  -  5  <  oo )  , 


0 
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where  A*,  B*,  L*  and  P*  eure  arBltrary  constants  .  All  transverse 
r  r  m  m 

field  coiQ>onents  can  be  obtained  from  equations  2.2-10  through  2.2-13 
using  equations  46  through  49. 

The  two  determining  equations  for  the  propagation  constant  are 


[i  1 

fi  .  1  ‘o  «cn 

/  2  2v/  o  2v 

(x  +  y  )(x  —  +  y  ) 

1 

inj 

X  y 


R*(a*  ) 

m  rs'  ms'  ’m  rs'  ms' 


=  0 


(50) 


and 

x^  +  y^  =  k^q^cosh^l^C-^  -  l)  ,  (5l) 

o 

2  2 

where  x  and  y  possess  the  same  definition  as  those  In  equations  39 
and  4o  and 

2« 


CO 

II 

ce^(q)  ce*(Tj)dT) 

X  S 

(52) 

(m'^ - — 

II 

DO 

$1^ 

se^(q)  se*(q)dTj 

(53) 

XSb-  / 

ce*U)  se*(Tj)dq 

(54) 

0 

2x 

6*  =  r 

ms  j 

se*'(q)  8e*(ii) 

(55) 

0 


These  integrals  are  evaluated  In  Appendix  A. 4. 

can  be  obtedned  In  the  same  way  as  outlined  for  \i(®i.g>®5jg)^(Pj.g>)^g) 


In  section  2.^a« 
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(c)  The  Mode 

'  '  o  nm 


The  field  coii^>onents  for  region  1  (O  ^  |  ^  |^)  are 


H  =  a  Se  (0  se  (n)  e 
m  in  m 


ipz 


13z 


ce  (tj)  e 

z^  mm  m 


- ^-r  Rse;(0  8e^(Tj)-  V®m(5)ce'(Ti)) 

p(k^-p  )  l““  “  3  mm  m'J 


13z 


13z 


_  = - 'o'  fa  Se  (0  8e'(Ti)+  — b  Ce'(|)ce  (r\)J  e 

\  p(k^-P  J  1“  “  “  ^  “  “  “  / 

(\  f"  ®®m^  +  \Ce^(  Dcejri)  |  e 


!  =  -  - 

p(kJ-3^) 


(56) 

(57) 

(58) 

(59) 

(60) 


!  ,  . . -..Ag- . 

'n,  _/,_2  „2 


p(k^-P") 


The  field  components  for  region  0  ( 6^  -  5  oo )  are 


oo 


=  E  Sj.Gel^r(S)  8e*(T])  e 
o  r  =  1 
odd 


13z 


(62) 


oo 


o  r  =  1 
odd 


(63) 


H  =  — _ 

\  ^/,-2  „2^ 


*0  p(k^3‘^)  r‘='U 


oo  r  flU€  /  I  4Q 

E!  ]  grGek^(Ose*(Ti)-  h^Fek^(  |) 065(11)  |  e 


odd 


(64) 


40  °o  r  ) 

""n  E,  U<^^(fi)se5'(i,).-^h^Fek;(6)ce5(n)  a 


13z 


p(k‘-3‘^)  r=l 
odd 


(65) 
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p(k; 


odd 


’Iq  p(k; 


i|-^  j^-g^  ^k^COsejCn)  +  h^ek^(|)ce*'(Ti)j  (63) 


odd 


^r*  ^r  arbitrary  constants.  Upon  matching  the 

bovindary  conditions,  one  ohteilns  the  dispersion  relations  for  the  propa¬ 
gation  constant 


Se'(g  )  ,  G«4'(l) 

1  m  o  -L  m  0 

”5  Se  C  S  )  Gek  ( £  ) 


and 


,2^  2.  .  2  *=0  2v 

(x  +  y  )(x  —  +  y  ) 

- inr^^ - 

X  y 


,71  )  R  (a  ,6  )  =  0 
rs'lns'  m'  rs'  ms' 


y^  =  coshes  (t^  -  1)  , 

^  o  ^ 


(68) 

(69) 


where  x^,  y^,  O,  t  ^  }  %  >  ^  }  Q_(6  R  (a  /6  )  have 

>  j  I  rs'  rs'  '^ms  ms'  ^  rs  ms'  m  rs  ms' 

been  defined  in  section  2.5A. 


Expressing  sill  arbitrary  constants  in  terms  of  a^  ,  we  have 


V  '  Gekjs,,)  s^(a^3) 


b 

m 


2  g 

(^)G€k'(|  )  N  (6  )  (— )+  Se'(g  ) 

'  2'  m'*o'  m'^^rs'  'a  '  m'^o' 

H _ L. 


m 


in  o 


r^x2  +  y2)(x2  +  y2  |1) -i 


T" 


Sn(Prs^4s) 


(70) 


(71) 


and 


(72) 
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B1 


DO  D 


m  Pek  ( i-)M  (a  )  m 
m  0  m  rs 


where  N  (P  )  and  M  (a  )  have  been  defined  in  section  2.5a  . 
m  rs  m  rs 


(d)  The  Mode 

0  Bn 

The  social  field  con^onents  for  region  1  (O  ^  |  ^  are 


00 

H  =5^  a*  Se  (l)  se  (tj) 
2  r  I  r  ' 

±  r  =  X 


Ipz 


00 

E  =  y  b*  Ce  (s)  ce  (n)  e 
z,  r  r'^'  r'  '' 

X  r  s  X 


Ipz 


(73) 

m 


And  the  axial  field  con^jonents  for  region  0  (|^^|<gd)  are 


0 


13z 


=  h*  Fek  (l)  ce*(Ti)  e 
z  m  m  m 
o 


l^z 


(75) 

(76) 


aj^  b*,  ^  and  h*  are  the  arbltreu:y  constants. 

The  dispersion  relations  for  the  propagation  constant  are 


ri 

1 

fi 

,  e  Fek'(|  ) -j 
1  0  ni'*o' 

y2  8el:jgj 

*  2  €,  Fek  (1  ) 

y  1  m'^o^J 

/  2  2.,  2  o  2v 

(x  +  y  )(x  ^  +  y  ) 

^1 

- im - 

X  y 


and 


+  y^  =  kf  cosh^|^(^  -  i)  . 


(78) 


2  2 


Where  x  ,  y  ,  a*^, 


given  In  section  2.^'b. 


Some  Interesting  observations  concerning  these  dispersion  relations 
for  the  principal  modes  can  be  made: 

The  presence  of  the  factor  R(a  ,6  iX)  or 

^  m'  rs'  ms' 

R*(ci^  )  In  the  dispersion  relations  Is  the 

m'  rs^  ms'  ^  rs'  ms' 

result  of  using  an  Infinite  series  to  represent  the  field  com¬ 
ponents  In  one  of  the  two  regions.  According  to  Karbovlak's 
approximations  (27),  l.e.,  oe^(T))  =  ce*(Tj),  =  se*(Tj), 

ce'(Ti)  =  -se  (n)m,  eind  se'(n)  =  ce  (n)m,  one  obtains  for  these 
m  m  m  m 

principal  modes  the  dispersion  relations  that  are  the  same  as 

those  given  above  except  that  the  factor 

or  R?(o*  replaced  by  -m^  . 

m  rs  ms  ni  rs  ms 

The  differences  between  the  dispersion  relations  for  the 

mode  and  the  mode  or  the  mode  and  the 

e  mn  e  mn  o  mi 

HE^^^  mode  are  the  factors  R  (a  ,6  )  ft  (P  )  and 

0  mn  m  rs'  ms  ^  rs'  ms 

*  N'-™®*^ically  spealcing,  these  factors 

are  In  general  not  identical  except  when  =  oo  , 

It  Is  therefore  expected  that  these  modes  are  not  degenerate 

modes,  In  general.  However,  It  was  found  numerically  in  the 

next  chapter  that  ■^en  m  =  1  and  n  =  1  within  a  certain 

arbitrary  region  In  the  x-y  plane*  ^l^^re'^s^  “* 

“l(°?e'®la)  «l(“?e’ne>  '  Thus  the  ^HE^^  mode  and- the 

mode,  or  the  He(^^  mode  and  the  mode  may  be  considered 

’  o  11  O  Xi. 

degenerate  within  this  region.  One  may  generalize  the  above 

*Here  x  and  y  are  not  the  coxnponents  In  the  rectangular  coordinates, 
x^  Euid  y2  are  defined  by  equations  2. 5a- 39  siKi  2.5a-40  respectively. 


0 


X 


Fig.  II-3*  Tbe  mode  and  the  mode,  or  the 

e  fflQ  /-\  e  oiQ  0  OQ 

mode  emd  the  HE^  ’  mode  are  almost  degenerate 
o  mn 

vlthln  the  shaded  region.  Tbe  houndeury  of  the  shaded 
region  Is  quite  eurbltrary.  It  Is  determined  accord¬ 
ing  to  the  allowahie  differences  of 

and  R*(<i*  ,6*  )Q*(p»  )  • 

m'  rs'  ms'^''^rs''  ms' 
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observatlon  by  stating  that  In  the  x-y  plane  there  exists  a  certain 

region  (shaded  In  Figure  II-2)  in  which  the  and  the  node 

e  mn  e  mn 

or  the  HE^^^  and  the  HE^^^  node  may  be  considered  degenerate.  The 
o  mn  o  mn  '  ** 

boundary  of  this  certain  region  is  determined  according  to  the  allowable 

differences  between  the  factor  R  (a  ,6  )  ft  (P  )  and  the  factor 

m'  rs'  ms' 

R*(a*  ,6*  )  ft*(3*  *y*  )  . 
m'  rs'  ms'  '^'^rs“ms' 


2.6  Cutoff  Frequencies  of  the  Principal  Dominant  Modes 

It  is  known  that  x  and  y  are  the  roots  of  the  dispersion  rela¬ 
tions.  Combining  equations  2.5-39  and  2.5-^  we  arrive  at  the  propeigation 
constant 


,1/2 


P  = 


q  cosh  I 


q  cosh  I 


2  €1  2  1/2 
r  X  +  ^  y  -1  ' 

ii.i 


(1) 


2  2  2 

In  order  to  have  a  guided  wave,  P  ,  x  and  y  must  all  be  real 

2 

and  positive*.  One  recalls  that  the  positive  and  real  values  of  y 

indicate  that  the  field  intensities  outside  the  dielectric  rod  decay 

2 

with  increasing  distance  from  the  surface  of  the  guide.  If  y  is 
negative  and  reed,  the  expressions  for  the  field  components  will  indicate 
the  presence  of  an  outgoing  radled  wave  at  a  large  distance  from  the  sur¬ 
face  of  the  dielectric  rod,  which  can  only  come  from  an  infinitely  long 


— - - 2 - - - 

■*The  fact  that  x  and  y  must  all  be  reed  and  positive  offers  a  way 
to  determine  the  upper  and  lower  bo\aids  of  the  propagation  constemt  p. 
According  to  equations  2.5-39  and  2.5-^^  P^  -  and  P  ■  • 

Thus  k^  6  p  ^  . 
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p 


(in  the  z  direction)  line  type  source  located  at  some  finite  5  . 

Such  sources  have  not  been  postulated  in  the  assumptions.  In  fact, 

2 

the  concern  here  is  with  tlie  source-free  problem.  Thus  y  must  be 

positive  real  for  all  surface  guided  waves  and  consequently  the  lowest 

2 

permissible  value  of  y  is  zero.  The  propagation  constant  and  the 

2 

frequency  corresponding  to  this  value  of  y  are 


and 


(y  =0) 


q  cosh  i 


(•ji-i) 

o 


(2) 


(y2=o) 


(3) 


q  cosh  Iq  -  1)/"*. 


respectively,  x  corresponds  to  the  root  of  the  characteristic  equation 
2 

with  y  =  0  .  The  frequency  defined  by  equation  3  is  called  the  cutoff 
frequency  of  the  wave,  since  below  such  frequency  the  mode  can  no  longer 
exist  on  tne  dielectric  guide.  Physically  it  means  that  below  this 
cutoff  frequency  the  structure  can  no  longer  support  such  a  wave  and 
thereby  ceases  to  be  a  binding  medium. 

The  approximate  expressions  of  the  modified  Mathieu  functions  for 
small  X  and  y  are  derived  in  Appendix  A. 2.  For  small  valiies  of  y 
we  have 


Fek' ( 6  ) 

1  m  o  _  ^  _ 

^o^  y^  L  8(m^-l)  cosh^l 


-41 


-m 


®  [(m+l)  + (m-l)e  °]  +0(y^) I  , 


W 


[for  m  ^  3  odd) ] 
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Fek'(|  )  ,  r  2 


11  -  O  Qt  ’o  -2|  ), 

1  I,,  ye  .^,eve  (jj 


^FekM°i  a  in(  °  y  °  ~  )[3-ae  '°]40(y‘*) 

y  l'®o'  y  I  ^  cosh  5-  ^  cosh  |  J 


Eind 


,  Gek’d  ) 
1  m  o 

^  Gek^j(§o) 


.2 


2| 


-m  - . 


8(m  -l)  cosh  I 


-4| 

2 — [(m+1)  +  (m-l)e  °]+0(y^)(, 


[for  m  ^  3  (ni  odd)]  (6) 


^  Gek'(5  )  T  f  a  ^o  -2f  y  '} 

^  ^  - )  [3+2e  °]  +0(y^)  f 

y  Gsk^^dg)  y  [  Bcosh^l  4  cosh  6  J 


(7) 


where  a  Is  the  Euler's  constant.  For  small  values  of  x  ,  we  have 


1 


tanh  0(x^)]  ,  [for  mil  (m  odd)] 


(3) 


and 


x^  Sejjj(  Iq) 


Tcoth  Gg  +  0(x^)j  ,  [for  m  i  1  (m  odd)]  (9) 
x 


where 


m-1 

2 


E  (-1)" 


r=0 


(^  +  r) : 

-  r)l 


(cosh  1^)^^ 

/D-l  ^ 

(-T 

4  ^  (2r)l 


m-1 

2 


,  +  r) ! 

E. 


(cosh  1^) 


2r 


r=0 


(V  -  y)i 


-  r) 

4  (2r+l): 


and 


m  . 
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^2  = 


m-1 


"T  +  r):  (slnh  i  ) 

E(-if  " 

r=0 


2r 


m-1 

2 


(^-r)l  (^-r) 

_ k  {2r)l 


+  r)  1  (sinh  | 

7  (.if  1.2 _ ^  ^ 

ho  r):  i^-r) 


(2r+l) 1 


It  can  be  shovn  that  for  small  values  of  x  and  y  to  the  first 
order  approximation, 


R  (a  ,6  ) 

m  rs  ms 


)  «  R*(a*  ,6*  )  Q*(p*  J*  )  «  -m" 
rs'  ms  m  rs  ms  ^  rs  ms 


(10) 


Substituting  the  above 'approximations  into  the  characteristic 
equation  2.5-41  or  2.5-50,  one  obtains  for  the  even  principal  waves 


2 

X  » 


[m(l  ^o  °1  ^o  *^2^^ 


2^0  -4| 

^  [(m+1)  +  (m-1)  e  °] 


(11) 


:osh  § 


[for  m  ^  3  (m  odd)] 


and 

€  e 

[  ( 1  +  -^)  +  ( —  tanh  I  +  coth  6  )  ]  8 
2  ^1  ^1  °  ° 

e  ^^o  -25  ~ 

S-^in(e«  - )[3-2e  °] 

1  cosh  I  2  cosh 

o  o 

[for  m  =  1] . 


(12) 


Upon  inspection  of  equation  11  we  may  immediately  conclude  that 
the  right  hand  side  of  the  equation  is  always  positive  and  non- zero  and 
is  not  necessarily  small  for  all  values  of  and  e^/ ,  thvis  x 

is  not  zero  and  is  not  necessarily  small.  In  other  words,  the  imposed 


small  X  approximation  is  net  valid  and  x  must  be  determined  from 
t!ie  original  characteristic  equation  2. 5-^1  or  2.5-50  with  y  -  0  . 

Tne  same  ocnclusion  may  be  reached  for  m  ^  2  (m  even)  even  princi¬ 
pal  modes . 

From  equation  12  it  is  noted  that  as  y  approaches  zero 

a  , 

to(e  2^cQsh  £  ^  approaches  -oo  ,  thus  the  right  hand  side  of  this 
equation  approaches  -t-0  .  In  other  words,  as  y  approaches  zero,  x 
also  approaches  zero  and  the  imposed  small  x  approximation  is  valid. 
Therefore  the  cutoff  frequency  of  the  mode  or  the  mode  is 

zero  (refer  to  equation  3) •  Ib  should  further  be  noted  that  the 
mode  and  the  mode  are  degenerate  at  zero  frequency. 

Substituting  the  above  approximate  expressions  4,  5;  6;  "7,  8,  9 
and  10  into  the  characteristic  equation  2.5-68  or  2.5-77,  we  arrive  at 
the  following  expressions  for  the  odd  principal  waves: 

€_  C 

-  coth  Gg) ] 

(li) 


o  €  e 

•4{m''-l)  rm(l  +  — )  +  (tanh  ~  ^2 


- 2 —  t(®+l)  +  (“-!)  e  ] 

cosh  r 


[ for  m  ^  3  (m  odd) ] , 


and 


2 

x 


o[(l  +  •^)  +  (tanh  l_  +  coth  E  )  ] 
1  0^0 


2L 


0  e 


e-,  ^2. 

1  cosh  6 


£n(e' 


^o  -2i 

“iL® - )[3-2e  ° 

cosh  L 


(for  m  =  l) 


(14) 


Similar  conclusions  as  those  for  the  even  principal  waves  are 
reached.  For  the  m  ^  3  (i“  odd)  odd  principal  waves,  the  right  hand 
side  of  equation  13  is  always  positive  and  non-zero,  thus  x  is  also 
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positlve  and  non-zero.  It  csin  be  shown  that  the  same  conclusion 
applies  for  the  m  ^  2  (m  even)  odd  principal  waves.  However,  for 
the  m  -  1  odd  principal  wave,  according  to  equation  l4,  as  y 
approaches  zero  x  must  also  approach  zero.  There  exists  no  cutoff 
frequency  for  the  or  mode  (refer  to  equation  i) .  And 

at  zero  frequency  the  and  modes  are  degenerate. 

The  results  of  the  analysis  in  this  section  are  summarized  as 
follows . 


i.  Along  an  elliptical  dielectric  rod  there  are  only  two  non¬ 
degenerate  modes,  namely  the  mode  and  the 

^  e  11  o  11 

mode,  which  possess  no  cutoff  frequencies. 

ii.  It  can  be  observed  from  equations  12  and  14  that  as  the 

elliptical  cross  section  of  the  dielectric  rod  gets  flat¬ 
ter,  X  approaches  zero  more  slowly  since  coth  is 
very  large  if  is  very  small.  This  fact  has  been 

verified  in  the  next  chapter  (see  Figures  III-l  and  TII-6) . 

iii. The  cutoff  frequencies  of  all  the  other  modes  are  higher 
for  flatter  elliptical  cross  section  rod. 


2.7  Transition  to  Circular  Cross  Section 

As  an  ellipse  degenerates  to  a  circle  its  semifocal  length  q 
tends  to  zero  while  |  approaches  infinity  so  that  the  product, 

q  cosh  I  or  q  sinh  |  or  qe  tends  to  a  constant  r  which  is 

the  radius  of  the  degenerated  circle.  The  degenerate  forms  of  the 
Matiiieu  and  modified  Mathieu  functions  are  given  in  Appendix  A.fi. 

Using  these  degenerate  expressions  one  obtains  the  following  de¬ 
generate  forms  for  the  factors  appearing  in  the  characteristic  equations: 
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in  o 

Ce  (i  ) 
m  *0' 


X  J' (x) 
a 


(1) 


m  o 

Se  (5  ) 
m' 


X  J'(x) 
a 

J  (x) 
a 


(2) 


Fek'(5  )  y  K' (y) 
m  o  p 

Fek  ( t  )“  "  K  (y) 
m  o  m 


(3) 


- 

Gek_(  ij 


B  O 


a  »  Of* 
rs  rs 


6  «  6* 
ms  ms 


R  (a  iY-  ) 

m'  rs'  ms' 


m'" 

K„(y) 

i 

9  * 

P*  -*  1 

when 

r 

=  s 

(5) 

rs 

rs 

-  0 

when 

r 

/  s 

■^ms" 

.p  -* 

ms 

m  when 

m 

=  s 

(6) 

0  when 

m 

/  S 

'rs'^s^ 

«  R*(a*  ,%*  )  ^(p*  yi-*  )  -  -m^ 
m  rs^'nis'  ^  rs  'ms' 

(7) 

2222  2222  » 
with  X  =  r^(k^  -  P  )  emd  y  *  r^(P  "  ^q)  •  Putting  these  degenerate 

expressions  into  amy  one  of  the  characteristic  equations  for  the  four 

types  of  principal  inodes  yields 


1  *0 1  1 

X  j^(x)  y  K,(y)  J 

=  m 


/  2  2w  2  S  2n 
(x  +  y  )(x  _g_  +  y  ) 

ITU 

X  y 


(8) 


which  is  the  characteristic  eqxiation  for  the  hybrid  HE^  mode  on  a  circu¬ 
lar  dielectric  rod.  (Compare  with  equation  III-40  in  reference  (13))* 
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One  can  easily  show  that  only  one  term  of  the  infinite  product 
series  used  to  represent  the  field  configurations  in  one  of  the  two 
regions  remains.  The  order  of  this  remaining  term  corresponds  to  the 
order  of  the  single  product  term  representing  the  field  configurations 
in  the  other  region.  In  other  words. 


L 

m 

when 

r  =  in 

(9) 

-♦ 

0 

when 

r  /  m  , 

p 

r 

- 

P 

m 

when 

r  =  m 

(10) 

0 

when 

r  /  m  , 

A* 

r 

- 

A* 

m 

when 

r  =  m 

(11) 

0 

when 

r  /  m  , 

B* 

m 

when 

r  =  m 

(12) 

0 

when 

r  ^  m  , 

®r 

®m 

when 

r  =  m 

(13) 

0 

when 

r  ^  m  , 

h 

r 

- 

h 

m 

when 

r  =  m 

(1^) 

- 

0 

when 

r  ^  m  , 

a* 

r 

- 

a* 

m 

when 

r  =  m 

(15) 

- 

0 

when 

r  /  m  , 

h* 

r 

-* 

b* 

m 

when 

r  =  m 

(16) 

0 

when 

r  ^  m  . 

It  thus  appears  that  the  and  HE^®^  modes  are 

e  mn  e  mn  o  mn  o  mn 

degenerate  when  the  elliptical  cross  section  degenerates  to  the 


clrcvilar  cross  section 
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CHAPTER  III  -  NUMERICAL  ANALYSIS  OF  THE  DOMINANT  MODES 

It  is  the  piirpose  of  this  chapter  to  investigate  in  detail  the 
propagation  characteristics  of  the  dominant  modes  on  a  lossless  ellip¬ 
tical  dielectric  rod. 

After  a  brief  review  of  the  method  for  computing  the  numerical 
values  of  the  Mathieu  and  modified  Mathieu  functions,  the  transcenden¬ 
tal  characteristic  equations  derived  in  the  previous  chapter  for  the 

HE^^^  mode  and  the  mode  are  solved.  Several  graphs  showing  how 

e  11  o  11  o  JT  a 

the  propagation  constants  vary  with  parameters  are  given.  Interpreta¬ 
tions  of  the  results  are  given.  The  field  configurations  and  the  axial 
electric  field  extent  of  these  waves  are  also  considered. 

3.1  Computation  of  the  Mathieu  and  Modified  Mathieu  Functions 

It  is  knovjn  that  the  periodic  Mathieu  functions  may  he  expanded  in 
terms  of  em  Infinite  series  of  trigonometric  functions,  and  that  the  cor¬ 
responding  modified  Mathieu  functions  can  he  expanded  in  terms  of  an 
infinite  series  of  products  of  Bessel  functions  (see  Appendix  A.l) . 

These  Bessel  function  product  series  converge  very  rapidly  [see 
McLachlan  (45),  p.257].  As  has  been  pointed  out  on  page  21,  Chapter  II, 
in  order  that  the  solutions  of  the  Mathieu  differential  equation  be 
periodic,  the  characteristic  number  c  or  the  sepeuration  constant  of 

the  wave  eqvtation  mvist  satisfy  a  certain  transcendental  infinite  con- 

2  * 

tlnued-fraction  equation  which  is  a  function  of  y  .  Furthermore,  the 

2 

coefficients  of  these  Infinite  series  are  functions  of  y  and  c  (42, 

45). 

*The  infinite  continued-fraction  equation  was  first  used  by  Ince  (48)  in 
calculating  the  characteristic  niuibers;  c  and  y^  are  defined  in  equa¬ 
tions  2.2-18  and  2.2-19. 
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Supposing  one  Is  Interested  In  obtaining  the  numerical  value  of  a 
certain  modified  Mathieu  fvinction  of  order  m  ,  he  must  first  determine 
the  val\ie  of  the  characteristic  number  which  is  the  root  of  an  infinite 

continued- fraction  transcendental  equation  and  then  find  the  coefficients 

2 

from  the  three-term  recurrence  relations  which  are  functions  of  Y"  and 
c  .  Substituting  these  coefficients  into  the  infinite  Bessel  fionction 
product  series  and  carrying  out  the  computations,  he  then  finally  obtains 
the  result. 

According  to  the  above  description,  it  is  qvdte  evident  that  the  task 
of  computing  the  numerical  values  for  a  great  numioer  of  Mathieu  and  modi¬ 
fied  Mathieu  fionctions  is  very  time  consuming  and  laborious.  Fortunately 

it  is  found  that  the  characteristic  numbers  euad  the  coefficients  for  a 

2 

certain  finite  range  of  y  ,  which  is  the  remge  of  Interest  for  this 
present  problem,  have  been  tabulated  and  published  by  the  National  Applied 
Mathematics  Laboratories  of  the  National  Bureau  of  Standards  (49).  These 
tabulated  values  are  used  in  our  coniputations . 


3.2  Solutions  of  the  Characteristic  Equations 

•  (1) 

The  solutions  of  the  characteristic  equations  for  the  dominant 

mode  and  the  mode  will  now  be  considered.  It  can  be  seen  that  all 

o  11 

these  transcendental  characteristic  equations  are  of  the  form 

y>  x)  =  0  .  (1) 

o  €3_ 

Knowing  |  which  determines  the  eccentricity  of  the  elliptical  cross- 
o 

section  and  e  /e,  which  is  the  relative  dielectric  constant  of  the  sur- 
o'  1 

roxindlng  medium  and  the  medium  of  the  rod,  equation  1  reduces  to 


g(y^x) 

6  /€^  = 
o'  1 


5o= 


const . 
const . 


0 


(2) 
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y  and  x  are  related  to  the  saajor  axis  of  the  rod,  the  frequency, 
the  propagation  constant  and  the  characteristics  of  the  medium  by  the 
relations 


and 


2  2  .2,  ,.2  ^2. 

y  =  -  q  cosh  P  ) 


2  2  ^2,  2  ^2v 

X  =  q  cosh  -p  )  , 


respectively.  In  order  to  have  propagating  waves  on  the  dielectric  rod 
X  and  y  m\ist  both  be  positive  and  real.  Furthermore,  for  these 
dominant  modes  as  y  varies  from  0  to  +00,  x  varies  from  0  to 
some  finite  positive  constant  which  is  a  function  of  €  /e  and  6  . 

0X0 

Eqxiation  2  can  most  readily  be  solved  by  the  "cut  and  try"  method. 
As8\aming  y  to  be  some  finite  consteuat,  say  y^  ,  the  first  root  of  x 


can  be  found  by  plotting  the  function  g 


(y  ,x)  versus  x 


Ipsconst . 

SQ/ej^sconst. 

as  X  varies  from  zero  and  up,  euad  obtaining  the  first  value  x^  where 
the  function  is  equal  to  zero.  Then  by  setting  y  to  be  another  cons¬ 
tant,  the  above  process  is  repeated. 

The  above  method  of  solution  will  now  be  applied  to  the  character¬ 
istic  equations  for  the  even  and  odd  dominant  modes. 


(a)  The  Even  Dominant  Mode,  the  HE^^^  Mode. 


Prom  equation  2.5a-4l  setting  m  =  1  and  n  =  1,  the  characteristic 

.(1) 


equation  for  the  mode  is 
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[7  ^  7  '1  oeki(lo)  . 


,2  2w  2  o  2. 
(x  +  y  )(x  —  +  y  ) 


^  VSrs'V  -°  •  (3) 


RT(a  ,6  )  and  Q,  (6  ^X,  )  are  given  by  (see  equations  2.5-36  and 

X  ITS  XS  X  ITS  X6 

2.5-37  ) 

R,  (a  ,8^  )  = 

1'  rs^  Is' 


and 


(4) 


where 
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\g  =  j  8ej^(Ti)  ce^(Tj)dti  , 
0 


23t 

“rs  "  /  ceg(n)dTi  , 

0 


2jt 

0 


T^lg  =  ce|(Ti)  Beg(Tj)dTi 

Q 


f 

^rs  "  J  seg(Ti)dn 

0 

2« 

0 


(6) 

(7) 

(8) 

(9) 

(10) 

(n) 


Upon  examining  the  solutions  of  equations  6  through  11,  the  foUoving 
deductions  can  be  made: 


(i) 

(11) 

(111) 


1^11 1  >  1^3!  > 
!^l>  1^13!  > 

Kii  >  1^131  > 

I«i3l 

|a  I  >  or  > 
”  l«3ll 

I  I 

|a  I  >  or  > 
l«5ll 


1^5!  ""  >  l^sl 

l^l5l>->  l^lsl 


a 


35 1 


>  a, 


3s  I 


l“35l 
or  > 

l“53l 


\^l\  >  >  l%l 


(iv)  IP^I  >  >  ...  >  16^J 

1^13! 

I  ^23 1  >  or  >1 P35 1  >  ’  *  ’  ^  1 P3S 1 

l^3ll 


r  =  1,3,5) 
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1^15!  1^35! 


It  is  found  that  three  terms  (i.e.,  s  =  1,3,5  and 
of  the  expansion  are  sufficient  to  approximate  the  value 
Ql(Prs,^s)  In  the  region  of  interest  (0^x^5  ,0-y-3)'  The 
validity  of  the  "3-term"  approximation  can  best  be  illustrated  by  the 
numerical  values  in  the  following  table: 


lo  ^  y 


3.0 

1.1 

2.1 

0,1 

2.9 

0.7 

1.4 

2.7 

0.1 

2.9 

0.2 

2.8 

4.7 

0.1 

1.7 

Exact* 


-0.9999 

-0.9999 


-0.99-99 

-0.8979 

-0.9991 

-0.7533 


3-Term  Approx. 

-0.9999 

-0.9999 

-0,9999 

-0.8984 

-1.0002 

-0.7781 


f  =  3.0  corresponds  to  an  ellipse  with  eccentricity  =  O.0932  (or  the 
*0 

ratio  of  semi-minor  to  semi -major  =  O.995)  which  is  very  close  to  a 

circle.  The  value  of  Ri(arg;\s)<ii(Pr3^^1g)  ®  circular  rod  is 

exactly  1.00000  for  all  values  of  x  and  y  .  As  gets  smaller, 

l.e.,  as  the' ellipse  gets  flatter  ,  the  dependence  of 

0  (0  .X  )  with  X  and  y  becomes  more  pronounced,  and  for  large 
^1  rs'  Is 

*By  the  exact  value  we  mean  the  value  obtained  by  using  five  terms, 
(i.e.,  s  =  1,3,5, 7.9  and  r  =  1,3,5, 7, 9*)  of  the  expansions  (equa¬ 
tions  6  through  U) . 
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valvies  of  X  or  y  the  value  of  R.  (a  ;6  )Q^(P  )  gets  quite 

small.  It  means  that  more  terms  of  the  expansions  (equations  6  through 
11)  are  required  to  obtain  the  value  of 

the  field  components  of  the  mode,  as  the  ellipse  becomes  flatter  and 
approaches  a  strip.  Conseqiiently,  the  conqputations  become  more  Involved. 
However,  If  we  limit  our  interest  to  only  the  region  oo  $  0.2  , 

0  ^  X  ^  5  and  0  -  Y  -  i  >  the  "three-term  approximation"  is  sufficient 
for  our  computations. 

According  to  the  discussion  in  section  2.5  it  is  known  that  the 
e^ll^  and  modes  cannot  in  general  be  considered  degenerate 

except  in  a  certain  region.  It  can  be  seen  from  the  following  table 
that  our  region  of  interest  is  within  this  region. 


^o 

X 

y 

3-Tenii  Approx. 

3-Term  Approx. 

3.0 

1.1 

0.1 

-0.9999 

-0.9999 

2.1 

2.9 

-0.9999 

-0.9999 

0.7 

1.4 

0.1 

-0.9999 

-0.9999 

2.7 

2.9 

-0.8984 

-0.8990 

0.2 

2.8 

0.1 

-1.0002 

-1.0004 

4.7 

1.7 

-0.7781 

-0.7805 

Since  in  this  region  Ri(a^s'^ls^^l^^rs'^ls^  ^  ^l^^Js^^ls^'^l^^s^^ls^^ 

one  may  assume  that  the  and  He!?^  modes  are  almost  degenerate 

e  11  e  11 

and  only  one  of  the  modes,  the  mode,  needs  to  be  considered. 

Now  we  eu:e  in  a  position  to  compute  the  roots  of  the  characteristic 
equation  3  for  the  mode. 
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The  con^>utations  were  carried  out  on  a  high  speed  electronic  com¬ 
puter,  the  IBM  7090*.  The  coefficients  and  prepared  hy 

r  r 

NBS  (49)  were  stored  In  the  con^iuter's  memory  cells.  A  three-point 

Lagreuiglan  Interpolation  (30)  sub-routine  was  used  to  Interpolate  the 

coefficients  and  B^™^  from  the  stored  values.  The  nvimber  of 

r  r 

decimals  for  the  various  coefficients  obtainable  with  a  meuclmum  error 
of  2.5  uni' 3  In  the  last  place  by  this  Interpolation  method  have  been 
tabulated  In  the  NBS  Table.  It  was  found  that  the  values  of  the  Mathleu 
fvinctlons  or  the  modified  Mathleu  functions  obtained  using  these  Inter¬ 
polated  coefficients  were  correct  at  least  to  the  third  significant 
figure.  The  roots  of  the  characteristic  equation  were  found  according 
to  the  method  outlined  on  page  60  of  this  chapter.  The  results  are 
shown  In  Figure  III-l  for  the  case  of  ~  various 

values  of  ranging  from  8^  .  3.0  to  -  0.£  . 

It  Is  required  that  the  propagation  constant  inside  the  rod  be  the 
same  as  that  outside  the  rod,  l.e.,  from  equations  2.5-39  and  2.5-IO  , 


2 

X 


2 

+  y 


2q  cosh  I  2  -  e 

>  ( — 7 — -)  « (i^  - 1)  • 

o  o 


(12) 


The  intersection  of  the  function  in  Figure  III-l  with  the  circle, 
determined  by  eqviation  12,  gives  the  values  of  x  and  y  required. 
The  propagation  constant  6  of  the  wave  Is  related  to  x  emd  y  by 


and 


2q  cosh  5  2  „ 

-  ( — ^ — -)  « 


—  -  (— ) 

,  €  ''  2«^ 
L  0 


y  = 


2q  cosh  12 

( - ; - 2) 


X  p  2 

(-2-)  -  1 
^  2rt' 


(13) 

(1*^) 


■••The  facilities  of  IBM  7090  were  provided  by  Western  Data  Processing 
Center  at  UCLA. 


Pig.  III-l.  Roots  of  the  characteristic  equation  for  ■ 

Minor  axls/MaJor  axis  =  tanh  i 
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Fig.  III-2.  Normalized  guide  wavelength  \/\^  of  the 
as  a  function  of  normalized  major  eucls. 


mode 

e  11 
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respectlvely,  where  is  the  free-space  wavelength.  Carrying  out 

the  method  outlined  above  graphically,  the  results  are  given  in  Figure 
III -2.  The  normalized  guide  wavelength  plotted  against  the 

normalized  major  axis  2q  cosh  £  /  X  for  various  values  of  4  in 
Figure  III-2.  It  is  noted  that  the  guide  wavelength  X  is  related  to 
the  propagation  constant  0  by  the  relation  P  =  2jt/x  and  2q  cosh  5^ 
is  the  major  axis  of  the  ellipse.  As  expected,  no  cutoff  frequency 
exists  for  this  dominant  mode.  For  small  values  of  2q  cosh  6  /X 

(i.e.,  the  size  of  the  major  axis  2q  cosh  4^  compared  with  the  free 
space  wavelength  X^  is  small)  the  guide  wavelength  approaches  that  of 
the  free  space  wavelength;  for  large  values  of  2q  cosh  i't 

approaches  asymptotically  to  the  characteristic  wavelength  of  the  rod 


material,  . 


For  small  valuies  of  y  which  correspond  to 


the  small  values  of  2q  cosh  I  / X  the  modified  Mathieu  functions 
describing  the  field  outside  the  dielectric  rod  decay  very  slowly; 
physically  it  means  that  the  field  strength  of  the  wave  falls  off  very 
slowly  away  from  the  rod  and  only  a  small  part  of  the  energy  is  trans¬ 
ported  within  the  dielectric  cylinder.  The  guide  wavelength  of  this 
hybrid  mode  actually  becomes  that  of  a  transverse  electromag¬ 

netic  plane  wave*  as  the  size  of  the  dielectric  rod  becomes  vanishingly 
small.  For  very  large  values  of  y  which  correspond  to  very  large 
values  of  2q  cosh  1^/^^  ^  the  modified  Mathieu  functions  describing 
the  field  outside  the  dielectric  rod  disappear  very  quickly,  so  the  field 
strength  of  the  wave  outside  the  rod  vanishes  very  fast  and  almost  all 


*Although  the  wave  is  propagating  at  the  plane-wave  velocity  of  a  medium 
it  does  not  follow  that  the  wave  is  entirely  transverse.  See  reference 
(51). 
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the  energy  is  transported  within  the  dielectric  cylinder.  The  guide 
wavelength  of  the  hybrid  mode  approaches  that  of  a  TEM  pleuae 

6  lx 

wave  propagating  In  a  uniform  medium  filled  with  a  dielectric  of 
dielectric  consteint  .  The  above  discussion  concerning  the  field 
decay  properties  of  the  wave  will  be  substantiated  later  In  this  chap¬ 
ter  with  numerical  results. 

It  may  be  further  observed  that  for  a  fixed  valxie  of  2q  cosh  |  /x 

^  o'  o 

as  the  ellipse  becomes  flatter,  i.e.,  as  becomes  smaller,  the 

guide  wavelength  become.s  closer  to  the  free  space  wavelength.  This 

effect  can  best  be  iilustrated  by  Figure  III-3  in  which 

plotted  against  for  various  fixed  values  of  2q  cosh 

The  fact  that  the  variation  of  the  curve  becomes  gentler  as 

2q  cosh  gets  smaller  is  expected,  since  at  very  low  frequencies 

most  of  the  energy  is  outside  the  dielectric  rod  thus  the  geometry  of 

the  cross-section  is  not  important  as  far  as  the  guide  wavelength  is 

concerned. 


It  is  also  noted  that  for  a  fixed  value  of  2q  cosh  |  /x  there 
is  more  binding  dielectric  material  in  a  circular  rod  ( =  oo)  than 
in  a  flatter  elliptical  rod,  therefore,  (x/x^)  is  smaller  for  larger 


I  .  However,  this  is  not  the  only  reason.  Supposing  we  plotted 

2q  cosh  |o\2 

X/X  against  the  normalized  cross-sectional  area,  ( - - )  temh  | 

O  Q  o 

for  various  fixed  values  of  in  Figure  III-4.  It  can  be  seen  for 

2q  cosh  I  2 

very  small  values  of  ( - )  tanh  |  ,  say  <0.05  >  that  x/\ 

X«  o  o 

°  ,2q  cosh  In  >2 

for  all  values  of  |  .  As  ( - r - -)  tanh  |  gets  larger,  the 

o  ^o  ° 


effect  of  varying  |^  becomes  more  noticeable.  For  a  fixed  value  of 
2q  cosh  lo  2 

( - )  tanh  I  ,  x/X  i®  smaller  for  smaller  |  .  This  behavior 

'  X  '  '  ^0  '  o  o 


Fig.  IJI-3.  Normalized  giilde  wavelength  of  the  mode  as  a  function  of 


NCSA 


2q  cosh  I  2 


(- 


-)  tanh  I 


k  area 


Fig.  III-4.  Normalized  guide  wavelength  of  the 

mode  as  a  function  of  normalized 

e  11 

cross'sectlonstl  area. 
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suggests  that  the  field  intensity  is  more  concentrated  in  a  circular 
rod  than  in  an  elliptical  rod  with  the  same  cross-sectional  area,  and 
that  more  energy  is  transmitted  inside  the  circular  rod.  We  conclude 

that  the  circular  dielectric  rod  is  a  better  binding  medium  for  the 

(1)  /2q  cosh  lo  2 

HE  mode  than  an  elliptical  rod.  As  ( - r - )  tanh  g  becomes 

6  i..L  A.Q  O 

very  large,  the  effect  of  varying  on  again  becomes  quite 

small,  since  most  of  the  energy  is  carried  inside  the  dielectric  rod; 

therefore,  the  geometry  of  the  cross  section  is  not  important.  When 
2q  cosh  §Q  2 

I  ->0  and  ( - r - )  tanh  -*00,  the  problem  can  best  be  handled 

O  A.  X 

O 

by  considering  the  case  of  a  IM  wave  propagating  along  a  thin  sheet  of 
dielectric  slab  in  space.-  Due  to  the  simple  geometry  of  this  equivalent 
problem,  it  can  be  easily  analyzed  (52).  The  results  will  not  be  given 
here . 

It  can  be  seen  from  the  above  numerical  results  that  the 

e  11 

mode  passes  smoothly  to  the  circular  mode  as  00 .  The 

=  3*0  curve  in  Figure  III-2  is  almost  identical  with  that  given  by 
Wegener  ( 17 ) . 

The  effect  of  the  variation  of  relative  dielectric  constant 
e^/e^  on  the  propagation  constant  can  be  seen  readily  from  Figtire  III-5' 
As  a  representative  example,  =  0.7  is  chosen  to  illustrate  the 
effect.  For  large  values  of  f  ^  approaches  to  the  characteris¬ 

tic  wavelength  of  the  rod  material,  =  X^/  ^  e^/ very  quickly;  6md 

for  quite  small  values  of  ^2,/^o  ^  ^l^^o  ^  *  ^^^o 

2q  cosh  I 

slowly  with  respect  to  - - -  .  One  may  therefore  deduce  that  for 

^o 

2q  cosh  Iq 

constant  - r -  and  constant  |  more  energy  is  carried  within 

° 


constant 


a  rod  with  higher  relative  dielectric  constant  and  that  the  field  out¬ 


side  the  rod  also  decays  faster  for  the  higher  dielectric  constant  rod. 
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It  may  then  seem  that  the  higher  dielectric  constant  rod  is  more 
desirable  as  a  transmission  waveguide.  Unfortunately  the  high  dielec¬ 
tric  constauat  material  usually  is  associated  with  a  large  loss  factor 
(63). 


(b)  The  Odd  Dominant  Mode,  the  Mode 

o  11 


.(1) 


Similar  procedures  as  those  used  for  the  HE,?'  mode  can  be 

e  11 

applied  here  to  analyze  numerically  the  characteristic  equation  for 

the  mode.  From  equation  2.5c-68,  setting  m  =  1  and  n  =  1, 

the  characteristic  equation  for  the  mode  is 

o  11 


1  .  1  ri  .  1 

Lx^  CeiU„)  /  Fek^(5„) 


^  y2  Qeii 


* - 5T - ^ <W> 

X  y 

where  (a  ,6  )  and  Q  (3  )  have  been  defined  by  equations  4 

X  X*S  XS  X  ITS  XS 

and  5  respectively.  The  discussion  given  in  section  3.2a  concerning 

the  validity  of  "3-term"  approximation  for  R..  (a  ,6  )  Q  (p  ,%  ) 

X  r*s  xs  X  P5  Xs 

also  applies  here.  F\irthermore ,  since  R,  (ot  ,6  )  Q  (P  ,X  )  a; 

X  1*5  Xs  X  X'S  XS 

R*(a*  ,6*  )  )  within  the  region  of  interest  (i.e.,  0  ^  x  ^  5^ 

X  1*5  Xs  X  iTs  xs 

0  -  y  -  3)/  the  qHE^^  mode  and  the  mode  are  also  considered  to 

be  degenerate. 

Equation  15  is  now  solved  according  to  the  method  outlined  in 
section  3'1»  Figure  III-6  shows  the  results  of  this  extensive  con5)uta- 
tion.  Again  y  is  plotted  against  x  for  various  values  of  ,  and 
a  constant  value  of  which  equals  2.5.  Combining  the  results 

shown  in  Figure  III-6  with  equations  12,  13  and  14,  the  guide  wavelength 
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Fig.  III-6.  Roots  of  the  characteristic  equation  for 
Minor  axLs/MaJor  axis  =  teuah  |  . 
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whlch  is  a  function  of  ,  uj  ,  the  size  of  the  guide;  and 

be  obtained.  The  normalized  guide  wavelength  x/x^  versus  the  nor- 

2q  cosh  Iq 

malized  major  axis  - - -  for  various  values  of  |  ,  and  a 

Xq  o 

constant  given  in  Figure  III-7.  The  pattern  of  the  curves 

is  very  similar  to  that  of  the  mode.  Again  as  expected,  no  cut¬ 


off  frequency  is  observed.  When  the  frequency  is  low,  i.e.. 


2q  cosh 


is  small,  the  guide  wavelength  becomes  that  of  the  free-space  wavelength 

and  most  of  the  energy  is  being  transported  outside  the  dielectric  rod; 
2q  cosh  Iq 

a  small  value  of  - - -  in^jlies  a  small  value  of  y  ,  which  means 

^o 

that  the  field  outside  the  rod  decays  at  a  rather  low  rate.  For  a  large 
2q  cosh  lo 

value  of  - r -  the  guide  wavelength  approaches  asymptotically  the 

o  ^ _ _ 

characteristic  wavelength  of  the  rod  material  =  X^/  \,  e^/ ,  almost 
all  the  energy  is  being  transported  inside,  and  the  field  outside  the 
rod  decays  very  rapidly. 

By  comparing  Figure  III-7  with  III-2,  it  is  noted  that  the  dif¬ 
ference  between  the  guide  wavelength  curves  for  the  mode  and  the 

mode  is  more  pronounced  as  |  gets  smaller.  The  normalized 
e  11  o 

guide  wavelength  of  the  mode  approaches  to  the  limit,  x/x^  = 


. - - —  2q  cosh  Iq 

1/  \/  €i  /  ,  faster.  For  example,  when  =  0.2  and - - -  =  0.9 

X  /x  for  the  HE^^^  mode  is  O.895,  while  X  /x  for  the  HEt^^  mode  is 

'  ^  /^ll  '  All 


=  0.9, 


o  11 


e  11 


0.987.  We  conclude  that  the  mode  binds  closer  to  the  dielectric 

o  11 

rod  than  the  mode.  When  5^  is  larger  than  3>  "the  guide  wave¬ 
length  for  the  mode  is  almost  .identical  with  that  for  the  HE,^^ 

o  11  i  e  11 

mode,  since  these  modes  are  degenerate  on  a  circular  dielectric  giaide. 

To  show  the  effect  of  the  variation  of  |  with  respect  to  X  A 

o  '0 

2q  cosh  lo 

for  a  fixed  valioe  of  - -  ,  Figure  III-8  is  introduced.  For  a 

2q  cosh  lo  °  (1) 

fixed  value  of  - r - —  the  curve  for  the  HE'/  mode  is  smoother 

Xq  0  11 

than  that  for  the  ^HE^^  mode.  It  is  again  quite  evident  that  when 
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Fig.  in -8.  Normalized  guide  wavelength  of  the  HEt-i  mode  as  a  function  of 
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2q  cosh  I 


is  very  small,  x/x  is  a  constant  vith  respect  to  the 

2q  cosh  I 

r  I  .  For  very  large  values  of  - - -  the  geometry- 


variation  of  I  .  For  very  large  values  of  - - -  the  geometry 

°  ^o 

of  the  rod  is  not  important  as  far  as  x/x^  is  concerned. 

In  Figure  III-9  the  normalized  guide  wavelength  is  plotted  against 


the  normalized  cross-sectional  area,  ( 


2q  cosh  6o  2 


)  tanh  6  for  various 


values  of  f  .  Unlike  the  case  for  the  HEi, '  mode,  it  seems'  that  the 

elliptical  rod  is  a  better  binding  geometry  for  the  mode  than  a 

circular  rod.  These  curves  for  various  values  of  are  quite  close 

to  each  other,  which  means  physically  that  the  field  lines  are  quite 

uniform  for  this  HE^^^  mode.  The  slight  differences  between  these 

curves  may  be  explained  by  the  fact  that  as  a  circular  rod  deforms  into 

an  elliptical  rod,  the  electric  lines  of  force  are  being  sqioeezed 

together  so  that  the  field  density  is  more  concentrated.  For  a  very  flat 

elliptical  rod,  the  electric  lines  of  force  are  almost  uniform  (the  field 

density  is  also  almost  uniform)  and  any  further  flattening  of  the  rod 

would  not  change  the  field  density  too  much.  Figure  III-IO  shows  the 

2q  cosh  Iq 

variation  of  the  X/X  versus  - r -  (with  |  =  const.)  curve  -with 

'  o  Xq  '  ^o  ' 

respect  to  the  various  values  of  •  The  behavior  of  these  curves 

for  the  mode  is  very  similar  -to  those  for  the  mode.  It 

o  11  e  11 

is  interesting  to  note  that  as  approaches  unity,  or  as 

2q  cosh  5 

approaches  infinity,  the  X/X  versus  - r -  curve  for  the 

o  Xq 

HEt^^  mode  becomes  identical  with  that  for  the  mode, 

o  11  e  11 

2q  cosh  Iq 

^  - r - -  00  and  |  -*  0  ,  this  problem  degenerates  to 


the  problem  of  TE  wave  propagation  along  a  thin  sheet  of  dielectric 


slab  in  space  (52). 


Fig.  III-9.  RormaUsed  guide  wavelength  of  the 

^HE^)  iBode  as  a  function  of  nornalised 

cross-sectional  curea. 
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3 « 3  Field  Configurations 

In  practice  the  field  configurations  are  most  quickly  found  by 
Inspection  of  the  mode  functions.  It  Is  found  that  the  patterns  of  the 
electric  and  magnetic  field  lines  are  quite  similar  to  those  knovn  In  a 
hollow  metallic  guide.  However,  owing  to  the  absence  of  the  metallic 
shield  around  the  dielectric,  the  field  Is  no  longer  confined  to  the 
Inner  space.  Furthermore,  d\ie  to  the  absence  of  conducting  walls  and 
therefore  the  absence  of  the  conduction  current,  all  the  electric  and 
magnetic  field  lines  must  form  closed  loops. 

Figures  Ill-lla  and  III-12a  show  the  transverse  cross-sectional 

field  distributions  of  the  mode  and  the  mode  respectively. 

e  11  o  11  X- 

The  longitxidinal  cross-sectional  views  of  the  field  distributions  of 
the  niode  and  the  mode  are  given  by  Figures  Ill-llb  and 

III-12b.  The  traveling  wave  patterns  are  shown  in  these  figures.  Solid 
lines  indicate  the  electric  lines  of  force;  dotted  lines  represent  the 
magnetic  lines  of  force.  Three  dimensional  sketches  of  the  field  con¬ 
figurations  for  these  two  dominant  modes  are  shown  in  Figures  III-I3  and 
III-14.  The  field  configurations  of  these  modes  are  quite  similar  to 
the  corresponding  dominant  modes  in  the  metal  tube  waveguide,  as  men¬ 
tioned  above;  a  sin5)le  method  of  excitation  is  thus  available.  The 
method  of  excitation  of  the  mode  and  the  HE^^^  mode  will  be  dis- 

cussed  in  greater  detail  in  Chapter  VI. 

3.4  Rate  of  Field  Decay 

The  dielectric  rod  waveguide  is  an  open  structure  and  hence  the 
field  is  not  confined  within  the  dielectric  rod.  Therefore,  the  guide 
is  susceptible  to  considerable  radiation  loss  when  it  is  mismatched  at 
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eljectrlc  lines  of  force  - 
magnetic  lines  of  force  -  - 


Fig.  III-ll. 


Field  configuration  of  the  mode, 

(a)  cross  section  * 

(h)  longitudinal  section 


Fig.  III-12.  Field  configuration  of  the  HE^  mode 
(a)  cross  section 
(h)  longitudinal  section 
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Pig.  A  sketch  of  the  electric  lines  of  force  for  the  ^HE^^mode 
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Input  and  output  ends,  when  it  is  curved,  or  when  extraneous  objects  are 

near  it.  The  knowledge  of  the  external  field  extent  and  the  rate  of 

field  decay  outside  the  rod  is  very  important.  It  is  known  from  the 

2q  cosh  I 

discussion  In  section  3*2  that  for  small  values  of  - -  most  of 

^o 

the  energy  is  transported  outside  the  dielectric  rod,  thereby  we  may 

expect  to  have  a  large  field  extent  and  a  slow  rate  of  field  decay.  For 
2q  cosh  lo 

large  values  of  - - -  most  of  the  energy  is  being  carried  inside 

0 

the  rod  so  that  the  field  extent  is  quite  moderate  and  the  rate  of  field 

decay  is  fast.  To  get  an  idea  of  the  variation  in  the  rate  of  field 

2q  cosh  Iq 

decay  and  the  field  extent  with  respect  to  the  cheuige  in  - r -  and 

^o 

nximerical  resiilts  for  the  longitudinal  electric  field  will  be  obtedned. 


(a)  The  Mode 

'  e  11 

According  to  section  2.5a,  setting  m  =  1  euid  n  =  1  ,  the 

equation  for  the  external  longitudinal  electric  field  of  the  mode 

e  il 


is 


00 


iPz 


r  =1 


(1) 


where 


pi 


Oek3(5^) 


13 

S3 

S3 

15 

S5 

S5 

17 

•  • 

St 

•  «  • 

s? 

•  •  • 

P33 

^53 

•  < 

^35 

•  1 

•  •  •  •  t 


I  •  •  •  •  • 


(2) 


-87- 


^33 

^13 

^73 

^35 

^15 

^75 

Gek^CSg) 

•  •  • 

•  •  • 

•  •  • 

^33 

^53 

•  •  • 

^35 

^5 

»  •  • 

•  ••  •••  ••• 


(3) 


The  above  infinite  determinants  may  be  solved  by  the  method  of  suc¬ 
cessive  approximations  (U7).  It  is  found  that 


»  P5  »  P7  •  •  • 


(M 


Therefore  the  external  may  be  approximated  by  only  considering  a 

few  terms  of  the  expansion.  The  normalized  external  longittidinal 
electric  field  is  given  by 


Gek^(l)  se*(Tj)  +  ^  Gek^(|)  se*(n)  +  ••• 
_ 

Gekj^(l^)  se*(Tj)  +^Gek^(6^)  se*(Ti)  +  ••• 


(5) 


where  E  is  the  axial  electric  field  intensity  at  |  and  E  is 

2  ''o 

the  intensity  of  the  axial  electric  field  at  I  =  •  Equation  5  is 

2q  cosh  lo 

coii5)uted  for  varioxis  values  of  and  - - -  with  ~ 

and  Tj  =  «/2  .  The  results  are  shown  in  Figures  111-15^  16,  17,  I8  and 
19  for  =  3.0,  1.0,  0.75,  0.5  and  0.3  respectively.  A  family  of 
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2q  cosh  Iq 

curves  for  various  values  of  - ^ -  are  shown  in  each  flgxire.* 

o 

These  figures  possess  similar  characteristics  as  far  as  the  variation 

2q  cosh  lo  , 

of  - r -  is  concerned.  The  axial  electric  field  decays  (not 


exponentially)  much  slower  and  extends  much  farther  for  smaller  values 

2q  cosh  Iq  2q  cosh  |q 

of - .  for  large  values  of  - r -  the  field,  decays 

Iq  O 

exponentially  quite  rapidly  and  its  extent  is  quite  small.  Physically 


it  means  more  energy  is  being  carried  outside  the  guide  for  smaller 
2q  cosh  Iq 

values  of  - - -  .  The  same  conclusion  was  reached  in  the  discus- 

o 

Sion  in  section  3*2 

To  observe  the  effect  of  axial  electric  field  extent  as  a  function 


of  frequency  for  various  values  of  eccentricities,  we  introduce  Figtire 

2 

III-20.  The  field  extent,  B/B^  ,  at  which  point  (E^/Ez^)  =0.1  is 

2q  cosh  Iq  2 

plotted  against  the  normalized  frequency  ( - - - )  tamh  1^  for 

^0 

various  |^  .  It  is  quite  evident  that  B/B^  is  larger  for  the  flatter 
elliptical  cross-section  rod.  As  frequency  becomes  very  high  the  nor¬ 
malized  field  extent  B/B^  approaches  to  unity  and  for  low  frequencies 
B/B^  can  get  very  leurge.  Since  B^  is  a  function  of  |^  ,  it  is 
somewhat  difficult  to  con5)are  the  absolute  axial  electric  field  extent 
of  a  circular  rod  and  that  of  an  elliptical  rod  having  the  same  cross¬ 


sectional  area  using  Figure  III-20.  Thus  Figure  III-21  is  introduced. 


B/x^  ,  the  normalized  absolute  field  extent,  is  plotted  against  the 

,2q  cosh  In. 2 

normalized  cross-sectional  area,  ( - )  tanh  |  for  various  |  . 

'  '  X  o  '*0 

0 


*In  each  figure  (E  /E^  is  plotted  against  B/B^  for  various  values 

of  2q  cosh  I  A  and  for  a  fixed  value  of  |  .  B  is  the  distance 

^  o'  o  o 

from  the  suds  to  the  point  of  observation;  B^  is  the  semi-minor  axis. 

(See  the  sketch  in  Figure  III -15) •  is  the  axial  electric  field 

strength  at  the  point  of  observation. 


bJitJJ 
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2<i  cosh  I  2 

( - tanh  -  NCSA 

o 

Pig.  III-20.  Pleld  extent  B/B„  where  (S  /B,  )^»  0.1  as  a  fynctlon 

'  o  s  lq  /i\ 

of  noraallsed  cross  sectional  area  for  the  Hl^'aode. 


fuactlon  of  noraBlizsd  cross-sectional  area  for  the 
Hl^^^eode  .  B  Is  the  distance  nsasured  free  the 

Q  Jba  2 

origin  to  the  point  of  observation  vhere  (I./E.  )  >0.1 

Z  Iq 


-9l^- 

Some  very  interesting  features  are  noted  in  this  figure.  For  the 
2q  cosh  Sq.o 

region  0.05  =  ( - i - )  tanh  5  ^  0.5  and  0.2  ^5  <00,  the 

°  ° 

veuriation  of  B/x^  is  quite  small;  it  varies  between  0.35  and  0.55* 

2q  cosh  Iq  2 

As  ( - - - )  teuih  I  approaches  infinity  so  does  B/X  ;  at 

A>  O  O 

O 

very  low  frequencies,  b/x^  approaches  zero.  As  the  cross  section  of 
the  rod  gets  flatter,  l.e.,  as  gets  smaller,  the  peaks  and  the 

valleys  of  the  curves  become  more  pronounced.  The  fact  that  within  a 
certain  frequency  range  the  absolute  axisU.  electric  field  extent  of  a 
flatter  elliptical  rod  is  actually  smaller  than  that  of  a  circular  rod 
of  the  same  cross-sectional  area  is  worth  mentioning.  As  a  typlcELl 
numerical  exan^ile,  we  choose  x^  =  3  cm  •  According  to  Figure  III-21 
the  absolute  field  extent  B  for  (E^/E^^)  =0.1  is  1.37  cm  for  a 
circular  rod  with  a  I.5  cm  diameter,  while  it  is  1.28  cm  for  ^  =  0.5 
elliptical  rod  with  the  same  cross-sectional  area. 

2 

Similar  curves  for  smaller  vedues  of  plotted. 

The  general  shapes  of  these  curves  remain  the  saxue,  only  the  peaks  and 
valleys  of  these  curves  are  more  pronounced. 

Figures  III-I5  through  III-I9  also  offer  a  convenient  way  of 
verifying  the  purity  of  the  mode  on  the  dielectric  guide.  (Experiments 
on  the  field  decay  properties  of  the  gHE^^  mode  have  been  carried  out 
and  the  results  are  reported  in  Chapter  VI.) 

Although  only  the  E  field  is  discussed  above,  it  may  be  shown 
that  the  other  electric  field  con^jonents  also  possess  similar  though  not 


identical  behaviors 
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(b)  The  Mode 

Setting  m  =  1  and  n  =  1  In  equation  2.5-63#  one  obtains  the 
ejqiresslon  representing  the  external  axial  electric  field  for  the 
mode; 

o  11 

?.  ,  (6) 
r  =t  1 
odd 


where 


!3 


“l  ■ 


“l3 

“53 

a?  3 

“15 

*  *  • 

“55 

•  •  • 

“75 

»  •  * 

“33 

“53 

*  •  • 

“35 

«  *  • 

“55 

•  »  • 

•  •  • 

•  »  • 

“33 

“13 

“73 

“35 

“15 

“75 

“37 

“l7 

a?? 

. . . 

•  •  • 

. . . 

ot^ 

ot^ 

33 

53 

“35 

“55 

(7) 


(8) 


The  above  Infinite  determinants  may  be  solved  by  the  method  of  succes- 
slve  approximation .  It  can  be  shown  that  h^  <  h^  «  h^  «  h^  •  •  •<  h^  . 
Therefore  the  value  of  the  external  axial  electric  field  can  be  approxi¬ 
mated  by  using  only  a  few  terms  of  the  e:q)anslon#  6.  The  normalized 
external  axial  electric  field  is 
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*  ^3  * 

E„  Fek,(|)  ce  (n)  +  Pek  (|)  ce  (tj)  + 

- i - - 3 - 3 -  _  (5) 

*o  Fek^(j^)ce*(T|)  +  ^  Fek^(5^)ce*(Ti)  +  ••• 

Numerical  computations  are  carried  out  assuming  t)  =  0  .  The  results 

are  shovn  In  Figures  III-22  through  III-25.for  6^  ranging  from 

2 

5^  =  3*0  lo  Iq  =  0*2  .  In  each  of  these  figxures  (E^/E^  )  Is  plotted 

2q  cosh  lo 

against  A/A  for  various  veilues  of  - -  .  2A  Is  the  major 

o 

axis  of  the  ellipse  while  A  is  the  distance  measured  from  the  origin 
to  the  point  of  observation  in  the  ij  =  0  plane  (see  the  sketch  in 
Flgiore  III-22) .  The  decay  characteristics  are  as  expected.  At  lower 
frequencies  the  field  decays  slower  since  a  larger  portion  of  the 
energy  is  carried  outside  the  rod,  and  at  higher  frequencies  the  field 
decays  faster  and  the  field  extent  Is  less,  since  more  energy  is  car¬ 
ried  inside  the  rod. 

Figure  III-26  which  is  similar  to  Figure  III-21  is  Introduced.  In 

this  figure  the  normalized  absolute  axial  electric  field  extent  2A/X^ 

2q  cosh  lo  2 

is  plotted  against  the  normalized  cross-sectional  area  ( - ^ - )  • 

o 

tanh  5^  ,  for  various  values  of  ;  the  point  of  observation  is  taken 

o 

to  be  the  point  where  (E^/E^^)  =  0.1  .  Again  as  o) -►  0  ,  2A/x.^ 

2q  cosh  lo  2  / 

approaches  to  zero;  and  as  ( - - - )  tanh  1^  ■*  oo ,  so  does  2A/X^  . 

However,  it  Is  Interesting  to  note  the  variation  of  these  curves  with 

respect  to  the  change  In  eccentricity.  Unlike  the  gHE^^  mode,  2A/x^ 

is  always  lEurger  for  flatter  elliptical  cross  section  rod.  This  is 

because  the  major  axis  of  a  flatter  ellipse  Is  always  longer  than  a 

rovinder  one  having  the  same  area. 

Flgiires  III-22  through  III-25  may  also  be  used  to  check  the  pvirity 
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Flg.  III-26.  Noraallzed  axial  electric  fieM  extent  as  a 

function  of  the  noroalized  cross-sectional  area  for 
the  ^HEj^^aode.  A  is  the  distance  asssured  froa  the 
origin  to  the  point  of  observation  vhere  (K  /S  )  •  0.1 

E  Zq 
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of  the  mode  on  the  elliptical,  dielectric  rod. 

3 . 3  Summary 

The  numerical  results  of  the  characteristic  eqiiatlons  for  the  two 

dominant  modes  are  obtained.  It  Is  found  that  for  the  mode  the 

e  J_L 

guide  wavelength  becomes  longer  as  the  elliptical  cross-section  becomes 
flatter,  and  for  the  mode  the  opposite  Is  tr\ie,  although  not  as 

pronounced.  As  expected,  there  Is  no  cutoff  freq\iency  for  these  two 
dominant  modes.  The  fact  that  these  two  modes  are  degenerate  when 
5^  -♦  00  is  also  demonstrated  numerically. 

Sketches  of  the  field  conflg\iratlons  for  these  modes  are  obtained. 
It  Is  observed  that  the  cross-sectional  views  of  these  modes  are  slmileur 
to  the  cross-sectional  field  pattern  of  the  dominant  mode  In  the  metal 
wavegxilde.  The  possibility  of  laxuachlng  these  dominant  dielectric  rod 
modes  by  mesms  of  the  metal  waveguide  Is  also  discussed. 

The  field  extent  of  these  modes  outside  the  dielectric  rod  Is 
considered.  It  Is  found  that  the  electric  field  extent  of  the 

o  11 

mode  Is  much  greater  than  that  of  the  HE^,  mode,  and  the  flatter  the 

0  A  li 

cross  section,  the  larger  the  contrast. 


IDl. 


CHAPTER  IV  -  AllfEHUATIOH  AHD  POWER  FLOW  CHARACTERIRTICS  OF  THE  DOMIWAHT 
MODES 


Having  o'btalnad  the  guide  vavelength  from  the  transcendental 
equation  and  Investigated  the  field  decay  characteristics,  It  wovild 
seem  appropriate  to  consider  the  attenuation  and  power  flow  properties 
of  these  dominant  modes.  Attenuation  Is  caused  by  inqperfectlon  of  the 
dielectric  material.  It  is  possible  to  include  the  dossy  cheLracterlstlcs 
of  the  dielectric  material  by  assuming  a  complex  dielectric  constant  to 
represent  the  permittivity  of  the  material.  The  effective  coiqplex 
dielectric  constant  is  given  by 


e 

e 


le 


+  ie 


le 


(1) 


Replacing  the  lossless  dielectric  constant  6^  (as  used  in  Chapters  II 
and  III)  by  the  effective  complex  dielectric  constant  g  and  substitut- 
Ing  into  the  equations  in  Chapters  II  and  III  where  applicable,  one 
notes  that  since  the  argiunents  of  the  Mathleu  and  modified  Mathleu  func¬ 
tions  are  complex,  the  roots  of  the  characteristic  equations  csin  no 
longer  be  real  and  must  be  complex.  Therefore  the  propagation  constant 
P  is  also  coiqplex  and  must  be  represented  by 


p  =  P'  +  la 


(2) 


where  P'  Is  the  new  real  propagation  constant  of  the  wave  on  the  lossy 
dielectric  guide  and  a  is  the  attenuation  factor  of  the  wave.  Of 
course  the  new  real  propagation  constant  P'  does  not  necessarily  equal 
the  propagation  constant  of  the  wave  along  the  lossless  dielectric  guide. 
The  numerical  solutions  of  these  coiqplex  characteristic  equations  are 
extremely  coo^licated  and  involved.  Even  the  anedytic  propexiiles  of  the 
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Mathieu  and  modified  Mathieu  functions  with  coniplex  eurgvunent  have  not 
been  well  investigated  and  understood  (45,^9)- 

However,  when  the  conductivity  of  the  imperfect  dielectric  is 
very  low,  in  other  words,  when  the  imaginary  part  of  the  effective 
dielectric  constant  given  by  equation  1  is  very  small,  i.e., 
it  can  be  assumed  (15)  that  to  the  first  order  approximation  the  dissi¬ 
pation  has  no  effect  on  the  field  configuration  of  the  wave,  which 
simply  remains  the  same  as  that  of  the  lossless  case.  Thus  the  propa¬ 
gation  constant  0  is  unchanged  by  the  presence  of  small  dielectric 
loss  and  therefore  the  analyses  carried  out  in  the  previous  chapters 
still  apply.  The  mode  functions  in  the  case  of  small  dielectric  loss 

differ  from  those  of  the  lossless  case  only  by  a  multiplicative  attenua- 
-Ctz 

tion  factor  e  ,  where  a  is  the  attenuation  constant  and  can  be 
calculated  by  a  perturbation  method  which  will  be  described  later. 

The  approximate  formula  for  the  attenuation  constant  a  will 
be  derived  by  the  Poynting's  vector  theorem.  The  problem  of  attenuation 
of  the  mode  along  a  slightly  lossy  dielectric  rod  will  then  be 

analyzed  analytically  and  numerically.  The  results  will  be-exhibited 
graphically.  Similar  considerations  concerning  the  problem  of  attenua¬ 
tion  of  the  mode  along  a  slightly  lossy  dielectric  rod  will  also 

be  made.  The. results  on  the  attenuation  properties  of  these  two  domi¬ 
nant  modes  will  be  discussed  and  compared  in  detail.  The  power  flow 


characteristics  of  these  modes  will  be  calculated. 
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k,l  The  Attenuation  Constant 

The  attenuation  constant  a  can  be  calculated  by  a  perturbation 
method,  provided  that  the  power  loss  per  wavelength  along  the  rod  is 
small  compared  to  the  power  flovlng  along  the  rod.  Since  there  is  no 
radiation  perpendicular  to  the  rod  at  large  distances,  the  power  flow 
is  only  in  the  z-dlrection,  i.e.,  only  along  the  axis  of  the  rod.  It 
has  been  pointed  out  earlier  that  the  fields  are  danped  exponentially  as 
they  propagate  along  the  rod,  and  if  their  attenuation  factor  is  a  , 
that  of  the  Poyntlng's  vector  is  2a  .  Therefore  the  attenuation 
constant  can  be  calculated  from  the  following  relation: 


2a 


1  ^ 
P  dz 


(1) 


where  P  is  the  time  average  transmitted  power  and  is  the  time 
average  power  loss  per  unit  length.  According  to  Poynting's  theorem 
(15)  we  have* 


S  +  =  -  J  •  E*  +  i(JD  (pH  •  H*  -  e  E  •  E*)  (2) 

oz  z  t  “t  —  “  “  “  ~  — 


where  S  is  the  longitudinal  component  of  the  Poyntlng's  vector  S  , 
z  " 

and  S.  is  the  transverse  part  of  S  .  Integrating  this  expression  over 

*  V 

a  cross-section  A  of  the  g\iide  (this  A  includes  the  cross-section  A^  of 
the  dielectric  guide  and  the  cross-section  A^  outside  the  dielectric  rod) 
we  get 


^  H  '  (ExH*)di  = 
c 


-  0 


d 


y  E  •  Ei 


E*  dA  +  4iaj 


[V  “e] 


(3) 


*  E*  or  H*  signifies  the  complex  conjugate  of  E  or  H  respectively. 
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whcre  P  is  the  time-average  pover  loss  per  vinit  length;  W  and 
oz  •  m 

W  are  the  time-average  magnetic  energy  and  electric  energy  per  unit 
e 

length  of  the  guide  respectively.  It  has  been  assumed  that  =  0 
outside  the  rod  and  J  =  o,E  inside  the  rod.  The  value  of  the  second 

—  d— 

integral  on  the  left  hand  side  of  equation  3  is  zero,  since  power  flows 
along  the  rod  only.  For  the  undisturbed  field,  “  '^e  '  thus  we  have 


^  ^  =  ■  ‘'d 


/  E  • 


E*  dA  . 


The  time-average  tremsmitted  power  is  given  by 


P  =  I  e^  •  (E^  X  i^)  dA  . 


E  and  H  are  the  transverse  components  of  the  electric  and  magnetic 

— t 

field  of  the  mode  under  consideration,  and  A  is  the  total  cross-sec¬ 
tional  area  of  the  guide.  Substituting  equations  4  and  5  into  1  we  get 
the  expression  for  the  attenuation  factor 


I  E  •  E-* 

A-' 


E*  dA 


2a  = 


where  the  unit  of  a  is  nepers/meter.  Changing  into  practical  units  we 


f  E'H 


E*  dA 


a  =  •  6 . 686  ■  ’  \  i“  - ^ 


(•db/meter) 


Hf)dA 


where  =  “>  ®l^d  ’  ^1  ^d  respectively  the  dielectric 
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constant  and  the  loss  tangent  of  the  rod.  It  may  he  noted  that  for  a 
plane  wave  propagating  in  an  infinite  homogeneous  medium  of  conducti¬ 
vity  0^  ,  the  ejjpression  within  the  absolute  value  signs  becomes 
1/  \j  €^/ where  is  the  dielectric  constant  of  the  surrounding 
medium. 


Let  us  now  consider  the  integrals  within  the  absolute  valxie  signs. 


In  elliptic  cylinder  coordinates  these  integrals  can  be  e:q?ressed  as 


follows : 


(E 


E*  +  E  E*  )q^(sinh^|  +  sin^T  )dT|d| 

Hi  Z;j_ 


(8) 


and 


(E.  X  H*  )dA 
~‘'o  ~^o 


(E-  H*  -  E_  H*  )q^(sinh^|  +  sinS)dHd| 

^1  ‘1  '1  ^1 


00  2k 

+  I  I  (e  h*  -  E  H*  )q^(sinh^5  +  sin^T))dild|  ,  (9) 

Q  *0  o  o 

^o 

where  the  subscript  1  and  the  subscript  0  represent  the  inside  and  outside 
regions  of  the  dielectric  rod  respectively,  and  q  is  the  semifocal 
length  of  the  ellipse.  is  the  boundary  surface  of  the  elliptical  rod. 


k.2  The  Attenuation  Factor  and  the  Power  Distribution  Characteristics  of 

the  HE^^^  Mode 
_ e  11 

The  field  components  of  the  HEt^^  mode  can  be  obtained  from  sec- 

e  11 

tion  2.5a  by  setting  m  =  1  and  n  =  1.  They  are,  for  region  1, 

(0^1^  lo), 
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\  (k^-p^)p 


(k^-P  )p 
and  for  region  0  ( §  ^  |  <  oo ) 
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<!— I ,  r  r  r 
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(7) 
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00 
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odd 

00 


E  =  y  T-  ^L/ek^(Oce*(Ti)+  P^Oek  (|)se*(Ti)l  (12) 

’fc  (k^-p)pr=l''  P 


odd 


2  2  2  2  2  2  1/2 

where  k^  =  co  k^=  o)  and  p  =  q(slnh  |  +  sin  T|)  . 

and  Pj.  are  the  arbitrary  constants  which  are  related  by  the  boundary 


-107- 


conditlons  and  are  given  by  (see  section  2.5a,  equations  43  through  45) 
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(14) 


-  ,  (15) 

A.  „ 
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The  symbols  used  In  the  above  expressions  have  all  been  defined  in 
Chapter  II . 

Upon  examining  the  second  integral  in  eqiiation  4.1-9,  it  can  be 
seen  that  it  would  be  extremely  laborious  to  carry  out  this  integral  if 
many  terms  of  the  expansions  representing  the  field  con5)onents  outside 
the  guide  are  required.  Fortunately  it  is  found  that,  within  the  region 
of  otir  Interest,  l.e.,  for  ^  0.2  ,  x  ^  5  ,  and  O-l  -  7  -  3  > 


P 

»  ,  so  that  the  expressions  representing  the  field  components 

outside  the  guide  can  be  approximated  by  only  two  terms  of  the  infinite 
series  e^qpanslon.  In  other  words,  the  infinite  series,  representing  the 
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field  components  in  region  0,  converge  rather  rapidly,  providing  that 

is  not  too  small.  For  Instance,  it  is  found  numerically  that  when 

5  1.0,  the  value  of  H  (see  equation  7)  computed  from  the  first  two 

o 

terms  of  the  infinite  series,  is  accurate  to  the  fifth  significant  figure 
at  the  worst;  when  1.0  -  -  0*5,  it  is  accurate  to  the  fourth  signi¬ 
ficant  figure  at  the  worst;  and  when  0.5  -  -  0.2,  it  is  acciirate  to 

the  second  significant  figure  at  the  worst.  It  can  be  seen  that  the 
accuracy  gets  better  as  gets  larger,  assuming  that  the  same  number 

of  terms  is  used. 

We  are  now  in  a  position  to  consider  the  Integrals  in  equations 
4,1-8  and  4.1-9.  Substituting  equations  2,  5  an<i  6  into  4.1-8,  one 
obtains 


Separating  the  angular  and  radial  integrals  and  simplifying,  one  gets 
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E*)dA 


B,  2  6 


v,2.  2.2  ^ 

cosh  I  q  Ai  f- 


o 
Bn 


2  e 

*  ‘IT>  > 


.  j  1  ^7'! 

eosh^S 


(21) 


2  2  fl  2  2 

X  +  y  e  X  +  y 

where  =  — j— — 2  ^  q  =  — ^r~2  ’  ^  "  q^cosh  L(k^-  ^)  ,  and 

X  (^  -  1)  +-^  y2  0  1 

0  o 

2  2  2  2 

y  =  -q  cosh  P  )•  through  are  the  angular  integrals  which 

can  be  integrated  analytically  and  are  given  in  Appendix  B.l;  Iq  through 
are  the  radial  integrals  which  must  be  evaluated  numerically  and  are 
defined  in  Appendix  B.2. 

Substituting  equations  3  through  6  Into  the  first  Integral  In  equa¬ 
tion  U.1-9  gives 


2n 


0 


(E.  H*  -  E  H*  )q^(slnh^|  +  sln^q)dT)d| 

5i  \  \  §2. 


p^a: 


2  /K.  2it 

li^o  f  f 

j  j 

''^1  ^  ^  0  0 


U) 


\A^ 


(ce^(l)cej^^(q)  +  CeJU)ce^(q)) 


€  (a. I €  B  2 

+  "T  “p“°  ®0l(  ’i)  + 

2 

+  (1  +  — ^)(Ce^(|)Se[(|)cej^(q)se^(q) 


-  Ce£(|)Se^(|)ce^(T))se|(ii)) 


dT|  d|  .  (22) 


Simplifying  this  e:q)ression  we  have 


where  C.  ,  C„,  x  and  y  have  been  defined  earlier.  The  angular  inte- 
grals  through  and  the  radial  integrals  Ig  through  are 

all  given  in  ;^pendix  B.l  and  Appendix  B.2  respectively.  The  above  ex¬ 
pression  23  also  represents  the  portion  of  the  total  transmitted  power 
being  carried  inside  the  dielectric  rod.  The  other  part  of  the  trans¬ 
mitted  power  which  is  carried  outside  the  rod  is  represented  by  the 
second  integral  in  equation  4.1-9.  Substituting  the  "two-term-approxi- 
mation"  of  equations  9  thro\agh  12  into  the  second  integral  in  equation 
4.1-9^  one  gets 


2jt 


^o 


(E,  H*  -  E  H*  )q^(sinh^|  +  sin^Ti)dTj  d| 
^o  \)  \)  'o 


p')"  i  i  ^ 

o 


.':v 


P^L^(l  +  -^)(B3)|  dri  di  , 


(24) 


in  whlc^ 
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(Bl)=  Fek^(|)ce*'^(Ti)+  (^)  Fek^(|)ce*^' (t))+  2(^)Fek^(  |)Fek^(  |)ce*'(Ti)  ce*'(Ti) 
+  Fek^ii)  ce^(Ti)+  (^)^Fek'^^(5)ce*^(,l)+  2(^)Fek^(  OFek^(  |) ce*( n) ce*( tj)  , 
P  2  ^ 

(B2)=  Gek^(|)se*'^(Ti)+  (^)  Gek^(  |)se*'^(n)  +  2(^)Gek^(  |)Gek^(  5)  se*'^(T|)se*'(Tj) 

P 

+  Gek^^(|)se*^(Ti)+  (^)^Gek'^(|)se^Ti)+  2(^)Gek|(  |)Gek^(  |)se*( r,)  se*(Tj) , 

and 

(B3)=  Fek^(4)Gek|(|)ce*’(Ti)se*(i^)  -  Fek^(  |)Gek^(  |)  oe*(Tj)  se*(  tj) 

r  T 

+  [Fek3(|)Gek^(Oce*'(Tj)3e*(Tj)  -  Fek^^(  e)Gek^(  0  se* '  (t^)  ce*(  t^) 

+  ^I^Fek^(^)Gek^(^)ce#(T))3e*(T|)  -  Fek^(  5)Gek^(  |)  ce*(T])  se* '  (ti) 

+  ij^)  (^)[Fek^(OGek^(|)ce*'(Ti)se»(T))-  Fek^(|)Gek^(  |)se*(tj)ce*(Ti)  . 

Simplifying  equation  24  results  in 


cosh^l  ^  ^ 

o  1  e 

o 


T  L.  P 


(A(Tr)(i  +  CJ  (Bill) 


where 


(BI)  =  dil-,.  1-1^2)"  (lJ)  ^^3^23"  '24^4)"  ^26^4^  ' 


Po2 


(BII)  -  (127^)"^  ^28^3^'*'  ^^29^9'*'  ^30^10^^  ^^31^11  ^32^12^  ' 


(25) 
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(Bill)  =  (I'ji'j-  (^)(i^,ii7- 


emd  C^,  x  and  y  have  been  defined  earlier.  Again  the  angular 

integrals  throxjgh  which  can  be  integrated  analytically,  and 

the  radial  Integrals  through  which  must  be  evaluated  numerl 

cally,  are  tabulated  in  Appendix  B.l  and  Appendix  B.2.  It  should  be 
noted  that  the  numerical  values  of  the  higher  order  terms,  neglected  in 
the  above  expressions  (Bl),  (BIl),  and(BIIl)  are  of  the  order  of 
P^/P^  or  •  Expression  25  represents  the  portion  of  the  power 

being  transmitted  outside  the  rod. 


Substituting  equations  21,  23,  and  25  into  the  expression  within 
the  absolute  value  signs  of  equation  4.1-7,  we  get 


and 
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R  Is  related  to  the  attenxiation  constant  a  in  db/meter  by  the  follow¬ 
ing  relation 

A  hAr  r~. 

(27) 


8 . 686  fa 

“  “  —  “dW^" 

V  o 


where  is  the  conductivity  of  the  dielectric  rod. 

Using  the  results  given  in  Chapter  III  regarding  the  relationship 

between  x  and  y  for  various  values  of  |  and  p,/e  ,  numerical 

o  1  o 

compatauion  of  R  can  now  be  carried  out.  Ail  radial  integrals, 
thrc-ogh  and  1^^  through  are  evaluated  numerically  using  Simp¬ 

son's  rij.e  (so).  Results  of  tliis  very  lengthy  computation  are  siaown  in 
Figure  IV-1*.  In  tliis  figure  the  value  R  ,  which  is  directly  propor¬ 
tional  to  the  attenuation  constant  (see  equation  27)  is  plotted  against 


the  normalized  major  axis  (NMA), 


2q  cosh  I 


,  for  various  values  of  | 


ranging  from  =  3-0  to  0.2.  It  is  assumed  that  =  2.5  . 

For  sufficiently  large  values  of  NMA,  R  tends  toward  the  plane-wave 

value  V\/^/^Q  values  of  ;  for  small  enough  values  of  NMA 

R  can  be  chosen  as  small  as  desired.  This  behavior  is  attributed  to 

the  fact  that,  when  NMA  is  sufficiently  large,  almost  all  of  the  energy 

of  the  wave  is  transmitted  inside  the  rod  ;  and  for  small  values  of 

NMA  almost  all  of  the  energy  is  outside  the  rod.  (it  has  been  assumed 

that  the  dielectric  suirrounding  the  rod  is  perfect.)  It  is  also  clear 

*It  takes  almost  2u  minutes  of  continuovis  computation  by  the  IBM7090 
computer  to  obtain  each  curve. 

**It  is  noted  that  when  NMA  is  very  large  the  attenxiation  factor  R  is 
numerically  identical  with  the  attenuation  factor  of  a  certain  wave¬ 
guide  mode  propagating  in  a  perfectly  conducting  metal  tube  waveguide 
filled  with  the  same  dielectric  material  as  that  of  the  dielectric  rod 
under  consideration. 
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^  fl^ 

Pig.  IV-1.  Attenuation  factor  R  for  the  mode  as  a  function  of 

normalized  major  axis. 
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that  R  tends  to  the  limit  1/  \J  much  slower  as  gets  smaller 

tind  that  flatter  elliptical  dielectric  guide  possesses  lower  loss  char¬ 
acteristics.  This  may  he  explained  by  the  fact  that,  according  to 
Figure  IV-1,  for  a  constant  value  of  NMA,  smaller  rod  has  less 

volume  of  dielectric  material  and  therefore  lower  dielectric  loss.  How¬ 
ever,  this  is  not  the  only  reason.  If  we  plot  R  against  the  normalized 

2q  cosh  I  2 

cross-sectional  area  (NCSA),  ( - r - —)  tanh  |  ,  for  various  values  of 

with  -  2.5  ,  as  in  Figure  IV-2,  the  same  effect  (i.e.,  lower 

loss  for  smaller  5^)  of  a  lesser  degree  can  still  be  observed.  As  the 

elliptical  cross-section  gets  flatter,  the  field  of  the  wave 

e  11 

spreads  out  more  so  that  the  total  integrated  effect  on  the  attenuation 

indicates  that  this  type  of  field  distribution  offers  less  loss.  The 

shape  of  these  curves  in  Figure  IV-2  shows  that  the  attenuation  factor 

can  be  made  extremely  small  if  a  very  flat  strip  is  used.  For  example, 

when  NCSA  =  0.15  the  attenuation  constant  a  of  the  mode  can  be 

e  11 

made  ten  times  smaller  if  =  0.3  elliptical  rod  is  used  rather  than 

a  circular  rod,  and  a  may  be  almost  90  times  smaller  if  e  =0.2 

o 

elliptical  rod  is  used.  It  is  interesting  to  compare  the  axial  electric 
field  extent  of  these  rods  corresponding  to  the  above  example.  According 
to  Figure  III-21,  when  NCSA  =  0.I5  the  axial  electric  field  extent 
where  (E  /E  )  =  0.1  for  ^  =  0.3  elliptical  rod,  is  0.47,  and  for 
=  3*0  elliptical  rod  it  is  0.435. 

The  fact  that  the  variation  of  slopes  with  respect  to  NCSA  in 
Figure  IV-2  is  smaller  for  flatter  rods  in  the  low  loss  region,  is  qviite 
significant.  It  means  that  a  small  imperfection  in  the  dimensions  of  a 
flatter  rod  would  induce  a  smaller  change  in  the  attenuation  factor  R  . 
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lt  is  interesting  to  note  the  distribution  of  the  transmitted 
power.  There  is  a  very  close  correlation  between  the  percentage  of 
power  carried  inside  the  rod  and  the  loss  factor  of  the  wave.  With 
the  help  of  Poynting's  vector  theorem,  one  can  easily  calculate  the 
percentage  of  power  transmitted  inside  the  dielectric  rod.  It  is 


V  S 


(28) 


where  f^  and  f^  are  given  in  equation  26.  Numerical  results  of 
equation  28  are  given  in  Figure  IV-3  in  which  P. /P^  is  plotted 
against  NMA  for  various  values  of  ;  e^/ is  assumed  to  be  constant 
and  is  equal  to  2.5  •  It  is  observed  from  Figures  IV- 2  and  IV-3  that 
a  higher  percentage  of  power  transmitted  inside  the  rod  corresponds  to 
a  higher  attenuation  factor  and  more  power  is  carried  inside  the  circu¬ 
lar  rod  than  an  elliptical  rod  of  identical  cross-sectional  area. 

Figure  IV-3  also  confirms  the  fact  that  more  power  is  carried  inside 
the  rod  as  the  frequency  gets  higher. 

The  I  =  3*0  curve  in  Figure  IV-l  corresponds  very  well  with  the 
o 

published  results  for  the  circular  dielectric  rod  (l8).  The  analytic 
expression  of  the  loss  factor  R  for  the  degenerate  circular  dielec¬ 
tric  rod  can  easily  be  derived  from  equation  26.  Noting  that  as 
-♦  oo,  q  -►  0  ,  q  cosh  -►  a  and  qe^/2  -►  p  ,  where  a  is  the 

radius  of  the  clrcule  and  p  is  the  radial  component  in  the  polar 
coordinates,  one  gets. 
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Q(P  ,X  )  -*  1  . 


Equations  21,  23  and  25  degenerate  to  the  known  expressions  associated 
with  the  loss  factor  of  the  dominant  mode  propagating  along  a  circular 
dielectric  rod  (21);  they  are,  respectively. 
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in  which  p  =  x  2-  and  n  =  y  -  .  The  attenuation  constant  of  the  HE, , 
a  a  11 

wave  on  a  circular  dielectric  rod  can  easily  he  obtained  by  substituting 
the  above  expressions  into  equation  4.1-7. 
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4.3  The  Attenuation  Factor  and  Power  Dtstrlbutlon  Characteristics 

of  the  Mode 

o  11 


The  case  of  the  HE'  wave  can  be  analyzed  in  a  similar  manner 
o  11 

to  that  of  the  wave.  From  section  2.5c  one  obtains  the  field 

e  11 

components  of  the  qHe|^^  mode.  They  are,  for  region  1  (0  ^  |  ^ 
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where  k^=  co  and  p  =  q(sinh  |  +  sin  n)  . 

and  are  arbitrary  constants  which  are  related  by  the  boundary 

conditions  and  are  given  by  (see  section  2.5c  equations  70^  71  and  72,) 
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(19) 


The  attenuation  factor  of  the  mode  can  he  calculated  in  a 

o  11 

similar  maimer  as  that  of  the  mode.  It  can  be  shov/n  that 

e  11 

gl/si  >  g^/a^  »  g^/a^  »  •  • .  »  /a^  ,  and  h^/a^  >  h^a^  »  h^/a^  » 

•••  »  4^/aQ_  )  provided  that  ^  0.2,  0  ^  x  ^  5  ^  and  y  ^  3*0  . 

The  first  two  terms  of  the  Infinite  series  are  sufficient  to  represent 


-125- 


the  field  components  in  region  0  . 

Substituting  the  appropriate  field  expressions  1  through  12  into 
4.1-7  one  finds  after  some  manipulation  that 

8.686  , 


a 


2  d 


^  R'  ,  (db/m) 


(20) 


in  vhich  R'  =  F^/  F^+F^  , 


where 


/(El*  E* 


E*)dA 
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u2,  2  2  u 
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b  2  e 
In  o 


C2(liili^+  IiqI^)  +  (a^)  ^  (I^lg+Iglj^) 


^  2V^2  Jr  ^We-  h2^5^ 


b  2  e 

+  (^) 


.  .  .  1  "8(^2“  ^7c) 

^l4_/2'^  2  .2. 


NW 


cosh  £ 


h;,)*  e  dA 
— tl  tl  — z 


F  = 
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=  C, 


^2  2  2  fr 

cosh  l^q  I  — 


(21) 


O  1 
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+  (-i) 

V\i 


6  e 
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£^\/^2^^^13^6"  ^12^5^ 


,  (22) 


and 


/  (E.  X  ^  ) •  e  dA 
J  -to  --to  -z 
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cosh^l  q^a^ 

°  ^  J^o 


^1 
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(BI)',  (BII)',  and  (Bill)'  are  respectively 


31 


^^i^2l''’  ^2^22^  ^^3^23"^  ^^25^5  ^26^61^  ' 


(^3^^1.4-^3343^  ^  (il'  (4646-  4?45^^  (^)"(4848-  4747^ 


where  through  1^  ,  Ig  throiigh  through  ,  1^^  and 

^14nN  given  in  Appendix  B  .  The  loss  factor  R'  for  the 

mode  as  a  function  of  NMA  is  con5)uted  for  various  values  of 

ranging  from  |^  =  0.2  to  =  3*0*  The  relative  dielectric  constant 

€^/€^  is  assumed  to  he  constant  and  equals  2.5.  The  results  are  plotted 

in  Figure  IV-4.  The  attenuation  factor  R'  for  the  mode  varies 

with  frequency  in  a  similar  way  as  that  for  the  mode,  viz.,  the 

attenuation  factor  R'  approaches  1/  /  e^/ as  frequency  approaches 

Infinity  and  R'  can  he  made  arbitrarily  smaJJ.  hy  lowering  the  frequency. 

It  shoxild  he  noted  that  the  slope  of  the  curve  'for  the  elliptical  rod  in 

the  low  loss  region  is  quite  steep,  i.e.,  a  small  variation  in  NMA  woiuld 

cause  a  rather  significant  fluctiiation  in  R' .  It  is  quite  obvious  that 

the  he!  4  node  is  more  suitable  than  the  He!4  mode  as  a  transmission 
e  11  o  11 

mode . 

The  distribution  of  the  treuasmitted  power  as  a  function  of  frequency 
can  easily  be  coaqputed.  The  percentage  of  power  carried  inside  the  dieleC' 


trie  rod  is 


0.16  0.4  0.6  0.8  1.0  1.2  1.4 

2q  cosh  i  /\ 

O'  O  ’ 

(1) 

Pig.  IV-4.  The  attenuation  factor  R'  for  the  '  mode  as  a  function  of 

the  normalized  major  axis. 


dielectric  rod  s\q>portliig 
;ed  major  axis. 
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^2 


(24) 


where  Fg  and  are  given  hy  equations  22  and  23  respectively. 

Figure  IV-5  shows  the  variation  of  Pj/P!  as  a  function  of  NMA  for 

•L  X) 

various  values  of  ;  ^l^^o  2*5  •  The  behavior  of  these  curves 

is  as  expected.  More  power  is  carried  inside  the  rod  as  the  frequency 
gets  higher.  Again  there  is  a  very  close  correlation  between  the 
amount  of  power  carried  inside  the  rod  and  the  value  of  the  attenuation 
factor . 


4.4  Summary 

A  detailed  analysis  on  the  attenuation  characteristics  of  the 
mode  and  the  mode  propagating  along  an  elliptical  dielec¬ 

tric  rod  are  carried  out  in  this  chapter.  Numerical  results  are  obtained. 
It  is  found  that  a  thin  elliptical  dielectric  rod  operating  in  the 

dominant  mode  is  a  better  guiding  structure  than  a  circular 

e  11 

dielectric  rod  operating  in  the  dominant  HE^^  mode,  because  the 

mode  has  much  lower  loss  on  a  flat  elliptical  rod  than  on  a  circular  rod 

of  identical  cross-sectional  area. 

It  would  be  interesting  to  compare  the  attenuation  constant  of  the 
mode  with  the  attenuation  constants  of  some  well  known  metallic 

e  11 

waveguide  modes  in  the  millimeter  wavelength  region.  The  values  of  the 
attenuation  constants  for  various  kinds  of  waves  are  tabulated  in  the 


following  table. 


Atteniaatlon 
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f 


The  Polystyrene  has  a  dielectric  constant  of  2,5  and  a  loss  tangent  0,  of  0.001. 
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CHAPTER  V  -  ELLIPTICAL  DIELECTRIC  ROD  RESONATOR 

To  conclude  the  theoretical  analysis  of  surface  wave  propagation 
along  an  elliptical  dielectric  rod,  we  Include  here  the  analysis  of  the 
Q  factor  of  the  elliptical  dielectric  rod  cavity.  The  earliest  work 
on  dielectric  resonators  was  carried  out  by  Richtmyer  (53)  in  1939-  He 
developed  the  theory  of  operation  for  several  interesting  dielectric 
resonators  of  simple  shapes,  such  as  the  spherical  dielectric  cavity 
and  the  "doughnut"  shape  dielectric  cavity.  The  dielectric  tube  resona¬ 
tor  was  first  used  by  the  group  in  the  Northwestern  University  (5^) • 
Later  in  1)^39  Becker  and  Coleman  (26)  made  use  of  the  dielectric  tube 
resonator  to  generate  millimeter  and  submillimeter  waves  and  to  operate 
as  a  frequency  meter.  Most  recently  Snitzer  (4o)  proposed  the  \ise  of 
dielectric  rod  cavity  as  a  mode  selector  in  laser  operation. 

In  the  present  problem  the  dielectric  rod  cavity  consists  of  an 
elliptical  dielectric  rod  suitably  terminated  at  its  ends  by  suffici¬ 
ently  large  flat  metal  plates  which  are  perpendicular  to  the  axis  of  the 
rod  (see  Figure  V-l) .  At  resonance,  the  length  of  the  cavity  L  must 


Figure  V-l.  The  Elliptical  Dielectric  Rod  Resonator 
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be  n  (n,  an  integer),  where  X  Is  the  guide  wavelength  of  the 

particular  ^de  lander  consideration  and  is  a  function  of  X^  ,  e^, 

and  the  size  and  shape  of  the  dielectric  rod.  The  relations  between 

X  and  the  mentioned  physical  c  nstants  are  determined  by  the  bo’jnd- 
g 

ary  conditions.  Only  the  and  the  modes  of  the  dielec- 

•'  e  lln  o  lln 

trie  rod  resonator  will  be  considered  in  this  chapter. 

The  Q  factor  of  a  resonator  is  indicative  of  the  energy 
storage  capability  of  a  structure  relative  to  the  associated  energy 
dissipation  arising  from  various  loss  mechanisms,  such  as  those  due  to 
the  imperfection  of  the  dielectric  material  and  the  finite  conductivity 
of  the  end  plates.  The  common  definition  for  Q  is  applicable  to  the 
dielectric  rod  resonator,  sind  is  given  by  (15) 


total  time-average  _ 

0  =  0)  energy  stored  W 

o  - ^ f 

average  power  loss  o  p 


(1) 


where  u)^  is  the  frequency  of  oscillation.  The  above  approximate  ex¬ 
pression  is  valid  when  Q  »  1  . 

In  our  case  the  time-average  power  dissipation  P  consists  of 
two  parts,  the  power  loss  due  to  the  dielectric  rod  and  that  due  to 
the  metal  end  walls 


^  ^dielectric  ^  ^wall 

The  power  dissipation  due  to  the  dielectric  rod  is  given  by  (15) 


(2) 


p  =  -o  I  I  (E-.  ’  E*  )dAdz 

aielectrlc  dj  J  —1  —  1' 

0 

while  the  loss  due  to  the  end  wall  is  (15) 


(3) 
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(4) 


All  the  symbols  in  the  above  expressions  have  the  same  meaning  as  those 
defined  in  the  previous  chapters.  The  arbitrary  constants 
and  which  are  related  by  the  boundary  conditions  are  given  by 

equations  4.2-13  through  4.2-19.  Expressions  1  through  4  satisfy  the 
boundary  conditions  on  the  surface  of  the  dielectric  rod  and  at  the  end 
z  =  0  .  To  make  them  also  satisfy  the  boundary  conditions  at  the  other 
end,  z  =  L,  we  restrict  ft  in  such  a  way  that  BL  =  nn  where  n  is 
an  integer,  (i.e.,  L  =  nX.  /2). 

Substituting  the  proper  field  expressions  into  5-5  >  carrying  out 
the  integrations  where  possible  and  retaining  enough  terms  of  the  expan¬ 
sion  to  give  the  same  order  of  approximation  as  obtained  in  Chapter  IV, 
one  finally  arrives  at  (after  some  rather  lengthy  algebraic  manipula¬ 
tions)  the  expression  for  the  energy  stored  in  the  cavity  for  the 
mode, 


W  = 


I  €  r 

I  (E, •  E*)dV  +  ^  (E  •  E*)dV 

J  — i  — 1  2  J  -O  -O 
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=  cosh  q  -jp  —  C, 


(5) 


€  2  € 
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2  2 

BI,  BII,  BIII|  X  and  y  are  given  in  Chapter  IV.  The  inte¬ 
grals,  i.e.,  I's  or  I'  's  ,  are  defined  in  Appendix  B. 


The  power  dissipation  due  to  the  dielectric  loss  is 

^o  2jt 


‘dielectric 


=  -  ^  J  /  (h'  ^ 

0  0 


E*)p  dT|d| 


.  '^d^  2  2  2  „ 

=  —  q  cosh 

o 


(6) 


where 


B  2  € - 

C  =  (1)  (-^) 
d  A^ 


I  . 

^  2  cosh^^  - 

o 


+  C, 


B  2  ^  B  r 

^  ^11^1+^  ^^^2  ^A^^  ^^12^5" 


Again  the  integrals  are  defined  in  Appendix  B . 

Another  source  of  power  loss  in  this  cavity  is  caused  by  the 
finite  surface  conductivity  of  the  reflecting  end  plates.  This  loss  may 
be  computed  from  equation 


P  =2 
wall 


2  ^s  J  («t- 


dA 


at  z=0 


=  R  A^cosh^l  q^C  (7) 

5  X  W 


where 

C..  =  C, 


B,  2  e_ 


(l8li+  ^2^A^^  '‘^10^3’^  ^11^4^ 


^  ^A^^  V  ^^12^5’  ^13^6^ 

4  f  L-  2  Pi  o  e  P.  /T"  ^ 

*  -  '(7^)(BI}+  Cj7i)^:;^(Bn)+  2(3)(^) 


Rearranging  expression  5-1^  ‘we  get 
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i  =  ^  _  ^dielectric  ^vall  _  1^1 

^  0)  W  0)  T7  CO  ^ 


Is  the  Q  factor  of  the  cavity  If  the  end  plates  are  perfectly  con¬ 
ducting,  and  is  the  Q  factor  of  the  cavity  If  the  dielectric  is 

perfect .  According  to  equation  6  we  have 


and 


CO  W 


dielectric 


o 


Q  -  W  L 

^  P  26  C 

wall  ^ 


(9) 


(10) 


where  0^  is  the  loss  tangent  of  the  dielectric  rod  and  6  is  the 
skin  depth  of  the  end  plates. 

The  expressions  and  C^/C^  are  evaluated  numerically  and 

the  results  are  siiown  in  Figinre  V-2  in  which  and  are 

plotted  against  the  normalized  cross-sectional  area  (NCSA)  for  various 
values  of  with  -  2.5*  For  small  values  of  NCSA,  C^C^  can 

he  very  large,  thus  can  also  he  very  large.  This  is  because  most 

of  the  energy  is  outside  the  rod.  As  NCSA  approaches  infinity,  C^C^ 
approaches  1/20^*  Again  one  notes  that 

the  flatter  the  elliptical  cross  section,  the  higher  the  factor. 

It  is  worthwhile  to  take  notice  of  the  behavior  of  as  a 

fiinction  of  NCSA.  For  an  ordinary  cylindrical  metallic  waveguide  of 
simple  cross-sectional  shape,  terminated  at  both  ends  by  short- circxxi ting 

*Incidentally,  the  Q  of  a  section  of  perfectly  conducting  metallic  wave¬ 
guide,  terminated  at  both  ends  by  perfectly  conducting  end  walls  and 
filled  with  a  dielectric  material  with  a  loss  factor  of  0,  is  also 
1/2^^. 
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plates,  the  Q  factor  resulting  from  the  Imperfection  of  the  end  plates 
Is  l/2&  .  L  is  the  length  of  the  guide  and  5  is  the  sld.n  depth  of 
the  end  plates.  This  Q  factor  is  independent  of  the  type  and  order  of 
the  mode  under  consideration  as  long  as  the  mode  is  either  of  TE,  or 

TEM  type  and  not  of  a  hybrid  type.  It  means  that  for  this  type  of 
cavity,  is  always  imity.  However,  is  no  longer  a  cons¬ 

tant  (see  Figure  V-2)  if  a  hybrid  wave  is  present.  This  characteristic 
is  probably  due  to  the  fact  that  the  TE  and  TM  waves  are  Inextricably 
coupled  to  each  other  on  a  dielectric  rod  except  for  the  circularly  sym¬ 
metric  waves . 

It  is  also  noted  that  is  independent  of  the  length  of  the 

cavity  and  is  directly  proportional  to  the  length  of  the  cavity. 

The  total  Q  of  the  cavity  can  be  computed  from  the  knowledge  of 

eind  using  equation  8.  For  a  very  long  cavity,  Qy  ^  therefore 

^total'^S  * 


5.2  The  Q  of  a  Cavity  Supporting  the  Mode 

For  the  sake  of  completeness,  we  include  here  the  analysis  of  the 
oHE^^  mode.  The  geometry  of  the  cavity  is  the  same  as  the  one  shown 
in  Figure  V-1.  The  analysis  in  this  section  follows  very  closely  that 
in  the  previous  section,  therefore  only  the  results  will  be  given  here. 

The  power  dissipations  due  to  the  dielectric  loss  and  the  end 


walls  loss  are  respectively. 
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and 


p'  =2 
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All  the  symbols  in  the  above  expressions  have  been  defined  either  in 
Chapter  IV  or  in  Appendix  B.  The  total  time  average  energy  stored  in 
the  cavity  is  given  by 
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P^  and  W^  are  related  to  the  Q  factor  by  the  following  relati 
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0,  is  the  loss  tangent  of  the  dielectric  and  6  is  the  skin  depth 
of  the  reflecting  end  plates. 

The  expressions  C,^/C^  and  evaluated  numerically. 

Results  are  given  in  Figure  V-3  in  which  ^.p/C^  and  are  plot¬ 

ted  against  the  normalized  major  axis  (NMA)  for  various  values  of 

with  e,/f  =  2.5  .  The  characteristics  are  similar  to  those  of  tne 
1  o 
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e^lln  large  as  desired  by  choosing  suitable 

values  of  NMA.  As  NMA  approaches  Infinity,  ^T^^d  ®PPro®ches  ^-J^q 

all  .  In  the  region  where  is  leurge,  the  slopes  of  these 

cvirves  are  very  large;  in  other  words,  a  small  variation  in  NMA  can 

cause  a  rather  larger  variation  in  ®  large  variation  in 

0^  .  The  behavior  of  C'/C  is  similar  to  that  of  the  mode. 

^d  T'  w  e  Un 

Similar  deductions  as  those  given  in  section  5*1  can  be  made  and  will 
not  be  repeated  here. 

5.3  Relation  between  Q  and  a 

In  19^^  Davidson  and  Simmonds  (4l)  derived  a  relation  between  the 
Q  of  a  cavity  con^josed  of  a  uniform  transmission  line  with  short- 
circuiting  ends  and  the  attenuation  constant  a  of  such  a  transmission 
line.  Later  in  1950>  Barlow  and  Cullen  (55)  rederlved  this  relation. 
These  authors  showed  that  this  relation  is  quite  general  and  is  appli¬ 
cable  to  sirbitrary  cross-section,  uniform  metal  tube  wavegiiides.  Since 
then  one  of  the  steindard  techniques  for  the  measurement  of  the  attenua¬ 
tion  consteint  a  is  the  use  of  the  cavity  method*.  This  method  offers 
an  excellent  way  of  measuring  the  attenuation  constant  of  the  guide  when 
the  loss  is  quite  small.  Later  on  this  method  is  generalized  smd  applied 
to  open  waveguides,  such  as  the  single  wire  line,  the  dielectric  cylinder 
guide  and  associated  guides,  by  various  authors  (6,7,19,21). 

However,  it  should  be  remembered  that  the  formula  by  Davidson, 
Simmonds  and  Barlow  is  derived  under  the  assun^ition  that  there  exists  a 
single  equivalent  transmission  line  for  the  mode  under  consideration. 

*The  procedures  of  this  method  in  general  are  the  following.  Short  the 
uniform  transmission  line  under  consideration  at  both  ends  and  measvire 
the  Q  of  such  a  resonator.  From  the  knowledge  of  the  measvired  Q  and 
other  constants  such  as  the  cut-off  frequency  of  the  giiide,  the  frequency 
of  oscillation,  etc.,  it  is  an  easy  matter  to  obtedn  a  from  the  formula 
derived  by  these  authors. 
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This  a38iiB?>tion  is  triie  for  a  piare  TE,  TM  or  TEM  mode,  but  it  is  not 
clear  that  such  a  single  equivalent  transmission  line  exists  for  the 
hybrid  waves.  This  suspicion  originates  from  the  fact  that  a)  the  TE 
and  TO  waves  are  intimately  coupled  to  each  other,  and  b)  the  charac¬ 
teristic  in^jedance  defined  by  Schelkunoff  (5b)  is  not  constant  with 
respect  to  the  transverse  coordinates.  It  is,  therefore,  very  difficult 
to  conceive^  the  possibility  that  there  exists  a  single  equivalent  trans¬ 
mission  line  for  this  hybrid  mode;  at  best  the  hybrid  wave  may  be  repre¬ 
sented  by  a  set  of  transmission  lines  coupled  tightly  with  one  another. 
Hence  tne  formula  by  Davidson,  Simmonds  and  Barlow  is  not  applicable  to 
the  hybrid  wave.* 

A  more  general  relation  between  Q  and  a  can  De  obtained 
without  using  the  transmission  line  equivalent  circuit,  provided  that 
a  is  very  small  compared  with  f>  (57).  The  propagation  constant  of  a 
guided  wave  with  small  attenuation  constant  at  is 

r(a)Q)  =  Q!(co^)  +  iP(cu^)  .  (1) 

At  resonance**,  the  following  relation  is  true 

r(u)  )  + -^  AO)  iiiip(a)  )  .  (2) 

o  dtu  o 

Combining  equations  1  and  2  we  have 

“(“o>  ■  •  ‘  ■  ^ 

According  to  the  definition  of  group  velocity  v^  which  is  and 

the  definition  of  the  Q  factor  which  is  (n^/2(^)  ,  we  finally  arrive 

TfBut  several  investigators  (19,21)  apparently  unaware  of  this  restric- 
tlon,  used  this  formula  in  their  investigations  of  the  hybrid  wave. 
**The  resonant  cavity  is  made  by  shorting  both  ends  of  the  guide  under 
consideration. 
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at  the  re3.ation 


a  = 


<u 


2Qv„ 


Vg  2Q 


This  is  the  general  relation  that  we  are  seeking.  Substituting  the 

values  of  v  /v  *  for  TE,  TM  or  TEM  into  equation  4,  one  gets  the 
P  S 

relations  derived  by  Davidson,  etc.  For  the  TM  or  TE  mode, 

V 


1-i^f 

c 


a  = 


1 _ P 

c 


and  for  the  TEM  mode. 


V  /v  =  1  ,  a  =  P/2Q  .  X  is  the  cut-off  wavelength. 

PS  ^ 

The  group  and  phase  velocity  of  the  dominant  modes  can  be 

obtained  easily  from  the  co-S  diagram.  A  sketch  of  the  co-3  diagram 

for  the  dominant  modes  is  shown  in  Figujfe  V-4.  It  can  be  seen  that 


at  low  frequencies  or  small  P's,  v, 


ph 


and  again  at  very  high  fre¬ 


quencies  or  large  P's,  .  Therefore,  the  relation  a  -  S/2Q 

is  applicable  only  at  very  low  frequencies  or  at  very  high  frequencies. 


*  V  =  v  ,=  the  phase  velocity  of  the  wave, 
p  ph 


•5- 


Fig.  W-k.  A  sketch  of  the  co-P  diagram  for  the  dominant  dielectric 
rod  mode. 
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chapter  VI  -  EXPERIMEin:AL  INVESTIGATIONS 


It  is  the  purpose  of  this  chapter  to  investigate  and  verify  the 
analytic  results  experimentally.  The  properties  of  a  certain  propagating 
mode  along  an  infinitely  long  uniform  waveguide  are  usually  specified  by 
three  characteristics^^  (a)  the  guide  wavelength  which  is  directly  re¬ 
lated  to  the  propsigation  constant  of  the  wave,  (b)  the  field  configu¬ 
rations  or  the  field  distributions,  (c)  the  power  loss  or  the  attenuation 


*In  order  that  a  good  matching  condition  may  be  obtained  so  that  in 
coupling  energy  into  and  out  of  a  dielectric  rod  guide  without  the 
presence  of  high  standing-wave  ratio  which  is  a  measure  of  the  relative 
intensities  of  reflected  and  incident  waves,  the  knowledge  of  the  charac¬ 
teristic  in^iedance  or  the  wave  impedance  is  very  important.  As  we  have 
pointed  out  earlier  the  characteristic  impedance  defined  in  the  usual 
manner  (i.e.,  the  ScheiJc’anoffi definition)  is  not  meaningful,  since  it  is 
a  function  of  the  transverse  coordinate  system.  A  mean  value  impedance, 
which  takes  into  account  the  energy  distribution  over  a  cross-section  of 
the  rod  was  first  suggested  by  Wegener.  He  divided  the  (circular)  dielec¬ 
tric  rod  into  four  sections  and  in  each  of  these  four  sections  he  assumed 
the  field  to  be  independent  of  0,  the  angular  variation,  so  that  in 
regions  I  and  III,  0  is  assumed  to  be-  zero  and  in  regions  II  and  IV,  0 
is  assumed  to  be  n/2.  [See  Figure  9,  ref.  (17)].  The  approximate 
expression  for  mean  impedance  is  therefore 


(^)  (ExH*)  e  dA 
^0  0=n/2 


f(E  X  H*) 


e  dA 
— z 


He  showed  that  Z/Z^  w  X  /x  .  Similar  approximate  mean  impedance  as 
defined  above  may  be  obtained  for  the  elliptical  dielectric  rod. 

The  fact  that  the  characteristic  impedance  of  the  hybrid  waveguide  is 
not  well  defined  shows  that  the  single  transmission  line  analog  is  at 
best  an  approximation.  Any  measurements  assuming  the  single  transmission 
analog  of  this  guide  are  therefore  approximate,  and  shoiild  not  be  consi¬ 
dered  as  precision  measurements. 

Since  at  present  we  are  only  concerned  with  the  problems  of  wave  propa¬ 
gation  along  an  infinitely  long  uniform  dielectric  rod,  the  "characteris¬ 
tic  impedance"  or  the  equivalent  circuit  network  of  this  guide  will  not 
be  considered. 
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constemt  of  the  vave.  Experiments  vd.ll  therefore  he  specifically 
designed  to  measure  these  three  quantities. 

After  a  detailed  description  of  the  experimental  apparatus,  the 
methods  of  measurement  for  these  various  quantities  are  discussed.  The 
experimental  results  are  then  con^iared  vith  the  theoretical  results.  A 
discussion  vd.ll  he  given. 

6.1  Experimental  Apparattis 

Figure  VI-1  is  a  photograph  of  the  general  physical  appearance  of 
the  experimental  set-up.  A  block  diagram  is  shovm  in  Figure  VI-2.  For 
the  sake  of  convenience  and  simplicity,  measurements  were  performed  in 
the  X-hand  frequency  range .  The  microwave  X-hand  power  was  obtained 
from  an  X-13  Varian  reflex  klystron  which  offered  a  maximum  power  output 
of  five  milliwatts  and  was  powered  by  the  Hewlett-Packard  power  supply. 
The  microwave  signal  was  modulated  vdth  a  1000  cps  square  wave.  The 
output  of  the  klystron  was  connected  to  an  isolator  followed  by  an 
attenuator,  a  cavity  resonator,  a  slotted  line  section,  and  a  section 
of  standard  X-band  rectangular  me.tallic  waveguide.  These  were  steindard 
X-band  coii5)onents .  The  other  end  of  the  rectangvilar  metallic  wavegviide 
was  connected  to  the  special  apparatus  specifically  designed  for  the 
present  experimental  investigation^  see  Figure  VI-2. 

The  follovd.ng  sections  are  devoted  to  a  detailed  description  of 
the  special  apparatus. 

A.  The  Launching  Device 

The  method  of  transferring  microwave  energy  from  an  ordinary  metal¬ 
lic  waveguide  into  a  dielectric  rod  was  not  very  difficult  or  complicated. 
Since  a  rectangiilar  metal  guide  operating  in  the  dominant  TE^^q  mode  had 


Fig.  VI -1.  Experimental  Apparatus 
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an  electric  field  whose  conflgviration  was  roughly  similar  to  the  trans¬ 
verse  coiqponent  of  the  electric  field  of  the  mode  or  the  ^HE^i ^ 

mode  on  the  dielectric  guide,  the  transfer  could  be  made  slnqply  by 
Inserting  the  dielectric  rod  longitudinally  into  the  metal  guide  for  a 
short  distance.  The  orientation  of  the  cross-section  depended  upon 

whether  the  mode  or  the  mode  was  desired.  To  Iviprove  the 

e  11  o  u 

matching  and  to  minimize  reflection  the  dielectric  rod  was  tapered  to  a 
point  within  the  guide  and  after  emerging  from  the  metal  guide  the  rod 
was  tapered  to  whatever  size  was  required  for  a  given  test.  Furthermore 
a  flare  pyramidal  horn  whose  flare  angle  was  adjusted  for  best  energy 
transfer  was  connected  to  the  rectangular  metal  guide.  (See  Figure  VI- 3) • 


Figure  VI- 3.  The  Launching  Device. 
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B.  The  Elliptical  Dielectric  Rod 

Slace  dielectric  rods  of  elliptical  cross-section  were  not  com¬ 
mercially  available,  they  vere  machined  from  available  rectangular 
Incite  strips  which  were  at  least  five  and  half  feet  long.  A  total 
of  twelve  rods  of  different  sizes  and  eUlptlcltles  were  made,  in 
order  that  the  e:q>erimental  data  would  cover  a  wide  range  of 
Emd  values.  A  picture  of  these  rods  is  shown  in  Figure  VI-4.  One 
end  of  each  rod  was  machined  very  flat  while  the  other  end  was  tapered 
as  described  in  section  6.1A  to  fit  into  the  meteO.  guide.  A  small  chunk 
of  Incite  was  taken  from  each  rod  in  order  to  measure  the  electxrical 

e 

properties  of  each  rod  individxially  by  Von  Hippel’s  method  (58)*  It  was 
fomd  that  the  dielectric  constant  of  these  rods  veurled  between  e  =  2.5 
to  e  =  2.6  and  the  loss  tangent  varied  from  tem  5  s  0.005  to 
tan  8  =  O.D03.  It  should  be  noted  that  due  to  the  resilient  property 
of  Incite  it  was  very  difficult  to  machine  such  a  reqviired  length  uni¬ 
formly.  A  special  and  rather  expensive  technique  was  developed  and  used. 
Although  extreme  cem^  was  taken  in  making  these  rods,  some  small  non- 
uniformities  which  might  attribute  to  experlmentEG.  errors  were  unavoid¬ 
able.  The  major  axis  and  of  these  elliptical  rods  ranged  from 
2A  =  1.5  in.  to  2A  =  0  S  in.  and  =  00  to  =0.3?  . 

C.  The  Shorting  Plate 

In  order  that  the  plate  could  be  a  good  shorting  device,  its  sur¬ 
face  had  to  be  very  flat  and  large  enough  to  Intercept  practically  all 
of  the  energy  outside  the  dielectric  rod  and  the  plate  had  to  be  made 
of  good  conducting  material.  A  l/4"x36"x36"  sdumlnum  plate  was  used. 
One  side  of  the  nmi  mun  plate  was  machined  flat  and  its  surface  was 
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cleaned  and  polished  In  order  to  assure  meuclmum  conductivity.  The 
plate  and  Its  support  can  he  seen  clearly  In  Figure  VI -1. 

D.  The  Probe  emd  Its  Carriage 

To  detect  and  measure  the  electromagnetic  field  on  the  dielectric 
rod,  a  small  electric  probe  was  designed  and  used.  The  probe  con¬ 
sisted  of  a  section  of  rigid  coaxial  cable  whose  outside  diameter  was 
about  1/8"  and  whose  length  was  about  1'3".  An  Inch  from  one  end  of 
the  cable  was  formed  Into  a  gradual  90°  bend  and  the  center  conductor 
protrxjded  about  1/8" .  The  other  end  was  connected  to  a  crystal  detec¬ 
tor  which  was  calibrated,  and  the  output  of  this  detector  was  connected 
to  the  HP  standing  wave  Indicator.  The  L  bend  was  Introduced  to 
reduce  the  amowt  of  metal  conductors  parallel  to  the  electric  field 
Indicator.  The  probe  and  the  detector  were  supported  by  a  stand  which 
was  fastened  to  a  HP  ceurrlage.  A  picture  of  the  probe  and  Its  s\^port 
Is  shown  In  Figure  VI-^.  The  whole  Instrument  was  so  designed  that  the 
probe  might  be  moved  and  down  radleiUy  with  respect  to  the  center 
axis  of  the  dielectric  rod  and  longitudinally  along  the  center  axis  of 
the  rod.  Furthermore,  the  probe  could  be  adjusted  to  detect  either 
Ej  or  E^  field.  The  longitudinal  movement  of  the  probe  could  be 
measured  from  a  scale  on  the  carriage;  and  a  dial  Indicator  was  used 
to  obtain  accurate  measurements  of  small  longitudinal  movements  of  the 
probe.  The  radial  movement  of  the  probe  was  measured  by  a  level  tele¬ 
scope  whose  movement  had  been  calibrated. 

With  the  help  of  a  transit  and  a  level  this  whole  esgperlmented 
set-v^  was  aligned  carefully.  The  dielectric  rod  had  to  be  very  straight 
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and  Its  axis  perpendicular  to  the  shorting  plate.  To  Insure  a  good  con¬ 
tact  between  the  polished  shorting  plate  and  the  flat  end  of  the 
dielectric  rod,  a  slight  pressure  was  asserted  on  both  ends  of  the  system. 
To  minimize  sagging  of  some  small  or  flat  rods,  very  thin  nylon  threads 
were  used  along  the  rod  to  provide  support.  Although  disturbances  due 
to  these  threads  were  vtnavoldable,  because  of  the  sizes  of  the  rods  used 
at  this  frequency  range  very  little  perturbation  was  observed. 

6 . 2  Method  of  Measurement 

In  general  there  are  two  most  commonly  used  methods  for  measiirlng 

the  propagation  characteristics  of  a  certain  mode  along  a  \iniform  low 

loss  waveguide.  The  first  one  Is  the  so-called  resonator  technique.  The 

guide  under  consideration  is  placed  between  two  peurallel  plates  with 

proper  coi^pllng  holes.  Resonance  occurs  when  the  length  of  the  cavity  is 

n  where  n  is  an  integer  eind  is  the  guide  wavelength  of  the  mode 

on  the  guide.  X  can  be  measured  easily  either  by  counting  the  number 

g 

of  minima  within  the  cavity  length  with  a  small  probe  or  by  moving  one 
of  the  reflector  plates  and  measuring  the  displacement  of  the  plate  for 
each  resoneut  peak.  By  measuring  the  Q  of  this  cavity,  the  attenuation 
constant  a  can  readily  be  calculated*.  This  method  is  psu:tic\ilarly 
useful  and  accurate  for  very  low  loss  transmission  lines'""'*’.  The  second 
method  is  the  standing  wave  measurement  technique.  The  guide  is  terminated 
.  by  a  perfectly  reflecting  plate  acting  as  a  short- circ\ilt  device.  The 
propagating  wave  Is  perfectly  reflected  by  the  termination  and  a  standing 

*See  equations  5-3-4 

**Thls  method  was  first  used  by  Chandler  (19)  on  the  measurement  of  attenua¬ 
tion  factor  for  the  HEn  mode  on  a  very  small  circular  dielectric  rod. 
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vave  pattern  is  fonoed  along  the  guide.  A  probe,  such  as  the  one  des¬ 
cribed  earlier  in  section  6. ID,  can  be  used  to  detect  the  microwave 
signal  along  an  open  waveguide.  By  measuring  the  distance  between  two 
adjacent  minima  of  the  standing  wave  pattern,  and  the  standing  wave 
ratio.  It  Is  an  easy  matter  to  calculate  the  guide  wavelength  and  the 
attenuation  factor  of  the  mode  on  the  guide'll. 


^*The  formula  relating  the  attenuation  factor  a  with  the  standing  wave 
ratio  can  be  derived  as  follows:  It  Is  well  known  that 

*  "  5  log^o  ^  db 

where  and  are  respectively  the  input  and  reflected  power  of 

P 

*  r  - 1  2 

—  =  — J-)  ,  where  r  is  the  standing  wave  ratio 

3  r  —  1  2  ‘ 

at  the  probe.  Therefore,  we  have  a  =  —  log  (— ^)  db/m  in  which  i 

JL  xv  IT  ^  X 

is  the  length  of  the  guide  as  indicated  In  Figure  VI-6. 


the  guide;  and 


elliptical  dielectric  rod 


^reflecting  plate 


a 


Figure  VI -6. 

To  take  into  account  the  loss  due  to  Imperfection  of  the  shorting  plate 
one  notes  that  the  attenuation  measured  at  point  a  Is  A^  =  ai^  +  T)  , 
suid  similarly  the  measured  attenuation  at  point  b  Is 
A^  =  al^  +  D  where  D  is  the  loss  of  the  x'cflectlng  plate.  Combining 
these  two  equations  and  eliminating  D  one  gets 

(A  -  A.  ) 

a  =  — 2 -  (db/meter) 

<'a  - 
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The  latter  method  was  used  for  our  measurements  since  It  presented 
a  simple  and  expedient  way  of  measuring  the  desired  quantities  with 
reasonably  good  accuracy.  To  avoid  perturbation  by  the  launching  device 
or  by  end  effects^  measurements  were  made  In  the  middle  section  of  the 
rod.  It  should  also  be  mentioned  that  throughout  this  whole  experiment 
the  coupling  between  the  probe  and  the  field  was  kept  at  a  minimum  In 
order  to  avoid  Interference  with  the  propagating  wave. 

6.3  Comparison  of  the  Theoretical  and  Experimental  Results 

The  resiilts  are  separated  Into  three  general  categories. 

A.  Gvilde  Wavelength 

Guide  wavelength  was  measured  according  to  the  procedures  described 
eeirller.  The  distance  between  adjacent  minima  of  the  stemdlng  wave  pat¬ 
tern  was  measured  at  several  sections  along  the  guide  and  the  average 
value  was  calcvilated  and  recorded  as  the  measxured  .  The  maximum 

difference  between  these  measurements  was  about  Wavelength  measure¬ 

ments  were  taken  from  nine  different  sizes  of  elliptical  dielectric  rod 

for  the  mode  and  the  mode.  Normalized  experimental 

e  11  o  11 

results,  together  with  their  corresponding  theoretical  results,  are 
given  in  Figures  VI-7  through  VI-15 •  The  physical  size  of  each  dielec¬ 
tric  rod  used  and  Its  measured  dielectric  constant  are  Indicated  In  each 
figure.  Excellent  agreement  was  obtained. 

To  Illustrate  the  agreement  between  the  analytic  and  experlmenteG. 
results,  we  introduce  the  following  table: 
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(circles  are  esqperlaental  points) 


2q  cosh 


(Circles  are  esqperimental  points) 


(Circles  are  ejq>erlnental  points ) 
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Mode  of 

Size 

^0 

2q  cosh  1 

(^) 

0  meas. 

(^) 

0  calc. 

c 

Operation 

(2Ax2B) 

X 

o 

0 

0.770"x0.645" 

1.22 

0.55 

0.845 

0.853 

2.55 

0.63 

0.80 

0.804 

o®®ll 

0.770"x0.645” 

1.22 

0.555 

0.83 

0.833 

2.55 

• 

0.635 

0.784 

0.787 

e  11 

0.769"x0.505" 

0.775 

0.56 

0.925 

0.925 

2.5 

0.785 

0.80 

0.802 

o  11 

0.769"x0.505” 

0.775 

0.55 

0.875 

0.882 

2.5 

0.695 

0.80 

0.80 

e  11 

1.005"x0.361" 

0.376 

0.725 

0.962 

0.964 

2.5 

1.00 

0.862 

0.862 

™11^ 

0  11 

1.005"x0.36r' 

0.376 

0.72 

0.87 

0.873 

2.5 

0.93 

0.79 

0.792 

It  was  found  that  the  wavelength  measurenents  were  rather  Insensi¬ 
tive  to  small  non-uniformity  of  the  rods  and  to  the  variation  of 
humidity  and  teii5)eratTire  in  the  laboratory.  Incidentally,  the  above 
e}q>erlment  also  suggested  a  rather  convenient  way  of  measuring  the 
dielectric  constant  of  a  certain  low  loss  dielectric  material. 

B.  The  Field  Distributions 

In  order  to  establish  the  degree  of  field  purity  an  examination  of 
the  radial  field  decay  at  a  fixed  axial  position  was  carried  out.  The 
decay  of  the  axial  electric  field  was  measured  since  it  can  most  easily 
be  detected  by  a  probe  pointed  in  the  axleil  direction.  For  maximum  sig¬ 
nal  strength,  the  probe  was  aligned  in  the  ij  =  x/2  plane  for  the 
mode  and  tj  •  0  for  the  mode.  The  general  method  of  measurement 
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has  been  outlined  In  section  6.2. 

Esqperlmental  results,  together  with  their  theoretical  results 
for  the  mode  are  shown  in  Figures  VI-16  through  VI-I8.  Six 

different  rods  ranging  from  =  00  to  =  0.376  were  used.  It_ 
can  be  seen  that  the  e:q>erlinental  results  corresponded  rather  well 
with  the  analytical  results.  The  largest  differences  were  found  among 
thin  rods.  This  effect  may  be  e^q^lained  by  the  fact  that  for  small 
valiies  of  NMA,  a  large  percenteige  of  energy  was  carried  outside  the 
dielectric  rod,  thvis  a  small  amoxint  of  c\irvature  or  sagging  may  have 
caused  some  errors  in  the  field  decay  measurements.  These  measurements, 
together  with  the  wavelength  measurements  verified  the  existence  of 
the  mode  along  an  elliptical  dielectric  rod. 

Similar  measurements  were  performed  for  the  mode.  Four 

rods  ranging  from  |  =  00  to  (  »  0.376  were  used.  Results  are 

o  0 

shown  in  Figures  VI-19  through  VI-20.  Again,  good  agreement  with 

theoretical  results  were  observed.  These  measurements  also  confirmed 

the  existence  of  the  HE,^^  mode. 

o  11 

The  above  discussion  shows  clearly  the  necessity  of  having  a 
structure  which  may  support  the  dielectric  rod  and  at  the  same  time 
will  not  interfere  with  the  desired  propagating  mode.  One  of  the  best 
ideas,  which  was  first  proposed  by  D.  D.  King  (20),  is  the  use  of  the 
image  plane.  He  took  eidvantage  of  the  symme triced,  property  of  the 
HE^  mode  and  mounted  a  haLLf -round  dielectric  rod  on  an  image  plane. 

It  can  be  seen  that  his  idea  can  very  well  be  extended  to  the  ellipti¬ 
cal  dielectric  rod.  This  image  plane  can  not  only  serve  as  a  support 
without  disturbing  the  fields,  but  edao  may  serve  as  a  poleurlzation 
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C.  Attenuation  Constemts 


Attenuation  measurements  were  made  by  probe,  using  the  steudard 
standing  wave  technlqvie  (39)*  For  each  experimental  point  two  mea¬ 
surements  at  two  different  locations  were  CEurrled  out  In  order  that 
the  loss  due  to  the  Imperfections  of  the  terminating  reflecting  plate 
could  be  eliminated,  (See  the  footnote  on  p.  155)*  The  percentage 
variations  of  the  loss  factor  of  these  dielectric  rods  were  found  to 
be  quite  large.  The  loss  factors  for  various  rods  were  found  to  vary 
from  tan8  =  O.OO5  to  tan6  =  O.OO3  •  Both  theoretical  and  experi¬ 
mental  results  for  the  and  the  waves  are  shown  In 

e  11  o  11 

Figures  VI-21  through  VI-27.  In  general,  the  agreement  Is  quite  good, 
and  It  Is  better  at  higher  frequencies  than  at  lower  f reqviencles .  One 
of  the  reasons  for  this  is  that  at  lower  frequencies,  more  energy  is 
distributed  outside  the  rod;  thus,  more  energy  Is  radiated  due  to  the 
slight  curvatiure  of  this  open  guide.  Furthermore,  the  disturbance  of 
the  field  caused  by  the  presence  of  the  siq)porting  threads  and  the 
probe  Is  more  pronounced  at  lower  f reqxxencies .  Since  the  attenuation 
Is  lower  at  lower  frequency,  the  standing  wave  ratio  Is  higher  and  the 
percentage  error  In  the* measurements  of  this  high  standing  wave  ratio 
Is  therefore  leurger.  Because  the  attenuation  constant  Is  a  meastire  of 
the  power  loss  as  coo^>ared  with  the  power  transmitted.  It  Is  q\ilte 
understandable  that  the  above  mentioned  factors  would  affect  the  ac¬ 
curacy  of  our  measxirements  more  at  lower  frequencies.  It  Is  for  this 
reason  that  the  resonator  method  Is  superior  for  low  attenuation  mea¬ 
surements.  At  higher  frequencies  most  of  the  energy  Is  carried  Inside 
the  guide;  very  little  disturbEuice  will  result  from  the  supporting 
threeids,  the  probe,  and  the  small  curvatures  of  the  rod.  The  accuracy 
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of  the  experimental  measurements  Is  thus  greatly  Istproved.  Another 
source  of  error  is  probably  due  to  the  approximations  used  in  confuting 
the  numerical  results  from  the  analytical  equations.  It  was  mentioned 
earlier  in  Chapter  IV  that  as  the  elliptical  cross  sections  become 
flatter,  l.e.,  gets  smaller,  more  terms  of  the  expetnslons  are 

required  to  obtain  more  acc\arate  numerical  results.  However,  s\iffl- 
clently  close  agreement  between  the  emalytical  and  experimental  results 
is  observed  to  warrant  the  verification  of  the  theoretical  predictions. 

It  is  noted  that  a  dielectric  ribbon  having  the  same  cross  sec¬ 
tional  area  as  a  circular  dielectric  rod  and  operating  In  the  domlnEuit 
mode  does  Indeed  offer  much  less  attenuation  than  the  circxilar 

e  11 

dielectric  cylinder. 

6.4  Conclvisions 

Special  e;q>erlmental  apparatus  was  designed  to  measure  the  propa¬ 
gation  characteristics  of  the  two  principal  dominant  modes,  namely  the 

and  the  modes.  Despite  the  mentioned  sources  of  eOTeri- 

e  11  o  11 

mental  errors,  measvired  results  were  found  to  be  in  good  agreement  with 
the  calculated  results.  Parti c\ilarly  good  agreement  was  observed  in 
the  gvilde  wavelength  category.  The  existence  of  these  two  dominant 
modes  was  verified. 

It  was  noted  experimentally  that  the  next  higher  order  mode  occurs 
at  a  higher  frequency  for  the  flatter  elliptical  cross-section  rod. 
Therefore,  the  flatter  rod  possessed  not  only  a  lower  attenuation  factor 
but  also  a  wider  bandwidth. 
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CHAPTER  VII  -  SUMMARY  AMD  COKCUJSIONS 

The  problem  of  electromagnetic  wave  propeigatlon  along  an  elliptical 
dielectric  rod  was  considered.  It  was  shown  that,  in  general,  no  pure 
TB  or  wave  might  exist  on  such  a  waveguide  except  when  the  eccentricity 
of  the  rod  was  zero.  In  order  to  satisfy  the  boundary  conditions,  an 
infinite  series  of  product  terms  of  Mathieu  and  modified  Mathieu  func¬ 
tions  were  used  to  represent  the  field  configurations  in  one  of  the  two 
regions.  The  field  coinponents  and  the  characteristic  equations  of  four 
principal  types  of  waves  were  obtained.  They  were  classified  as  the 

node,  the  node,  the  HE^^^  node,  euad  the  HE^^^  node.  These 

emn  emn  omn  onn 

nodes  were  degenerate  and  became  the  well  known  HE  mode  'vAien  the  eccen- 

mn 

trlclty  was  zero.  The  mod^  with  m  =  1  and  n  =  1,  called  the  dominant 
modes,  possessed  no  cut-off  frequency.  The  propagation  characteristics 
of  the  dominant  modes  were  considered  in  detail  both  analytically  and 
experimentally.  Extensive  numerical  coD^)utatlons  on  the  properties  of 
the  guide  wavelength,  the  rate  of  field  decay,  and  the  attenuation  charac¬ 
teristics  of  the  domlneuit  modes  were  carried  out. 

E:q)erlments  were  designed  and  performed  using  various  sizes  of  elllp- 
ticsJ.  luclte  rod  to  verify  the  analytic  results  of  the  dominant  modes. 
Meeisured  data  were  co]iq)ared  with  theoretical  results  and  it  was  shown 
that  very  good  agreement  was  obtained. 

The  Q  of  a  dielectric  cavity  resonator  operating  in  either  one  of 
the  dominant  modes  was  also  coiqputed.  It  was  found  that  a  very  high  Q 
cavity  could  be  made  using  proper  size  dielectric  strip. 

It  was  demonstrated  aneilytlcally  and  experimentally  that  a  thin 
elliptical  rod  operating  in  the  gHE^^^  mode  typically  had  considerably 
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lower  loss  than  did  a  circular  rod,  operating  In  the  mode,  having 
the  sane  cross-sectional  area. 


Other  advantages  are  listed  below: 


a.  A  flat  elliptical  dielectric  rod  has  larger  surface  area, 
thus  it  would  be  easier  to  handle  at  very  high  frequencies, 
svich  as  in  the  mm  wavelength  range. 

b.  The  mode  possesses  greater  bandwidth,  since  the  cut¬ 

off  frequency  of  the  next  higher  order  mode  is  higher  for 
flatter  eUlptical  cross  sections. 

c.  Depolarization  effects  are  minimized  because  the  guide  wave¬ 
lengths  differ  for  the  even  and  odd  nodes.  It  is  known  that 
internal  stredn,  non-uniform  dimensions,  and  bends  of  the 
circular  rod  cause  the  mode  to  change  polarization. 

d.  A  flat  elliptical,  dielectric  rod  which  may  be  approximated 
by  a  strip  is  easier  to  fabricate  since  wider  dimenslonad. 
tolerances  are  permitted.  This  is  because  the  guide  wave¬ 
length  and  the  attenuation  constant  sure  slower  vau:7lng 
functions  of  the  dimensions  for  the  strip  than  for  the 

i'od. 

It  shovild  sQso  be» noted  that  the  mode  can  be  launched  as  easily 

as  the  HE^  mode. 


The  advantages  and  disadvantages  of  using  the  dielectric  tape  line 


as  a  transmission  line  in  coo^arison  with  the  conventional  metal  tube 
waveguide  (at  frequencies  above  30  kmc)  eu:e  discussed  below. 


a.  Keeping  the  spread  of  the  field  outside  the  dielectric  tape 
within  a  reasonable  distance  from  the  surface  of  the  dielec¬ 
tric  guide,  the  loss  factor  of  the  gHE^^^  mode  on  a  dielectric 
tape  line  can  be  made  smaller  than  that  of  the  dominant  mode 
in  a  rectanguJ.ar  metal  tiibe  guide,  but  still  somewhat  lEurger 
than  that  of  the  mode  in  a  clrculEO*  metal  tube  guide. 


b.  In  order  to  avoid  mode  conversion  the  surface  of  the  cir¬ 
cular  metal  waveguide  must  meet  very  close  tolerances. 

These  kind  of  tolerances  are  not  required  for  the  elliptical 

•  dielectric  rod;  thus,  11}  is  easier  to  fabricate. 

c.  Unlike  the  metal  tube  waveguide,  the  field  of  a  dielectric 
tape  is  not  entire^  confined  within  the  strip.  Consequently 
it  can  be  subjected  to  Interference  due  to  nearby  foreign 
objects  or  foreign  signals.  The  presence  of  curvature  or 
discontinvilty  of  the  strip  will  also  cause  energy  loss  by 
radiation.  It  is  rather  difficult  to  8tq)port  the  dielectric 
rod  without  dlstturblng  the  field.  The  vise  of  image  plane, 

as  proposed  by  Klisg,  as  a  supporting  device  appeeurs  to  be 
quite  suitable. 

The  greatest  attraction  of  a  surface  wave  dielectric  tape  line  as  a 
millimeter  wave  transmission  line  is  in  its  simplicity  of  construction, 
its  low  cost  and  ease  of  meuaufacture,  and  its  flexibility. 

The  analytic  method  of  solving  this  elliptical  dielectric  wave¬ 
guide  problem  should  prove  applicable  to  others  involving  the  use  of 
the  elliptical  coordinate  system  and  the  Mathleu  functions.  For 
Instance,  one  may  apply  this  technique  to  the  problemsobf 

a.  the  propagation  of  electromagnetic  waves  along  an  ellip¬ 
tical  dielectric  tube,  or 

b.  the  surface  wave  propagation  along  an  elliptical  Sommer- 
feld  or  Goubau  wire,  or 

c.  the  scattering  of  electromagnetic  wave  by  an  elliptical 
dielectric  cylinder. 

Other  analogous  mechanical  problems  can  also  be  solved  in  a  similar 
manner. 

We  did  not  consider  the  source-present  problem.  However,  it  is 
emphasized  here  that  it  is  not  possible  to  express  any  arbitrary  field 
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dlstribution  In  terms  of  the  propagating  modes  alone  for  any  open 
boundary  problems.  The  problem  with  source  present  must  then  be  for¬ 
mulated  In  terms  of  Green's  function  In  the  form  of  a  contoxir  Integral. 
The  residues  at  the  poles  of  the  Integrwd  will  give  rise  to  modal 
type  waves  which  are  eilso  called  guided  waves.  The  contribution  of 
the  Integral  around  the  branch  cut  will  give  rise  to  a  radiated  wave. 


APPENDIX  A 


Mathematical  Relations 

Many  of  the  formulas  ^and  eiqpresslons  which  are  used  throughout 
this  report  are  given  in  the  following  sections  for  convenient  refer¬ 
ence.  Those  relationships  wUch  are  considered  well  known  are  stated 
with  only  a  reference  to  their  origin.  Others,  which  are  considered  not 
so  well  known  are  discussed  in  more  detail.  Some  relations  which  are 
given  nere  for  the  first  time  are  derived. 


A.l  Series  Representations  of  Mathleu  and  Modified  Mathleu  Functions 
The  Mathleu  differential  equation  may  be  written  In  the  form 


2 

+  (c  -  2-ir^co8  2tj)@  =  0 

dTj 


(A.1-1) 


where  y  is  a  consteint  which  may  be  positive  or  negative,  and  c  Is 
the  sepeuration  constant  or  the  characteristic  number.  The  periodic 
solutions  of  A.1-1  which  may  be  expanded  in  terms  of  trigonometric 
functions  are  given  below  (U2,U5): 

For  ^  0  , 


/  2\  S  .(2n) 

ce2^(il,r  )  =  ^  t^2T 
r  =  0 

2rTj 

(A.1-2) 

°®2n+l^’*^^  ^  "  E  ^+1 

r  =0 

cos(2r+l)Tj 

(A.1-3) 

sin(2r+l)T] 

(A.1-4) 

2,  S  ,(2n+2) 

®®2n+2^’*^^  ^  ■  L.  ®2r+2 
r  =  0 

sin(2r+2)  tj 

(A.1-5) 
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cx> 


1)  =  E 


r=  0 


cos  2rT) 


(A.  1-6) 


“  E  (-1)^'^“  co8(2r+l)T)  (A.1-7) 

r  =  0 

se*^^l(Ti,|r^|)  =  Y,  (-1)^^“  4^1^^  siii(2r+l)Tj  (A.1-8) 

r  =0 

*  8e|n+2^nJr^l)  =  f  (-!)''■'“  8in(2r+2)Ti  (A.1-9) 

r  =  0 

The  e^qpanslon  coefficients,  A^^,  ®^2^^  “‘® 

I  2 1 

fiinctlons  of  |r  I  and  have  heen  tabulated  by  NBS  (^9)  for  various 
1  21  I  21 

values  of  |  r  I  to  1  r  I  =  25  . 

The  modified  Hathleu  differential  eqxiatlon  be  vrltten  in  the 


form 


9 

~  -  (c  -  2r^co8h  2|)R  =  0 

dr 


(A.1-10) 


2 

where  y  is  a  constant  which  may  be  positive  or  negative  and  c  is 
the  separation  constant  or  the  cheuracterlstlc  number.  The  stable  solu¬ 
tions  of  A.1-10  which  correspond  to  the  periodic  solutions  of  A. 1-1  can 
be  expressed  in  terms  of  Bessel  function  product  series.  CompEurlng  this 
with  other  ways  of  ei^resslng  the  solutions  of  A.1-10,  the  Bessel  func¬ 
tion  product  series  converge  the  fastest  and  therefore  are  best  suited 
for  computational  purposes  (45). 

Let  u  =  It  I  e”^  and  v  =  |y|  e^  j  the  set  of  stable  solutions 
of  A.1-10  with  i  0  is  (45), 


r  (2n) 
^r 


Jr(^)Ji.(v) 


(A. 1-11) 
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(A.1-12) 
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(A.1-14) 


And,  the  set  of  stable  solutions  of  A. 1-10  with  f  ^  0  is  (45), 
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0 
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1  (A. 1-16) 
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panslon  coefficients  A^^^ ,  } 

bin*  ooemclnnts  - 

suid  have  etll  been  tabulated  in  a  table  prepared  by  NBS  (49).  The 
1  are  absolutely  and  uniformly  convergent.  The  analytic 
connecting  coefficients  are 

Ben  ■  4^^  (A.1-19) 
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“2q+1  ■  /  If  (A.l-21) 

=2a*2  =  “k*2‘°>’-">«  W5  ” 't")  /  ■  (A-1-22) 


The  prime  signifies  the  derivative  of  the  function  with  respect  to  t)  . 
The  normalizations  Introduced  by  Goldstein  (60)  will  be  used. 


They  are 

(A.l-23) 

2 

It  is  also  defined  that  when  T  =  0 

[*£f 

[bL^^T  -  [BLlf>]^  -  1 

when  r.s  n 

=  0 
when  T  ^  n. 

(A.  1-25) 


A. 2  Approximate  E3q)ression8  for  the  Modified  Mathleu  Functions  Suitable 
for  Small  Values  of  |  r^l 

Mathleu  (61)  first  derived  the  expressions  for  the  expansion  coeffi¬ 
cients  A^J  4^2^^  I  r^l  is  small.  They  are 


m+2r  \rt-2r  '  '  r!  (m+r)  {  'TT^  * 

Md  b(“>  -aW  , 

m-2r  m-2r  r!(m  -l) !  ^  4 '  * 


(r  ^  0,  m  >  O) 


(A.2-1) 
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vhere  »  means  approximate  equality. 

a.  Approximate  Expressions  for  1 1 )  8^  ^^2n+l^  Ir^l ) 
when  ly^l  — >  0  .  The  approximate  expression  for  Pekg^^^C 1  TT^j )  from 
A. 1-16  Is  • 


(-1)^ 


2n+l 


{■■■  *®(£li)-2 
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(n  i  1)  (A.2-3) 

(A.2-4) 

(A.2-5) 

(A. 2-6) 
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0{x  )  means  that  the  next  term  In  the  expemslon  Is  of  the  order  of  y  . 


The  value  of  0{x  )  In  A.2-5  can  easily  be  obtained  from  the  normallza- 
tlon  f  Thus  . 
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(A.2.7) 


The  following  approximations  for  the  modified  Bessel  functions  (62) 
will  be  used,  when  u  and  v  are  small: 
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and 


Kjv).  ,0  (-D^^yv)  M^)  (A.2.9) 
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(A.2-10) 


where  a  Is  Euler's  constant,  a  =  0.3772  •••  . 

Substituting  these  expressions,  A. 2-4  and  At2-6  throxigh  A.2-10 
Into  equation  A. 2-3  and  remembering  that  u  =  ye"^  >  v  =  ye^,  after 
considerable  eOgebredc  manlpxilatlons  one  obtedns  the  approximate  es^res- 
slon 


Fek^.,(|, 


2n+l' 


2n+l 
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r  (-2n+l)6 
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(A.2-11) 


Putting  the  approximate  expression  for  into  A.2-11,  one 

finally  arrives  at  the  equation 


(-1)“8 


2n+l 


^®^2n+i^  I  I )  =  ^L+y^+i)i  •  sIcLi)  t 


+  2e"^^  -  ne"^^]  +  0(r^)|  .  (n  ^  1)  (A.2-12) 


For  the  special  case,  n  =  0 


Pek^(|,|r^|)  =  +  I|-  jn(^  ye^)  -  Ig— ^  +  0(y3) 
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The  approximate  ejqpression  for  the  C}ek2j^^2^(  |,lr^  | )  can  be  obtained 
In  a  similar  manner.  Carrying  out  the  tedious  algebraic  steps,  one 
arrives  at  the  equation  for  n  ^  1 

■’1 


+  0(r  )  >  (A.2-14) 


and  for  n  =  0 


^  Gek^(|,|r^l)  =  ^  n(~  re^)ri+e"^^]-  ■  +  O(r^)  .  (A.2 

^1  re  ^  ^ 


•15) 


o  2 

b.  Approximate  Expressions  for  ^  and  Se2^^^(t,r  ) 

12 1  2  2 
_ r  l->0  «  It  is  possible  to  expand  ^  ^®2n+l^^'^  ^ 

♦ 

in  a  series  of  Bessel  functions  (45).  They  are 


and 
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These  series  are  absolutely  and  uniformly  convergent  (45).  Substituting 
eqxistlons  A. 2-1  and  A. 2-2  and  the  small  argument  approximation  of  the 
Bessel  function  into  equations  A.2-16  and  A.2-17,  and  keeping  the  first 
order  approximation,  ■we  get 
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Taking  the  derivative  of  equations  A. 2-12  through  A. 2-15  with  res¬ 
pect  to  I  ,  comhinlng  these  derivatives  with  the  required  fvmctlons, 
and  keeping  only  the  second  order  term  in  the  approximation,  we  arrive 
at  the  following  expressions: 


-  "w  *=  -(2n+l)  -- — ; - -  ■{(n+l)+  ne  > 

^ek2n+i(5^lr  I)  4n(n+l)  [  J 


+  ne‘^^)  +0(rS  >  (n  ^  1) 


(A.2-20) 


Fek|(5,|r^l)  2  25  a  p  _2  1 

^ =  -1  +  in(^  re^)  [3-26  +  0(r  )  ,  (n  ^  1) 

Fek^(5,lr^|)  2  2  L  J 


(A.2-21) 


C  Ltl  k 


(A. 2-22) 


Gek^(6jr  1)  2  2|  a  .  r  _2n  k 

■  ■■  .  g-  =  -1  +  ^  in(^  re^)  I  3  +  2e  +  0(/)  ,  (n  5.  1)  . 

Gek,(|,  r  )  *■  ^ 

^  (A.2-23) 

Taking  the  derivative  of  equations  A.2-18  and  A.2-19  with  respect  to  |  , 
combining  these  derivatives  with  the  required  functions,  and  keeping  only 


the  first  order  term  in  the  approximation,  we  arrive  at  the  following 
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ei^resslons: 


Ce* 

^  - - 5-  =  tanh  g  { 

ceg^^id^r  ) 


=  coth  g  / 


t  (-1)' 

r=:0 

(n+r)  1 

(cosh  g)^^ 

(n-r)  I 

4“"^(2r)I 

“  r 

(n+r)  I 

(cosh  g)^ 

(n-r)  1 

1— 1 
+ 

CVI 

1 

13 

r  =  0 

(n+r)  I 

(sinh  g)^ 

(n-r) ! 

4’^"'^(2r)l 

00 

E  ■ 

r  =  0 

(n+r) ! 

(sinh  g)^ 

(n-r)l 

^““^■(ar+i)  I 

)  +  o(r  )  , 

(n  ^  0)  (A.2-24) 


)+  o(r^)  , 

(n  ^  0)  (A.2-25) 


It  fihoiild  be  noted  that  all  the  approximate  expressions  derived  In 
this  section  reduce  to  the  well  known  approximate  expressions  for  the 
Bessel  functions  when  the  ellipse  degenerates  to  a  circle. 

Similar  approximate  expressions  for  the  even  order  modified  Mathleu 
functions  can  he  obtained  in  the  same  manner. 


A.  3  Degenerate  Forms  of  Mathieu  tuad  Modified  Mathieu  Functions 

When  the  ellipse  tends  to  a  circle,  i.e.,  as  the  semlfocal  lengths 
q  and  g  tend  to  zero  and  infinity  respectively,  such  that 

q  cosh  g  «  q  sinh  g  — >  r  , 

where  r  is  the  radial  congwnent  of  the  circle,  all  A^°^  and  tend 

P  P 

to  zero  except  that  =  B^°^^  ->  1  .  Therefore  the  degenerate  forms  of 

m  m 

the  Mathieu  functions  are  (4^) 

2  ■ 

cem(tl,r  ) ->  cos  mil  (A. 3-1) 

sejjj(Ti,r^) sin  mil  (A. 3-2) 


-188- 


with  m  ^  1  ,  and  when  m  =  0  . 

ce^(Ti;r^)  - •  (A. 3-3) 

The  degenerate  forms 'of  the  modified  Mathieu  functions  are 


ce^^Cs.r")  -  (-1)“  p^J2„(x  ^)  , 

o 

(A.3-4) 

^  P2n+l’^2n+l^^  ^ 

o 

(A.3-5) 

^  ®2n+l'^2n+l^^  ^ 

0 

(A.3-6) 

®®2n+2^  ^  ®2n+2'^2n+2^*  r“^  > 

o 

(A.3-7) 

Oe'  (!,r^)  -  (-1)“  P2^(x  f-)j;  (x  f)  , 

(A. 3-8) 

O  0 

(■«"f2n+l‘=‘ f->  r>  ' 

(A.3-9) 

O  0 

r' i-> ' 

(A. 3-10) 

O  0 

(A. 3-11) 

0  o 

-  ‘  >  Pa.  Vi'r  )  ' 

O 

(A.3-12) 

^®^2n+l^^'^  ^  ■'  «  ®2n+1^2n+l^^  r  ^  ^ 

o 

(A. 3-13) 

i’  > 

(A.  3- lit) 

■"  n  ®2n+2^2n+2^^  r  ^  * 

o 

(A. 3-15) 

Fe4^(s,r  )  .  P2„(y  ^  )  K^(y  ,  )  , 

0  0 

(A. 3-16) 

^“^n+lfS'P  >  '  ‘  “ai*l  (S'  r  )  ‘‘anxl'i'  r  >' 

(A.3-17) 
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i_iL_  p  (y  £_)  K'  (y  ~) 
jf  ^2n+l  r  ^  2n+l^^  r  ^ 

o  o 


If  ®2n+2  r  ^  ^^+2^^  r  ^ 
o  o 


(A.  3-18) 
(A. 3-19) 


where  ^2n+l*  ®2n+l^  ^2n+2  connecting  factors  which  have  been 

2  22  22  2  22  22 
defined  in  section  A.l,  and  x  =  kfr^  -  Pr  ,  y  -  &  r  -  k-r^  in 

which  r^  is  the  red.ius  of  the  circle.  The  prime  on  the  modified  Mathieu 

function  represents  the  derivative  of  the  function  with  respect  to  |  , 

while  the  prime  on  the  Bessel  or  modified  Bessel  function  represents  the 

derivative  of  the  fianction  with  respect  to  its  argument. 
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J  “  %5o 


For  r  -  0 


•  2x 


f  ce|^(il,lr^i)dT,  =  2*  +  n  f  [4?^]' 

0 


(A.l^-6) 

(A.4-7) 

(A.4-8) 

(A.4-9) 

(A.4-10) 

(A. 4-11) 
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APPENDIX  B 

Tabulation  of  Integrals  Involving  Mathleu  and  Modified  Mathlfeu  Functions 

Integrals  resulting  from  the  attenuation  constant  and  Q  factor 
calculations  are  tabulated  In  this  appendix.  Integrals,  are  Integrated 
analytically  whenever  possible,  and  the  results  are  given.  The  Inte¬ 
grals  are  divided  Into  tws  categories;  those  Involving  the  Mathleu 
functions  are  called  the  angle  Integrals  since  they  correspond  to  the 
trigonometric  Integrals  of  the  circular  guide,  and  those  Involving  the 
modified  Mathleu  fvinctlons  are  called  the  radial  Integrals  since  they 
correspond  to  the  Bessel  Integrals  of  the  circular  guide. 

B.l  Angle  Integrals  Involving  Mathleu  Functions 

The  definite  integrals  involving  Mathleu  functions  can  usually  be 
integrated  analytically.  They  are 

r  =  0 

ce^(^)dT,  =  X  r  ^ 

r  »  u 

se  (i,)dTi  =  K 

r  =  0 

se'2(T,)dij=  K  f  B^^^(2r+1)2 

ce'_(n)se3^(,)dTi=.-  «  f  4rll 

r  =  0  . 

cej^(il).sej^(il)di|  =  -  I^ 


-192- 


1,  ./  se=(,)  =08  2,1,  .  K  [-  i  f  ^  B^ij 

I72  =  j  =e^(ti)  cos  2iidi| 

0 

l[  =  r  ce*'^(tj)dTi  =  «  E  (2r+l)^ 

r  ss  0 

!■=/'  ==.=(, )d,.. 

S  “/  =  *  f  ^  B*^3)^(2r+l)2 

0  •  .  » 


^  ce*^(Ti)dTj  =  n  5]  ® 


0 

2k 


r  =0 


♦(3)‘ 

2r+l 


0 

2k 


^  =[  ce*'(Ti)  ce*'(T))dT,  =  -  K  B|^j|(2r+1) 


^6  =  J  «e5(n)dn  =  -  «  Z 

•L  r=  0 


(3) 

2r+l 


0 

2k 


I.  .  r  8e'*2(,)d,  .  ,  2 

■i,  r  =  0 


00  2 

^9‘f  “r^fD'i'i  - « 

^io'/ 


/ 
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2k 

^11  =  /  8e5(n)<iil  ‘  ^  ^arii  A^+^^r+l) 


0 

2n 


’^ia  ■  /  ”  E 


0 

se*(Ti)dTj  = 

0 
2k 


r  =0 


00 


2r+l  2r+l 


K  V  B*^^j(2r+1) 

2r+l  2r+l  ■ 


III,  =  f  ce*(Tl)  8e*(Ti)dTj  =  - 

J 


0 

2i 


I.' 


r 


r  =  0 


13 


00 


2r+l .  2r+l 


15  "  I  sej(T|)dTj 


16 


17 


■f-i 


■f~i 


'(.T|)  ce*(Ti)d^ 


'(ii)  8e*(Tj)dij  = 


r  =  0 


■  ^15 


00 


0 

2k 


1£q  =  J  ce*(7i)  se*'(ii)dT}  =  - 
0 

=  J  ce^'(Tl)  8e*(n)dii  = 

■o' 

,^20  "  f  ce*(Ti)  8e^t(^)dTj  = 

0 


r  =  0 


^i7 


GO 


n  y  B*^^J(2r+l) 

2r+l  2r+l'  ' 

r=  0 


^19 


2k 

=  j'  8e*^(Ti) 
•0 


cos  2Tjdn  =  -  K 


-  +  y  A*^^^  A*^^^ 

2  ^  ^  ^  2r+l  *2 


r  =0 


2r+3 


2n 

r  2 

=  j  se*  (Tj)cob  2T]dT)  = 


,1  ^  .E„ 


r  =  0 


2r+3 


4 
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I se*(,).e.(,)=CB  2,1,  .  I f  *  A^3);^•(l), 


0 

2« 


;=J  <=*f(>i) 


0 

2tr 


,  =  J  ce*^(n) 


COB  2TldTJ  =  K 


COS  2TjdTJ  =  It 


1  ?  .•(!)  B»(l) 

'®-  ■  t„®2rtl®2r+3 


2  “1 


0 

¥ 


r  =  0 


.2  oo 


1  -o^CS)  ^  -ni 

-B*  ®2r+l.®2r+l 


b*(3)  ■n*(3) 


2  “1 


r  =  0 


r  =Jce*3(„)ce*(t,)coB  2T^t,  = -|  f 

o  r  =u 


The  prime  slgalfles  the  derivative  of  the  function  with  respect  to  t)  . 

2 

The  *  Indicates  that  the  argtunent  x  Qf  the  Mathleu  fmctlons  or  the 

2222 

expansion  coefficients  is  (p  -  cosh  1^/4. 

B.2  Reidlal  Integrals  Involving  Modified  Mathleu  Functions 

All  these  radial  Integrals  are  Integrated  numerically  by  Sliqpson's 
Rule  (^O).  The  normalized  dimensionless  variables 


^1  = 


xe 


5 


2  cosh  § 


^2  = 


ye- 


2  cosh  I 


2  2  2  2  2  2  2  2  2  2 
where  x  =  q  cosh  ^  y  =  q  cosh  i^CP  "  1^^)  are  used.  , 

■  •  '  •  •  r  ’  r 

As  the  ellipse  approaches  a  circle,  z^^-^x  —  Euid  ^2  ~  ^  ~ 
limits  of  Integration  instead  of  being  from  0  to  6  and  |  to  00 
will  be  from  - - -  to  .  and  t — ; —  to  00  respec- 


2  cosh  I  2  cosh  | 

®o  *0 


2  cosh  g 


tlvely.  Assuming 


a,  = 


1  2  cosh  I  ^  ^2  2  cosh  |  ^  "^1  2  cosh  f 

O  O  *0 


xe 


b,  = 


_  ye 


and  b»  =  00  , 


we  have  the  following  radial  Integrals 
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r‘^  2 

^8  “J 

®2.  ■; 

^ia‘J  ^C^h'>\^h 

®1  :  ■ 

r  p  dz,  . 

^u'i  ®^^^i>  -zr 

‘i-  ,  .  - 

“i<^i>  s=i(zi)dzj^ 

“i 

®2 

^13*/  Cej^(*i)  Se^(z^)dz^ 


r  o  r  ^  cosh  I  z, 

>.  'j  T~Tr' 

N  a^  L  4  cosh  6_z. 

J  Ce^(zjj^) 


2.  2  2 


-  X 


dz. 


HN  a. 


[  4  cosh^l^z^  -  dZi 


4  cosh'll  z-x 
o  1 


dz. 


r 

[^  =  J 

\  ^ 


r  ' 

1^2  -J  F=4‘ 


(z  )z  dz 
'  2^  2  2 


dz. 


(i,  "  /  P«1:?(0  — 

23  3V2'  zj 
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=24'/ 


2  2 

.  ^^2^ *2^*2 


^25 


4‘ 


dz. 
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26 


r 

=  J  Fek^ 


dZj, 

(z^)  Fek'Ckj)  — 


r  2  ^^2 

^27-]  ®'^i<'2>-t: 


28 


-  r  0e4= 


.v( 
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3'  2^  Zg 
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/■ 
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u 
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U 
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33  3v 
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fz 

Pekj[(z2)  Gek^(z2)dz, 
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Pek^Ug)  Gek£(z2)dZ2 
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Pek^CZg)  Gek^(z2)dz2 


\  Fek^Czg^  Gek'(z2)dz2 
^2 

I"  Pek^(z2)  Gek^(z2)dz2 


o, 

r 


Fek^iz^)  Gek^(z2)dz2 
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-2 

Gek^ 
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p  P  P 
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2  '  '  ■ 
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P  P  P 
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2 
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1 


■/ 
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\  ^14  cosh'll  Zpy  .  2 

^2.  0  2' 


r 

ie-l 


(Zg)  Fek^(z2) 


1.  w2.  2  ,2 
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4  cosh  I  z«y 
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The  prime  on  the  modified  Mathleu  functloh  indicates  the  derivative  of 
the  function  vith  respect  to  its  argument. 
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